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The  description  of  human  pilot  dynamic  characteristics  in 
mathematical  terms  compatible  with  flight  control  engineering 
practice  is  an  essential  prerequisite  to  the  analytical  treat¬ 
ment  of  manual  vehicular  control  systems.  The  enormously 
adaptive  nature  of  the  human  pilot  makes  such  a  description 
exceedingly  difficult  to  obtain,  although  a  quasi-linear  model 
with  parameters  which  vary  with  the  system  task  variables  has 
been  successfully  applied  to  many  flight  situations.  The 
primary  purposes  of  the  experimental  aeries  reported  are  the 
validation  of  an  existing  quasi-linear  pilot  model,  and  the 
extension  of  this  model  in  accuracy  and  detail. 

To  this  end  the  influences  of  controlled-element  dynamics 
and  system  forcing  function  characteristics  on  the  pilot's 
dynamic  characteristics  are  investigated  using  a  five-stage 
process:  (l)  Pre-experiment  analyses  are  conducted  using  the 
existing  model  to  predict  the  outcome  of  experiments  which  are 
especially  contrived  to  exercise  the  model  to  its  limits ; 

(2)  controlled-element  dynamics  which  are  both  crucial  task 
variables  (per  the  pre-experiment  analyses)  and  idealizations 
of  aircraft,  booster,  and  space  vehicle  dynamics  are  delineated 

(3)  describing  function  and  remnant  measurements  are  taken  in 
an  extensive  experimental  program  involving  eight  different 
controlled  element  forms  and  three  forcing  functions; 

(it-)  analytical  abstractions  of  the  data  are  made  by  curve- 
fitting  procedures;  (5)  variations  in  the  pilot's  oh,,  . ucter- 
istics  due  to  controlled  element  and/or  forcing  function 
changes  are  described  in  terms  of  the  parameters  of  the  curve 
fits.  The  outcome  of  this  process  is  a  substantially  refined 
and  extended  adaptive  and  optimalizing  model  of  human  pilot 
dynamic  characteristics.  Models  corresponding  to  three  levels 
of  precision  and  complexity  are  developed,  the  Beveral  aspects 
of  pilots '  adaptation  to  controlled  element  and  forcing  func¬ 
tion  changes  are  detailed,  the  selective  variability  nature  of 
human  pilot  dynamics  is  presented,  key  remnant  sources  are  dis¬ 
covered,  and  many  other  aspects  of  human  pilot  dynamics  are 
treated  with  a  combined  experimental-analytical  approach. 
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The  investigations  reported  here  are  an  element  in  a  U.  S.  Air  Force  | 
research  prottmm  to  explore  and  exploit  the  concept  that  vehicle  dynami  c  i 
handling  qualities  are,  to  a  large  extent,  dependent  on  the  action  of  /  * 
the  pilot  as  a  control  element  in  the  pilot-vehicle  closed-loop  system,  | 
In  thin  concept,  control  analytic  tcchniquca  and  pilot  dynamic  response  - 
models  are  used  in  the  study  and  optimization  of  man-vehicle  systems.  1; 
Such  procedures  promise  to  greatly  enhance  the  processes  Involved  in  the  $ 
design  of  manned  vehicles.  * 


This  report  documents  an  analytical  and  experimental  investigation 
of  human  pilot  dynamics  accomplished  under  Contract  AF  33(61 6) -7501 , 

Project  No.  8219,  TaBk  No.  821905,  sponsored  by  the  Flight  Control  Divi¬ 
sion  of  the  Air  Force  Flight  Dynamics  Laboratory.  The  researcn  was  per¬ 
formed  by  Systems  Technology,  Inc.,  at  both  its  Hawthorne,  California, 
and  Princeton,  New  Jersey,  offices,  and,  under  subcontract,  by  The 
Franklin  Institute  Laboratories  for  Research  and  Development,  Philadelphia, 
Pennsylvania.  The  project  principal  investigators  were  D.  T.  McRuer  and 
D.  Graham,  of  STI,  and  E.  S.  Krendel,  of  FIL.  The  Flight  Control  Division 
project  engineer  through  most  of  the  program  was  R.  J.  Wasicko,  succeeded 
in  the  last  phases  by  R.  H.  Smith  and  P.  E.  Pietrzak. 

As  in  most  team  efforts,  many  have  contributed  to  the  results  reported 
here.  The  major  contributors,  who  all  participated  in  the  analytical, 
planning,  experimental  execution,  data  interpretation,  and  reporting 
phases,  are  listed  aB  authors.  An  indispensable  portion  of  the  program, 
the  design  and  development  of  analysis  apparatus,  waB  accomplished  meticu¬ 
lously  by  R.  A.  Peters  and  K.  A.  Ferrick  of  STI.  Important  contributions/ 
were  alGo  made  by  R.  J.  Wasicko  in  experimental  planning,  and  by  R.  E. 
Magdaleno  of  STI  in  the  interpretation  of  the  resultB.  Special  mention 
and  thankB  are  due:  R.  P.  Harper,  Jr.,  of  Cornell  Aeronautical  Laboratory; 
G.  E.  Cooper,  of  Ames  Research  Center,  NASA;  and  Lt.  Comdrs.  M.  Johnson  and 
T.  Kastner,  Capt.  B.  Baker,  and  Lts.  G.  Augustine,  F.  Hoerner,  and  J.  Tibbs, 
of  the  Naval  Air  Test  Center,  Patuxent  River,  Maryland,  for  their  interest 
and  assistance  as  subjects.  The  authors  would  also  like  to  thank  their 
co-workers  Diana  Fackenthal,  S.  H.  Greene,  and  M.  M.  Solow,  of  FIL,  for 
their  contributions  to  running  the  experiments,  reducing  the  data  and 
assisting  in  its  analysis;  and  H.  B.  Grudberg,  A.  V.  Fhatak,  and  D.  B. 
McElwain,  of  STI,  for  their  assistance  in  pre-experiment  predictions, 
and/or  postexperiment  data  analysis  and  interpretation.  Acknowledgment  is 
due  to  Boll  Beranek  and  Newman,  Inc.,  of  Cambridge,  Massachusetts,  for  the 
amplitude  distribution  data  processing  and.  some  goodness  of  fit  analyses. 
Finally,  the  report  has  been  substantially  improved  by  the  Incorporation 
of  many  suggestions  due  to  the  extremely  careful  review  by  R.  0.  Anderson, 

P.  E.  Pietrzak,  and  R.  J.  Woodcock  of  FDCC. 

This  technical  report  has  been  reviewed  and  In  appr^' 


W.  A.  SLOAN, J»{ 

Colonel,  USAF 

Chief,  Flight  Control  Division 
AF  Flight  Dynamics  Laboratory 
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IlfRROICJ  OTIOS 

A.  ntOJlOI  lAQUKOUB  A*  FURKM 

The  effective  use  of  manned  flight  vehicles  has  always  required  a 
satisfactory  match  of  vehicle  characteristics  (which  include  vehicle 
dynamics,  control  manipulator,  display,  ete.)  with  the  human  pilot's 
characteristics  as  a  flight  controller.  An  agreeable  mating  is  not 
inherent  in  the  design  process,  and  the  provision  of  proper  vehicle 
handling  qualities  has  often  posed  serious  problems  which  the  vehicle 
system  designer  must  solve. 

Classically,  handling  qualities  concepts  were  based  on  engineering 
knowledge  of  vehicle  characteristics,  leavened  by  pilot  opinion  ratings. 
The  opinion  ratings  were  subjective  expressions  of  the  over-all  suita¬ 
bility  of  the  manual  control  system  consisting  of  the  pilot  and  the 
vehicle.  A  convenient  way  to  express  the  handling  qualities  was  as 
catalogs  of  vehicle  dynamic  parameters  given  as  functions  of  pilot  rat¬ 
ings.  In  Bpite  of  their  reliance  on  the  subtleties  of  subjective  pilot 
ratings,  such  catalogs  of  handling  qualities  characteristics  have  been 
evolved  in  the  past,  and  will  continue  to  be  expanded  in  the  future. 

But,  in  a  fundamental  Bense,  these  catalogs  are  only  reports  of  specific 
test,  results  —  they  fall  to  adequately  explain  the  mutual  interactions 
between  the  pilot  and  the  vehicle,  and  they  are  difficult,  if  not  impos¬ 
sible,  to  extrapolate  to  new  slbuatlonB  and  novel  vehicle  characteristics. 

To  achieve  understanding  and  the  capability  to  extrapolate  requires 
a  mathematical  theory  which  can  be  used  to  explain  old  findings  and  to 
predict  new  ones.  For  handling  qualities  a  theory  of  this  kind  has  been 
in  the  process  of  construction,  refinement,  and  successful  application 
for  some  years.  It  is  based  on  the  methods  of  control  engineering,  and 
treats  the  pilot-vehicle  system  as  a  closed-loop  (in  general,  a  multi¬ 
loop)  entity.  The  sine  qua  non  of  the  theory  is  a  model  of  pilot  dynamic 
characteristics  in  a  form  suitable  for  application  using  relatively 
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conventional  control  engineering  techniques.  Further,  the  applications 
to  which  the  theory  can,  with  confidence,  be  employed  is  limited  funda¬ 
mentally  by  the  level  of  pilot-model  knowledge. 

An  adequate  description  of  a  pilot's  dynamio  response  characteristics 
is  not  easily  obtained  because  of  the  pilot's  inherent  adaptability  and 
capacity  for  learning.  Nevertheless,  suitable  experimental  techniques 
have  been  devised  and  employed  in  the  past  to  provide  a  limited  amount 
of  data.  The  quality  of  some  of  these  data,  however,  has  been  open  to 
question  because  of  real  or  imagined  deficiencies  and  uncertainties  In 
the  experimental  task  variables  and  in  the  necessary  analysis  equipment 
such  as  spectral  and  cross- spectral  analyzers.  In  many  cases  the 
analyzer  calibrations  were  insufficient  to  establish  fully  the  accuracy 
of  the  data  reduction  methods  (Ref.  51).  Furthermore,  most  experimenters 
had  not  made  certain  measurements  which.  It  turned  out,  were  of  crucial 
Importance  in  the  context  of  pilob-vehicle  control  system  analysis.  In 
spite  of  suoh  imperfections,  the  collation  of  all  the  existing  data, 
expanded  by  ultraconservative  extrapolations  based  only  on  the  limited 
high- reliability  data,  yielded  a  data  base  from  which  a  serviceable,  but 
inconplete,  mathematical  model  was  evolved  (Refs.  28,  34  -  36).  The 
extensive  use  of  this  pilot  dynamics  model  in  handling  qualities  and 
pilot-vehicle  system  analysis  (Refs,  1-4,  8,  11,  12,  17,  27,  33-36,  43, 
*+9,  55..  56,  and  58)  heightened  the  desire  for,  and  increased  the  potential 
importance  of,  a  more  complete  understanding  of  the  mathematically  describ- 
able  aspects  of  human  dynamics  in  vehicle  control  ByBtems.  The  nature  of 
desired  improvements  in  the  model  was  fundamentally  one  of  degree  rather 
than  kind,  i.e.,  increased  scope  and  precision.  Such  an  expanded  view 
could  not  be  evolved  from  existing  data,  which  had  been  used  in  the  con¬ 
struction  of  the  model j  instead,  new  data  were  required.  Fortunately,  by 
i960  enough  effort  had  been  devoted  to  model  building  and  application  to 
give  a  more  definitive  notion  of  how  the  data  were  to  be  used.  This, 
coupled  with  a  very  much  better  appreciation  for  the  data  reduction  prob¬ 
lem,  made  possible  the  planning  and  initial  execution  of  a  research  program 
to  meet  the  real  needs. 

Work  on  the  project  was  initiated  1  July  i960,  and  almost  two  years 
were  spent  in  the  design,  construction,  and  calibration  of  the  data 


reduction  machines .  The  results  of  this  effort  have  been  reported 
elsewhere  (Refs.  9>  10,  1 6,  44-46).  Tests  with  human  operators  were 
initiated  in  May  1962.  It  is  primarily  to  the  task  of  recording  and 
Interpreting  the  data  from  some  of  these  experiments  that  this  report 
is  addressed. 

The  primary  results  desired  from  the  program  from  a  vehicle  oontrol 
standpoint  are,  of  course,  better  pilot  models  to  use  in  handling  quali¬ 
ties  and  manual  control  system  analyses.  Although  not  treated  here,  the 
applications  of  these  models  are  expected  to  be  far-reaching  in  the 
future.  One  reason  for  this  anticipation  are  the  extensive  uses  of  the 
much  less  precise  and  more  vague  circa  i960  model.  For  example,  in  the 
references  cited  above  this  model  has  been  used  to: 


1 .  Estimate  human  pilot  and  over-all  pilot-vehicle 
system  dynamic  response,  stability,  and  average 
performance. 

2.  Determine  barely  controllable  vehicle  dynamics  and 
controllability  boundaries. 

3.  Delineate  those  features  of  the  vehicle  dynamics 
which  are  most  likely  to  affect  the  vehicle 
handling  qualities . 

4.  Indicate  the  type  of  additional  system  equalization 
(to  be  achieved  via  the  display,  the  manipulator, 
or  by  vehicle  modifications)  desirable  to  achieve 
better  pilot  control — as  well  as  the  effects  on 
the  pilot  characteristics  of  such  modifications. 

5.  Find  the  maximum  forcing  function  bandwidth  com¬ 
patible  with  reasonable  control  action  on  the  part 
of  the  pilot. 


The  new  models  developed  here  are  Intended  to  be  used  for  the  same  sorts 
of  things,  but  with  far  greater  confidence  and  considerably  better  pre¬ 
cision.  The  refined  models  are  expected  to  be  very  useful  in  other  ways 
also.  For  example,  the  characteristics  of  the  human  pilot's  "actuator" 
and  "sensor"  dynamics,  which  were  previously  lumped  into  a  mid-frequency 
approximation  to  lower  and  higher  frequency  effects,  are  distinctly- 
reflected  in  the  new  data  and  models.  This  new  knowledge  should  have 
significant  impact  both  on  the  content  and  nature  of  the  information 
displayed  to  the  pilot  and  on  the  design  of  the  manipulative  devices 
with  which  the  pilot  imparts  his  desires  to  the  vehicle. 
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OUTLXSt  OF  Sal  MOTORS! 

As  noted  above ,  the  major  purpose  of  the  Gxpoi  imental  program  lu  be 
discussed  here  was  to  provide  data  for  the  extension,  and  validation  of  the 
existing  pilot  model.  Therefore  a  necessary  preliminary  is  to  consider  the 
status  of  the  mathematical  model  of  the  human  operator  at  the  time  the 
experiments  were  planned.  This  is  accomplished  In  the  next  chapter,  which 
comprises  a  statement  of  the  best  that  was  known  about  human  operator 
behavior  in  compensatory  tracking  in  i960  together  with  a  discussion  of  ' 
those  areas  in  which  specific  knowledge  oi’  operator  behavior  and  its 
measurement  was  either  totally  lacking  or  substantially  incomplete. 

Previous  knowledge  of  human  operator  behavior  in  compensatory  tracking  is 
summarized  in  the  first  section-  of  Chapter  II  entitled  "The  Analytical- 
Verbal  Model,"  while  the  discussion  of  what  it  appeared  necessary  to  find 
out  follows  in  a  treatment  entitled  "Objects  of  the  Experiments." 

A  major  part  of  the  plan  for  the  model  validation  aspects  was  to  make 
analytical  predictions  of  human  operator  performance  in  compensatory 
tracking  and  then  to  observe  whether  or  not  these  predictions  could  be 
confirmed  with  experimental  results.  To  this  end  the  analytical-^verbal 
model,  reviewed  in  Chapter  II,  la  bravely  put  to  use  in  Chapter  III. 

Although  it  was  recognized  that  additional  data  were  clearly  required  to 
substantiate  some  of  the  conjectures  on  which  the  model  was  founded, 
predictions  critical  to  model  validity  were  made  which  could  later  be 
compared  with  experimental  results.  These  pre- experiment  analyses  also 
provide  excellent  concrete  examples  of  the  techniques  involved  in 
applying  the  human  pilot  models. 

The  really  essential  portion  of  this  program  is  experimental,  for  the 
potential  of  all  existing  data  as  sources  for  model  buildl””  and  elaboration 
had  been  exhausted.  Past  experimental  efforts  have  sell  ,v  '  ,  irthy 

examples  of  precision.  In  the  new  experimental  serif  ■  .  -  impor¬ 
tant  desired  feature  was  the  provision  of  supplement  ,,  .•  ■„  ,ion 

techniques  and  methodc  which  would  maximize  flexibility  ^lze  the 

chances  of  experimental  error.  The  apparatus  web  deeig  make  these 

objectives  feasible,  but  changing  feasibility  to  actual.  requires  excep¬ 
tional  experimental  and  data  reduction  procedures  whir  almost  amount  to 
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trade  secrets.  Some  of  these  are  described  in  Chapter  IV,  together  with 
a  description  of  the  experimental  and  analytical  apparatus. 

The  experimental  objectives  of  Chapter  II  plus  the  pre-experiment 
i  analyses  of  Chapter  III  lead  directly  to  a  desired  experimental  program 

plan.  An  outline  of  this  program  is  presented  in  the  first  section  of 
Chapter  V. 

Chapters  V  and  VI  present  the  experimental  data — deBoribing  funotiona 
in  Chapter  V  and  remnants  in  Chapter  VI.  The  data  are  aggregated  in  vari¬ 
ous  ways  to  illustrate  the  groBS  and  detailed  effects  of  changes  in  the 
forcing  function  and  controlled  element  task  variables.  Other  aggregations 
are  used  to  illustrate  the  effects  of  intra-  and  intersubject  variations. 
The  describing  function  chapter  also  includes  statistical  treatments  of 
certain  major  conclusions.  A  special  feature  of  the  remnant  discussion  is 
a  presentation  of  data  which  indicate  the  likely  major  remnant  source,  as 
veil  as  data  tending  to  indicate  what  the  remnant  is  not. 

Chapter  VII  ie  devoted  to  detailed  analyses  of  the  data  and  to  the 
development  of  interpretations  and  rationalizations.  In  extending  the 
existing  analytical-verbal  model  it  is  concluded  that,  by  and  large, 
the  hypotheses  and  extrapolations  made  from  the  limited  data  previously 
available  were  generally  correct,  and  reasonably  explained.  Further, 
the  updated  models  developed  here,  which  subsume  the  old  model  and  con¬ 
form  to  the  new  data  as  well,  answer  many  other  questions  concerning 
human  behavior.  These  models  are  of  three  levelB  of  precision.  The 
first  is  relatively  crude,  and  ie  intended  for  use  in  the  region  of 
crossover  only.  The  second  is  sufficiently  precise  to  be  suitable  for 
most  handling  qualities  analyses,  which  primarily  emphasize  vehicle 
dynamics,  including  those  involving  statically  unstable  vehicles. 

Finally,  the  third  is  a  precision  model  which  is  capable  of  representing 
the  high  and  very  low  frequency  actuator  and  sensor  dynamic  characteris¬ 
tics  of  the  human  pilot.  Thus,  the  three  models  provide  a  range  of  com¬ 
plexity  and  utility  which  is  analogous  to  the  several  degrees  of  model 
complexity  used  in  ordinary  automatic  pilot  design. 

Finally,  Chapter  VIII  summarizes  the  general  conclusions  and  findings 
of  the  study. 
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a.  A  QUIDS  VQR  THE  MEADE* 


Because  this  le  a  very  Ions  report,  with  many  involutions,  it  is 
appropriate  to  give  some  words  of  advice  for  readers  with  various 
interests.  For  the  casual  reader,  this  introduction  and  the  general 
summary  of  Chapter  VIII  give  the  gist  of  the  effort.  If  just  a  bit 
more  is  desired,  add  Chapter  II. 

The  applications -oriented  engineer  should  stafct  with  Chapter  II, 
Section  A,  for  a  review  of  the  existing  model,  followed  by  Chapter  V, 
Section  A  for  the  experimental  plan  outline  and  Section  D  for  the  grand- 
average  describing  function  data,  and  then  on  to  Chapter  VIII  for  the 
general  summary  of  results.  He  will  then  probably  wish  to  absorb  more 
details  on  the  new  models,  as  developed  in  Chapter  VII,  and  may  deBire 
to  examine  the  pre-experimental  analyses  of  Chapter  III  aB  concrete 
examples  of  applications. 

Those  who  are  interested  primarily  in  data  can  proceed  directly  to 

Chapters  V  and  VI,  although  the  experimentalist  will  also  wish  to  cover 

Chapter  IV .  , 

I 

The  model  builders  can  turn  directly  to  Chapter  VII,  and,  if  they 
have  models  to  test  against  the  realities  of  data,  to  Chapters  V  and  VI. 

Finally,  for  the  reader  who  is  interested  in  the  entire  effort  the 
way  is  directly  through  the  report  as  laid  out,  although  many  of  the 
details  can  be  sloughed  over  lightly  on  a  first  reading.  These  diligent 
souls  have  both  the  authors'  sympathy  and  blessings! 
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1 .  Oensrml  fetus*  of  tfas  Nodal 


The  primary  objective  of  most  of  the  past  experimental  and  analytical 
programs  to  develop  mathematical  descriptions  for  pilot  response  charac¬ 
teristics  has  been  to  achieve  reasonable  descriptions  of  the  pilot  as  a 
component  In  engineering  systems.  Major  efforts  in  model  building  have 
thus  been  placed  on  the  evolution  of  models  which  can  predict  pilot 
dynamic  response,  characteristics  of  engineering  significance,  but  which 
are  otherwise  of  minlmiim  analytical  complexity*  Such  models  are  concep¬ 
tual  descriptions  of  the  human.  They  are  intended  to  exhibit  analogous 
oause-and-effeet  behavior  rather  than  to  be  analogs  in  any  structural 
sense.  The  models  are  valid  to  the  extent  that  their  behavioral  proper¬ 
ties  resemble  the  performance  of  the  actual  human  operator.  They  gain  in 
scope  and  descriptive  power  if  oartaln  of  their  features  oan  also  be 
identified  structurally,  although  they  cannot  be  rejected  because  of  any 
failure  to  satisfy  this  desirable  quality. 


As  a  control  component  the  human  exhibits  a  bewildering  variety  of 
nonlinear  and  time-varying  behavior.  Table  I  lists  some  of  these.  An 
appropriate  type  of  engineering  mathematical  description  for  nonlinear 
control  elements  of  this  nature  is  some  hind  of  quasi-linear  system. 
This  is  an  equivalent  system  in  which  the  relationships  between  some, 
but  not  necessarily  all,  pertinent  measures  of  system  Input  and  output 
signals  have  "linear- like"  features  for  fixed  input  conditions  in  spite 
of  the  presence  of  nonlinear  elements. 


The  quasi-llne&r  system  concept  originally  evolved  from  the 
observation  that  a  great  many  nonlinear  systems  have  responses  to 
specific  inputs  which  appear  similar  to  responses  of  equivalent  linear 
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TABLE  I 


HUMAN  OPERATOR  NOHUKSARITIES  ADD  NONSIATIONARITIES 


1 .  Adaptation  and  Learning 

The  adaptive  human  responds  to  changes  in  external  environment 
by  modifying  characteristics  so  as  to  improve  performance  over  that 
which  would  be  achieved  if  the  characteristics  used  in  the  original 
environment  were  maintained  in  the  new  one. 

The  learning  human,  after  operating  experience  in  a  given 
external  environment,  modifies  characteristics  to  achieve  better 
performance  than  previously  exhibited  in  the  given  environment. 

2.  Set  Changes 

"Bet"  characterizes  the  temporary  operating  points  or  baseline 
conditions  of  the  human  subsystems  involved  in  the  control  tasks. 
Changes  in  these  internal  conditions  facilitate  a  certain  more-or- 
less  specific  type  of  aotivity  or  response  in  the  adaptation  and 
learning  processes.  Set  changes  include: 

Internal  system  topography  changes,  l.s.,  uss  of  different  feed¬ 
back  and  feedforward  paths 

Variation  in  steady-state  muBcle  tension 

Variation  in  force  ranges 

Variation  in  indifference  threshold 

3.  Series  Nonlinearities 

Sensory  thresholds 

Maximum  force  and  displacement  limits 

4.  Fluctuations  in  Attention,  Motivation,  Etc.;  General  Drifting  of 
Characteristics 


system  to  the  same  specific  inputs.*  For  such  combi  nation?  of  specific 
inputs  usd  nonlinear  systems!  the  response  of  the  nonlinear  system  can  be 
divided  into  two  parts— one  component  which  corresponds  to  the  response  of 
the  equivalent  linear  element  driven  by  the  particular  input,  and  an  addi¬ 
tional  quantity  which  represents  the  difference  between  the  output  of  the 
actual  nonlinear  system  and  the  equivalent  linear  element.  This  second  compo¬ 
nent  is  called  the  "remnant"  because  it  is  left  over  from  the  portion  of  the 
system  response  representable  by  a  linear  element.  Quasi-llnear  equivalents  to 
the  nonlinear  system,  for  the  specific  input  of  Interest,  are  characterized 
mathematically  by  a  describing  function  (which  is  the  equivalent  linear 
element)  and  the  remnant.  An  essential  feature  is  that  the  quasi-linear 
system  has  a  response  to  the  input  in  question  identical  to  that  of  the 
original  nonlinear  system,  so  the  quaBi-linear  system  is  an  exact  cause- 
effect  representation  of  the  nonlinear  system  for  the  speoiflc  kinds  of 
Inputs  and  responses  considered.  When  the  inputs  are  changed,  the  quasi- 
linear  model  also  changes.  If  the  device  were  such  that  its  quasi-linear 
system  remained  the  same  for  all  kinds  of  forcing  functions,  and  if,  further, 
the  remnant  were  zero,  then  the  system  would  have  a  constant-coefficient 
linear  nature. 

The  most  oommon  quasi-linear  system  element  in  engineering  use  is  the 
sinusoidal- input  describing  funotlon,  which  is  of  such  great  value  in 
stability  studies  of  nonlinear  servomechanisms.  Here  the  action  of  the 
describing  funotlon  on  a  sinusoidal  input  results  in  an  output  which  is 
the  fundamental  of  the  output  of  the  aoWal  nonlinear  system.  The  remnant, 
which  must  be  added  to  the  output  fundamental  to  achieve  equivalence  with 
the  nonlinear  system,  Is  made  up  of  s 11  the  higher  harmonics  resulting  from 
the  passage  of  the  sinusoid  through  the  nonlinearity.  Describing  functions 
aan  also  be  defined  for  transient  Inputs,  such  as  step  functions,  and  for 
random,  inputs.  In  principle  the  systems  can  be  time-varying  as  well  as 
nonlinear.  Random-input  quasi-linear  systems  representing  the  human  oper¬ 
ator  for  certain  conditions  are  the  type  pertinent  to  the  operator  data 
treated  here. 

"Many  texts  on  nonlinear  control  theory  treat  aspects  of  quasi-linear 
systems.  Chapters  3-6  of  Ref.  21  are  especially  pertinent  as  background 
in  the  context  of  this  report. 
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6.  9»riw  of  Bslsttaj  Quael-IlfiMF  iystoa  Poioriting 
Funotlon  Flu*  Snout  Modal  for  Sl&cl*-Loop  Tasks 

In  tarma  of  Idealizations,  the  simplest  manual  control  system  Is  that 
shown  in  the  block  diagram  of  Fig.  1.  This  is  also  the  most  commonly 
occurring  system  structure  in  practical  manual  control  systems— ■  either 
as  a  total  system  by  itself  or  as  a  component  part  of  a  more  complex 
multiloop  system.  To  control  engineers,  this  is  a  single-loop  feedback 
system  (except  for  feedbacks  Internal  to  the  operator);  to  engineering 
psychologists,  it  is  a  compensatory  system.  The  important  single-loop 
features  are  the  solitary  stimulus  (the  error)  and  the  random-appearing 
nature  of  the  forcing  function  (the  system  input) .  If  either  of  these 
is  changed,  the  whole  complexion  of  the  task  also  changes.  For  instance, 
if  the  pilot  is  shown  the  system's  input  and  output  directly  (a  pursuit 
display)  instead  of  their  differences,  he  is  often  able  to  take  advantage 
of  the  additional  information  and  thereby  to  improve  over-all  system 
performance.  Also,  if  the  forcing  function  is  not  random-appearing,  but 
perhaps  periodic  over  relatively  short  time  Intervals,  the  pilot  can 
often  detect  and  anticipate  the  repetitive  or  deterministic  nature  of  . 
the  Input  and  adjust  his  response  accordingly.  Both  of  these  higher 
order  types  of  behavior  amount  to  the  presence  in  the  system  of  further 
signal  paths  and  a  more  complex  than  single-loop  structure. 

For  the  system  shown  in  Fig.  1  there  are  three  task  variables  that 
have  a  major  effect  on  the  pilot's  dynamics — the  forcing  function  char¬ 
acteristics,  the  controlled-element  dynamics,  and  the  manipulator.  Many 
other  factors  are  implicitly  involved.  These  include  operator-centered 
variables  such  as  training,  fatigue,  and  motivation,  and  external  environ¬ 
mental  characteristics  such  as  ambient  illumination  and  temperature. 
Ideally,  all  of  these  implicit  (or  "procedural")  variables  should  be 
taken  into  aocount,  and  someday  perhaps  they  will.*  But  for  the  present 

^Chapter  VII  of  Ref.  34  presents  a  preliminary  discussion  of  the 
effects  of  a  few  implicit  variables  on  pilot  response  measures.  Refer¬ 
ence  48  provides  an  excellent  detailed  example  of  pilot  describing  func¬ 
tion  changes  accompanying  ohanges  in  a  typical  environmental  variable — 
in  this  caBe,  the  pilot's  effective  "g"  field. 
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Figure  1  -  Single-loop  Manual  Control  System 


the  central  aim  la  to  explore  operator  dynamics  in  specialized  situations 
wherein  the  vast  majority  of  these  procedural  variables  are  held  constant. 
This  is  simple  enough  to  do  for  the  external  environmental  factors,  but  for 
the  operator-centered  factors  the  best  that  can  be  hoped  for  is  the  estab¬ 
lishment  of  reasonably  stable  levels  of  stationarity .  Stationarity  in  the 
over-all  experimental  situation  is  enhanced  by  confining  the  forcing  func¬ 
tions  to  signals  possessing  stationary  characteristics ,  and  by  using  highly 
trained  and  motivated  subjects  drawn  from- a  narrowly  limited  population  for 
which  high-grade  skill  in  manual  control  is  an  essential  feature. 

TuBtin  (Ref.  59)  first  noted,  in  a  formal  way,  that  operators  in  manual 
control  systems,  responding  to  random-appearing  visual  forcing  functions, 
exhibit  a  type  of  behavior  which  ie  analogous  to  the  behavior  of  equalizing 
elements  inserted  into  a  servo  system  to  improve  the  over-all  dynamic  per- 
formance.  Since  then,  a  number  of  measurements  of  human  response  to  visual 
inputs  have  been  made  in  situations  such  as  the  one  Illustrated  by  the  con¬ 
trol  system  block  diagram  of  Fig.  1 .  For  the  actual  measurement  situations 
the  human  being  Is  represented  by  his  quasi-linear  model,  that  is,  as  a 
describing  function  plus  a  remnant.  The  dynamics  of  the  display  and  other 
system  elements  are  lumped  into  a  "controlled  element,"  and  the  Bystem  forc¬ 
ing  funotlon  Is  modified  (if  necessary)  into  an  equivalent  forcing  function. 
The  equivalent  block  diagram  then  takes  the  form  shown  in  Fig.  2.  The  con¬ 
trol  loop  signals  are  represented  as  time  functions  and  their  Fourier  trans¬ 
forms,  e.g.,  e(t)  and  E(jcu),  and  also  as  powsr  spectra  or  power  spectral 
densities,  e.g.,  ®ae(ai) .  The  linear  constant-coefficient  controlled  element 
is  totally  defined  by  its  transfer  function,  Y0(jo>),  whereas  the  nonlinear 
time-varying  human  requires  the  describing  funotlon,  Yp,  and  the  remnant 
power  speotrum,  Qnna >  to  provide  an  adequate  dynamic  description  in  the 
sense  that  the  power  spectral  densities  of  the  signals  in  the  actual  and 
the  quael-llnear  equivalent  system  are  the  same. 

The  number  of  conditions  studied  by  the  principal  investigators  of 
human  operator  describing  functions,  through  the  year  I960*  are  summa¬ 
rized  in  Table  II.  The  table  is  organized  with  respect  to  task  variables. 
The  most  influential  of  these  turned  out  to  be  the  forcing  function  and 
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the  dynamic  a  of  the  controlled  element#  By  comparing  the  results  from 
all  these  experiments,  the  Influence  of  the  manipulator  was  shown  to  to 
unimportant  for  the  ranges  of  frequency  measured  and  manipulators  tested. 

The  dependence  of  the  human  operator  describing  funotlon  on  forolng 
funetlon  and  controlled  element  dynamics  has  tended.  In  the  actual 
experimental  situations,  to  be  obscured  by  the  limited  frequency  range 
(bandwidth)  and  the  run- to-run  variability  of  the  measurements.  Refer¬ 
ence  34,  however,  shows  how,  by  considering  averaged  data,  it  is  possible 
to  formulate  a  fairly  simple  analytical  describing  function  form  which 
can  be  adjusted  to  describe  the  main  features  of  human  behavior.  Although 
the  Hall  data  (Refs.  24  and  25)  were  not  available  at  the  time  Ref.  34 
was  written.  Ref.  36  represents  an  effort  to  bring  these  data  into  conso¬ 
nance  with  the  results  considered  earlier,  and  the  circa  1 960  statement 
of  the  model  below  takes  account  of  those  results.  Thus,  when  the 
\  describing  function  data  for  all  the  experiments  represented  in  Table  II 
''are  considered  as  a  whole,  they  serve  as  the  data  base  for  evolution  of 
a  'servo  model  describing  human  operation  and  adaptation  for  compensatory 
tracking  with  a  visually  presented,  random-appearing  forolng  function. 

This  model  is  the  key  element  in  a  dynamic  description  of  the  human  oper¬ 
ator's  capabilities  in  such  tasks  (the  other  element  being  the  remnant). 

It  characterizes  the  predominant  majority  of  all  the  experimental  results. 

The  model  comprises  two  elements: 

a.  A  generalized  describing  function  form 

b.  A  series  of  "adjustment  rules"  which  specify  how 
to  "set"  the  parameters  m  the  generalized 
describing  function  so  that  it  becomes  an 
approximate  model  of  human  behavior  for  the 
particular  situation  of  interest 


The  most  extensive  and  generalized  describing  function  form  for  one- 
and  two-dimessional  compensatory  control  tasks  developed  in  Ref.  34  is: 


(TlJcj+  1) 
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where 


Kp 

t 

(1^0}+ 1) 

(TxJcb+0 

[($+*&* ']<*»*') 

The  describing  function  is  written  in  terms  of  the  frequency  operator, 

Joo,  Instead  of  the  Laplace  transform  variable,  s,  to  emphasize  that  this 
describing  function  is  only  valid  .  the  frequency  domain  and  only  exists 
under  essentially  stationary  conditions.  For  instance  it  cannot  be  used 
to  compute  the  system  response  to  a  discrete  input,  e.g.,  a  step  response. 

The  Indifference  threshold  effect  was  derived  as  a  serial  member  In 
the  operator's  characteristics  in  the  Goodyear  studies  (Refs.  18  and  19), 
and  was  shown  (Ref.  J4-)  to  be  compatible  with  Elkind ' s  variable  amplitude 
resultB  (Ref.  13).  However  the  ratio  of  threshold  to  itB  input,  awji/ cep, 
is  quite  small  relative  to  one  for  input  signal  levels  conventionally 
used  in  tracking  tests,  so  the  describing  function  Kj  is  near  unity.  The 
Indifference  threshold  is,  therefore,  a  second-order  effect  that  can  be 
ignored  here  (although  it  can  be  of  importance  in  other  applications) . 

The  third-order  neuromuscular  system  description  Bhown  in  Eq  1  is 
based  primarily  on  high  frequency  data  available  from  step  funotion  inputs. 
None  of  the  describing  function  measurements  referred  to  in  Table  II  give 
a  direct  indication  of  the  complete  third-order  representation  because 
there  is  no  forcing  function  power  at  the  high  frequencies  about  ujj. 

The  describing  function  data  do,  however,  reflect  the  low  frequency 
effects  of  the  neuromuscular  system.  In  fact,  these  data  support  a 
first  approximation  to  the  neuromuscular  system  consisting  of  a  first- 
order  lag.  Thus  the  uBual  version  of  the  describing  funotion  model 


-  gain 

-  reaction  time  delay  (transport  lag) 

-  equalization  characteristic 

■  indifference  threshold  describing 
funotion  (1  -  y'S/rt  a^/oj  when 
a r/ap  «  1 ) 

-  neuromuscular  system  characteristic 
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develops!  In  Ref.  34  ties  only  e  first-order  neuromuscular  lag  term, 

(%)«  + 1  r1 ,  where  %  *  Tty  +  (2JN/<%) . 

With  the  simplification  discussed  above,  the  general  low  frequency 
describing  function  of  Eq  1  becomes  approximately 

Ip  ’  (Ii>  +  I  )(%)!■  + 1)  ^ 

While  Eq  2  can  be  shovm  to  have  a  very  considerable  range  of  validity 
for  a  variety  of  operators,  forcing  functions,  controlled  elements,  and 
manipulators,  only  the  form  of  the  describing  function  has  this  over-all 
validity.  Moat  of  the  parameters  in  the  describing  function  are  adjust¬ 
able  as  needed  to  make  the  system  output  follow  the  forcing  function  — 
i.e.,  the  parameters  as  adjusted  reflect  the  pilot's  efforts  to  make  the 
over-all  system  (including  himself)  stable  and  the  error  small. 

The  pure  time  delay  represented  by  the  e”^  term  is  due  to  sensor 
excitation  (the  retina  in  the  visual  case),  nerve  conduction,  computa¬ 
tional  lags,  and  other  data  processing  activities  in  the  central  nervous 
system.  It  is  closely  related  to,  but  not  identical  with,  certain  kinds 
of  classical  reaction  times,  t  is  currently  taken  to  be  a  constant 
because  it  appears  to  be  essentially  Invariant  with  forcing  function  and 
oontrolled  element  dynamlas  for  either  single  or  dual  random-appearing 
input  tasks.  However  both  inter-  and  intra-sub Jeot  t  variations  occur. 
Observed  t'b  run  as  low  as  about  0.1  sec  and  as  high  as  0.2  sac. 

The  neuromusoular  lag,  Tjj,  is  partially  adjustable  for  the  task. 

The  nature  of  the  adjustment  is  somewhat  obscure  due  to  the  leak  of  high 
frequency  data,  although  the  general  trend  is  a  monotonic  decrease  in  TN 
with  Increasing  forcing  function  bandwidth  (see  Table  13,  Ref.  34).  The 
observed  variation  of  Tjj  with  forcing  function  bandwidth  ranges  from 
less  than  0.1  seo  to  somewhat  greater  than  0.6  sec.  Because  the  details 
of  the  Tfl  variation  with  forcing  funotion  bandwidth  are  not  known,  this 
Important  variation  haB  often  been  ignored  in  applications  and  typioal 
values  of  Tjj  near  0.1  sec  have  been  used. 

The  equalizing  characteristics,  (T^Jm*  1 )/(TIJ<o  + 1 ),  coupled  with 
the  gain,  Kp,  are  the  major  elements  in  that  adaptive  capability  of  the 
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human  Which  allows  him  to  control  many  differing  dynemlo  devtnen.  ThMr 
function  is  the  modification  of  the  stimulus  signal  into  a  suitable 
nsuromusauler  aommand  which  is  properly  scaled  and  phased  for  Suitable 
over-all  man-raaohine  system  operation.  For  given  input  and  controlled 
element  characteristics,  the  form  of  the  equalizer  is  adapted  to  compen¬ 
sate  for  the  controlled  element  dynamics  and  the  pilot's  reaction  time 
delay. 

The  describing  function  adjustment  rules  are  not  Bimply  stated  since 
they  depend  intimately  on  interactions  of  the  elements  in  the  man-machine 
system.  In  general,  the  adjustments  can  artificially  be  divided  into  two 
categories  —  adaptation  and  optimalization.  Broadly  speaking,'  adaptation 
is  the  selection  by  the  operator  of  a  specific  form  (lag-lead,  lead- lag, 
pure  lead,  pure  lag,  or  pure  gain)  for  the  equalization  characteristics j 
and  optimalization  is  the  adjustment  of  the  parameters  of  the  selected 
form  to  satisfy  some  internally  generated  criteria.  The  result  of  the 
adaptation  process  is  fairly  well  understood,  since  the  form  selected  is 
one  compatible  with  good  low  frequency,  closed-loop  response  and  the 
absolute  stability  of  the  system.  The  internal  optimalizing  criteria 
are  not  known,  although  they  appear  to  be  generally  compatible  with  the 
minimization  of  the  rms  error  (Refs.  51  and  54). 

The  known  adjustment  rules  for  the  human  operator's  describing 
function,  in  decreasing  order  of  their  certainty,  can  be  summarized  as 
follows  (Refs.  54,  56,  and  57) s 

1 .  Stability:  The  human  adapts  the  form  of  his  equalizing 
characteristics  to  achieve  stable  control. 

8.  Tons  Saltation— Low  frequency :  The  human  adapts  the  form  of 
his  equalizing  characteristics  to  achieve  good  low  frequency 
closed-loop  system  response  to  the  forcing  function.  A  low  frequency 
lag,  Tp,  is  generated  when  both  of  the  following  conditions  apply : 

a.  The  lag  would  improve  the  system  low  frequency 
characteristics. 

b.  The  controlled  element  characteristics  are  such  that 
the  Introduction  of  the  low  frequency  lag  will  not 

result  in  destabilizing  effects  at  higher  frequencies  which  cannot  be 
overoome  by  a  single  first-order  lead,  Tp,,  of  somewhat  indefinite  but 
modest  size. 
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3.  fens  Saltation  —  Lttdt  After  good  low  frecuency  characteristics 
are  assured;  within  the  above  conditions;  lead  is  generated  when 
the  controlled  element  characteristics  together  with  the  reaction  tine 
delay  are  such  that  a  lead  term  would  be  essential  to  retain  or  improve 
high  frequency  system  stability* 

Parameter  Adjustment :  After  adaptation  of  the  equalizing  form; 
the  describing  function  parameters  are  adjusted  so  that: 

a.  Closed-loop  low  frequency  performance  in  operating  on 
the  forcing  function  ie  optimum  in  some  sense  analogous 

to  that  of  minimum  mean-squared  tracking  error. 

b.  System  phase  margin,  9^,  lies  somewhere  between  0°  and 
50°  when  the  forcing  function  bandwidth,  is  much  less 

than  the^ system  crossover  frequency,  and  from  50°  to  110°  when  is 
near  o^.  This  strong  effect  of  forcing  function  bandwidth  on  the  phase 
margin 'is  associated  with  the  variation  of  with  the  same  task  variable .** 

3.  Invarianoe  Properties 

a.  —  Kq  Independence :  After  initial  adjustment,  changes 
in  controlled  element  gain,  Ko,  are  offset  by  changes  in 

pilot  gain,  Kpj  i.e.,  system  crossover  frequency,  o^,  is  invariant  with  Kq. 

b.  %  —  Independence:  System  crossover  frequency  does  not 
depend  on  forcing  function  bandwidth  for  flu*  <  ti^,  .  (cdc0  is 

that  value  of  o\,  adopted  for  ojj.  «  <%.)  0 

0.  Wq  Regression:  When  nears  or  becomes  greater  than  a^0,  the 
crossover  frequency  regresses  to  values  much  lower  than  2^, 


■^"Forcing  function  bandwidth"  1b  a  vague  term  unless  the  forcing  func¬ 
tion  spectrum  is  rectangular.  For  other  spectral  shapes  an  effective 
rectangular  forcing  function  bandwidth  must  be  defined.  Several  possible 
baseB  for  this  exist,  but  recent  results  by  Elkind  (Ref.  Ik)  can  be  used 
to  support  the  selection  of  "effective  degrees  of  freedom"  as  the  basis 
for  defining  a  rectangular  bandwidth  equivalent  for  a  nonrectangular 
spectrum,  i.e., 


[fa  gjldaf 


fn  (®ii)£ 


den 


(3) 


The  phase  margin  adjustment  rule  noted  here  has  undergone  many 
changes  in  the  course  of  time.  Initially  (Ref.  3k)  a  range  of  0-30°  was 
suggested;  this  was  based  almost  entirely  on  extrapolation  of  low  fre¬ 
quency  data  to  crossover.  Subsequent  re-examination  of  some  of  the  Elkind 
data  having  forcing  function  power  in  the  crossover  region  led  to  exten¬ 
sion  of  the  upper  value  to  about  6o°.  The  Hall  data  almost  all  spanned 
the  crossover  region  (nee  Table  II) .  Phase  margins  extracted  from  these 
(Ref.  36)  were  highly  variable,  but  ranged  from  30°  to  110°.  The  criterion 
will  again  be  modified  in  this  reportl 
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It  turns  out  that  the  operator  describing  function  ewlaptwu  xor  a 
given  teak  is  very  similar  to  the  one  that  e  oontrol  engineer  would  select 
If  ha  vara  given  an  element  to  oontrol  together  with  e  controller  "black 
box/'  having  within  it  element*  making  up  the  describing  function  given 
by  Bq  2,  and  knobs  on  the  outside  for  the  adjustment  of  Tj,  Tj,,  and  Kp. 

Thus,  the  adaptation  of  equalization  form  covered  by  Adjustment  Rules  2 
and  3  may  be  foreign  to  the  reader  who  le  not  thoroughly  grounded  In 
control  theory  since  it  is  analogous  to  operations  which  hav<^  an  artiBtic 
flavor  even  in  conventional  control  system  synthesis.  Examples  often  help,  so 
Table  III  is  presented  to  illustrate  the  equalizer  form  taken  for  several 
simple  controlled  elements. 

TABLE  III 

EQUALIZER  FORMS  FOR  DIFFERENT  CONTROLLED  ELEMENTS 


WHEN  TEE 

CONTROLLED  ELEMENT 
TRANSFER  FUNCTION,  Y0,  IS 

THE  OPERATOR'S 
EQUALIZER  FORM 

ADAPTED  IS 

Kc 

Ja> 

Kb 

Kq 

(Jo)2 

KC 

(Jo)2  +  25%(Jo)  +  Ota 

Pure  gain,  Kp 

Lag-lead,  Tj  »  Tl 

Lead-lag,  T^  »  Tj 

|  Lead-lag  (if  %  «  S/t) 

(  lag-lead  (if  %  »  2/v) 

This,  then,  is  the  analytical-verbal  describing  function  model  of 
the  human  operator  which  existed  prior  to  the  ourrsnt  study.  The 
analytical  portion  of  tha  model  is  the  expression  of  Eq  2,  and  the 
verbal  portion  comprises  the  adjustment  rules  given  in  the  numbered 
statements  above. 

While  the  describing  function  is  the  critical  factor  in  determining 
operator-system  stability,  the  uneorrelated  —  in  a  linear  sense — portion 
of  the  operator's  output,  the  remnant,  oan  be  important  for  estimates  of 
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system  performance,  aiioh  •«  meen-erpierert  error,  This  Situation  ariSOS 

because  tbe  remnant,  represented  as  a  power  spectral  density,  ®nn(co) , 
is  one  of  two  additive  terms  in  the  expression  for  the  epeotral  density 
of  the  operator's  output. 


Considering  that  no  disturbances  are  present,  the  mean- squared  error 


35  "  hL  |f®| 


To(d®> 

+  YpYe(jcn) 


where 


^nncf®)4"3 


and 


E(jco) 

TO) 

*nn0  (“) 


~~Y  Y  J j to)'  *  tlie  errorAnPu,t  describing  function 
the  power  epeotral  density  of  the  input 

the  power  spectral  denBity  of  the  remnant  expressed 
as  an  "equivalent"  open-loop  input  applied  at  the 
pilot ' s  output 


The  first  term  in  Eq  4  derives  from  the  describing  function  portion  of 
the  pilot  model  operating  on  the  forcing  function.  At  the  operator* • 
output  this  component  can  be  represented  by  a  power  spectre  1  density 
|Ht2«ii(cn) ,  where  H  ■  Yp/(1  +  YpYc)  is  the  closed-loop  describing  func¬ 
tion  relating  pilot's  output  to  the  system  forcing  function.  The  power 
spectral  density  of  the  remnant  can  also  be  expressed  in  closed- loop  form 
as  ©nn  ■  |l/(l  +  YpYn) |£ftnnc.  Then  the  total  power  spectral  density  of 
the  operator's  output  is 

0Oo(®)  ■  I  H|  2©ii(oi)  +  ^(cd)  (5) 


The  ratio  of  the  linearly  correlated  pilot- output  power  to  the  total 
pilot-output  power  is  the  square  of  the  "correlation  coefficient,"  p: 


*nn 

$00 


(6) 


The  meaning  of  £  in  a  specific  instance  is  dependent  on  the  data 
analysis  apparatus  and  on  the  nature  of  the  system  forcing  function. 
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For  our  experiments  p  Is  found  using  an  analyzer  which  mechanizes 
spectral  and  cross-spectral  measurements  using  multiplications  and 
very  low  pass  filters  (Chapter  IV) .  If  the  forcing  function  is  a  sum 
of  sinusoids ,  9a  will  he  a  sum  of  delta  functions  (i.e. ,  a  series 
of  line  spectra  which  exist  only  at  the  frequencies  of  the  individual 
forcing  function  sinusoidal  components) .  Then,  in  general,  the  output, 
$QC,  will  be  a  sum  of  delta  functions  at  the  same  frequencies  as  those 
in  $££,  plus  delta  functions  at  other  discrete  frequencies  (if  non- 
linearities  or  constant  rate  sampling  are  present),  plus  a  continuous 
power  spectral  density  component  representing  random  fluctuation  in 
the  output*  At  the  frequencies  for  which  they  exist  the  delta  function 
components  will  generally  overpower  the  random  component,  and  the  p 
measured,  at  forcing  function  frequencies  will  generally  he  1 .0  unless  low- 
frequency  time  variations  in  H  result  in  additional  power  within  the 
measurement  filter  bandwidth.  In  fact,  p  will  be  1 .0  even  in  the  pres¬ 
ence  of  many  kinds  of  system  nonlinearities .  At  other  i  frequencies  p 
will  be  undefined  since  4>j_i  is  zero. 

For  forcing  functions  which  are  samples  of  random  processes,  the 
power  spectral  densities  in  Eq  6  are  all  continuous.  The  meaning  of  p 
for  this  case  is  quite  different;  i  ts  value  relative  to  1 .0  becomes 
primarily  a  measure  of  the  relative  importance  of  thu  remnant*  Near- 
unity  p  values  indicate  a  linear  constant-coefficient  system,  whereas 
lesBer  values  imply  nonlinearity  and/or  nonstatiorarity  and/or  "noise" 
injection  (Ref.  34).  Thus  p  for  the  random-input  case  is  not  as  dis¬ 
criminating  a  measure  as  when  the  Input  is  made  up  of  simple  sinusoidB. 

All  past  human  operator  data  for  which  p  values  exist  used  forcing 
functions  whioh  were  samples  of  random  processes  (Refs.  24,  29,  34, 
and  JO),  or  very  many,  closely  spaced  in  frequency,  sinusoids  (Ref.  13) 
which  were  not  separable  in  the  analysis  technique  used.  Therefore 
past  p  data  fall  into  the  random-input  category.  Much  of  the  existing 
data  show  correlation  coefficients  near  1 .0,  although  considerably  lower 
values  were  not  uncommon.  In  general,  the  larger  the  correlation  coeffi¬ 
cient,  the  smaller  was  the  remnant  and  the  mean-squared  tracking  errors. 
The  smallest  observed  remnants  ocourred  in  connection  with  controlled 
elements  whioh  had  the  least  energy  storage,  i.e.,  YQ  ■  Kg.  Here 
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the  form  of  the  operator's  equaliser  characteristic  is  either  a  low 
frequency  lag  or  a  pure  gain.  Finally,  the  least  remnant  i s  associated 
with  forcing  functions  having  the  least  bandwidth  or  high  frequency 
components. 

Unfortunately,  mechanisms  for  the  description  of  the  remnant  are  not 
nearly  as  well  "understood"  in  an  equivalent  mathematical  model  form  as 
those  for  the  describing  function.  Observed  remnants  in  one-dimensional 
control  tasks  have,  however,  been  "explained"  In  three  possibly  equally 
likely  ways  (Ref.  34) I 

a.  Random  "noise,"  with  a  mean-squared  value  proportional 
to  mean-oquarcd  linear  output,  superimposed  on  the 
operator's  linear  output 

b.  Nonsteady  operator  behavior,  that  is,  variation  of  the 
operator' e  describing  function  during  the  course  of  a 
measurement  run 

o.  Nonlinear  anticipation  or  re lay -like  operation,  super¬ 
imposed  on  the  operator's  linear  output 

Remnant  data  are  far  more  sparse  than  describing  function  information; 
and  the  data  available  are  not  especially  reliable.  Consequently,  the 
summary  statements  above  are  about  all  that  should  be  said  in  connection 
with  the  circa  i960  model. 

3.  Rationale  for  the  Ibdatlng  Model 

Zt  may  be  well  to  remark  that  the  analytical-verbal  model  is  in  part 
a  hypothesis  which  derives  from  what,  after  all,  can  only  be  described 
as  limited  data.  The  rationale  of  the  model,  however,  does  not  rest 
exolueively  on  observation.  The  adjustment  rules  are  partly  an  expres¬ 
sion  of  praotieal  synthesis  procedures  for  inanimate  servomechanisms 
(Rei’s.  7  >  22,  and  26). 

Desirable  properties  oi  a  "good"  feedback  control  system  are  to: 

a.  Provide  specified  command-response  relations!,!.;^ 

b.  Suppress  unwanted  inputs  and  disturbances 

c.  Reduce  effects  of  variations  and  uncertainties  in 
elemenbB  of  the  control  loop 
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It  can  be  shewn  (Ref.  7,  for  example)  that  these  three  functions  are 
accomplished  in  a  single-loop  manual  or  automatic  feedback  control 
system  by  making  the  complex  gain  of  the  open-loop  system,  YpY0,  very 
large  over  the  range  of  frequencies  In  vhloh  the  eorafland  Input  and  load 
disturbances  have  substantial  components,  and  very  small  outside  this 
range. 

For  the  unity  feedback  system  of  Fig.  2,  the  desired  closed-loop 
transfer  function  is  the  transfer  function  of  a  lew  pass  filter.  The 
entire  range  of  positive  real  frequencies  may  then  be  divided  into  three 
regions  of  principal  interest  in  terms  of  the  magnitude  of  the  open-loop 
system  transfer  function.  These  are  displayed  in  Table  IV. 


TABU!  IV 

OPEN-  AND  CLOSED-LOOP  TRANSFER  FUNCTIONS 
OF  A  "GOOD"  SERVOMECHANISM 


FREQUENCY 

OPEN-LOOP 

CLOSED-LOOP 

TRANSFER  FUNCTION 

TRANSFER  FUNCTION 

S.  «  ! 

“b 

IVcl  »  1 

-Ml.  i  , 

1  +YpY0 

iL  A  1 

■  row 

±  »  1 

IVcl  <<:  1 

Wo  '  . 

1+ Vo  "  YPY= 

The  region  near  the  crossover  frequency,  o^,  (where  | YPY0 |  A  1)  has 
a  surpassing  importance  in  the  synthesis  of  feedback  systems.  First, 
for  "good"  performance  the  crossover  frequency  must  exceed  the  largest 
frequency,  dj.,  at  which  there  are  appreciable  components  of  the  commands 
or  external  disturbances,  <i^  is,  in  this  sense,  somewhat  greater  than 
an  upper  bound  on  the  frequency  region  wherein  |YpY0|  »  1,  a  condition 
whioh  proyipea  good  following  of  the  system  input  by  the  output  and 
suppress! dm  of  disturbance  effects.  Second,  the  shape  of  YpYc  at  and 
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near  crossover  frequency  determines  the  dynamics  of  the  dominant  modes 
of  By  B  Lew  itia^uuDc »  The  uunuitiun  for  uvubml  ouoiULabory  stability  is 

|YpY0|  *  1  and  <YpYc  *  -«  when  sgn  YpYo(0)  «  +  (7) 

S 

Hence,  neutrally  stable  or  unstable  dominant  modes  are  moot  often  avoided 
by  adjusting  the  system  so  that 

|YpYc|  <  1  when  4 YpYc  «  -it 

and  (8) 

*YpY0  >  -it  when  |YpYc|  -  1 

These  are  the  often  quoted  conditions  of  positive  gain  margin  and  phase 
margin  (e.g.,  Ref.  22),  They  are  an  egression  for  single  minimum  phase 
or  minimum  phase  plus  transport  lag  systems  of  the  Nyquist  stability 
criterion.  Since  stability  is  quite  literally  the  sine  qua  non  of  a 
feedback;  control  system,  these  requirements  for  gain  and  phase  margin 
nay  severely  restrict  the  choice  of  the  crossover  frequency. 

The  phase  angle  at  a  frequency  <%  of  a  transfer  function,  YpY0,  which 
contains  a  transport  lag,  t,  but  which  is  otherwise  minimum  phase,  in 
terms  of  amplitude  ratio  slopes  is  (RefB.  7  and  £t) 


where  the  slopes,  d|YpY0|/d[ln(co/u*j)]  ,  are  expressed  in  db/decade.  As 
illustrated  by  Fig.  3,  the  In  coth  (l/2) |ln(aj/c^) |  term  in  the  Integral 
applies  a  large  weighting  to  slope  changes  in  the  immediate  vicinity  of 
and  greatly  attenuates  the  effeots  of  the  integrand  of  slope  changes  else¬ 
where.  Consequently,  the  phase  at  is  affected  primarily  by  tc%,  by  the 
local  db  amplitude  ratio  slope  (the  second  term  in  the  expression),  and  by 
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Figure  5*  Weighting  Function  for  Bode's  Amplitude  Ratio  Slope: 
Phase  Relationship 


local  changes  in  this  Blope  (the  integral  term) .  If  the  db  amplitude  ratio 
slope  is  essentially  constant  over  a  wide  region  about  ci^,  the  expression 
reduces  approximately  to  the  second  term  alone  plus  the  transport  lag's 
contribution.  In  this  event  the  phase  associated  with  a  constant  amplitude 
ratio  slope  of  -20n  db/decade  will  be  simply  —rut  —  nx/2  rad. 

For  low  pass  open-loop  transfer  functions  the  amplitude  ratio  slope  at 
gain  crossover  is  negative,  so  a  positive  phase  margin  can  usually  exist 
only  when  d[YpYc^.  /d[ln  (cn/a^.)]  in  the  immediate  vicinity  of  crossover  is 
less  (numerically)  than  -40  db/decade,  the  local  changes  in  Blope  are  moder¬ 
ate,  and  the  to\,  contribution  is  minor.  The  available  crossover  regions  for 
most  transfer  functions  are,  therefore,  confined  to  areas  where  the  local 
amplitude  ratio  elope  fulfills  theoe  conditions,  and  the  choice  of  crossover 
frequency  ie  delimited  accordingly. 
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These  considerations  are  embodied  in  the  Primary  Rule  of  Thumb  of 
fi-equexxuy  rc ayuuoa  synthesis  i  "Find  ux  uica tc  a  falx-  it  txt>  Lull  of  — £0  ub/ueutule 
slope  for  the  amplitude  ratio,  and  then  make  It  the  crossover  region  by 
putting  the  o  db  line  through  it"  [i.e.,  make  the  gain  such  that  | YpYc ( jo>) |  -  i 
where  YpY0(jffi)  *  1  /j(<b/c»c)] .  This  orude  prescription  for  stability  and  good 
response  is  generally  adequate  for  minimum  phase  systems »  It  can  be 
extended  directly  to  the  transport  lag  case  by  adding  a  prescription  for 
a  positive  phase  margin.  Typically,  the  phase  margin,  flfyj,  for  a  well- 
adjusted  regulator  is  approximately  30-40°;  a  somewhat  higher  value  is 
customarily  used  for  servomechanisms,  which  must  follow  commands  as  well 
as  suppress  disturbances.  The  human  has  his  own  ideas  (Adjustment  Rule  4b) l 

Finally,  then,  for  a  "good"  feedback  control  system  the  operator'^ 
describing  function,  Yp(j<u),  must  be  adjusted  so  that  the  crossover 
frequency,  u\,f  exceeds  the  highest  important  frequency  in  the  input,  ay.; 
and  so  that  lpYc( Jcc)  conforms  rt  low  frequencies,  high  frequencies,  and 
crossover  frequency  to  the  requirements  noted  in  Table  IV. 

Many  of  the  above  remarks  about  the  rationale  of  equalization  adopted 
by  the  pilot  can  be  made  more  concrete  and  understandable  by  the  defini¬ 
tion  of  an  approximate  "crossover  model"  for  manual  control  systems. 

This  has  been  done  in  Ref.  38,  where  it  iB  pointed  out  that  considera¬ 
tion  of  the  requirements  of  "good"  feedback  system  performance  leads 
directly  to  the  conclusion  that  the  pilot  adjusts  his  describing  function 
so  that  the  open-loop  function,  YpYc,  in  the  vicinity  of  the  gain  cross¬ 
over  frequency,  (%,  is  closely  approximated  by 


Vc 


o^e 


-Jutte 


J» 


(10) 


This  crossover  model  is  not  a  replacement  for  the  analytical-^verbal  model, 
but  is  instead  a  convenient  approximation  suitable  for  many  engineering 
purposes.  While  it  is  a  better  description  of  amplitude  ratio  character¬ 
istics  than  of  phase  characteristics,  it  often  describes  the  most  signif¬ 
icant  features  of  operator  behavior  adequately.  This  is  because  the 
actual  shape  of  the  open-loop  function  away  from  the  gain  crossover 
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frequency  is  usually  almost  irrelevant  to  the  closed-loop  performance. 

It  is,  therefore,  often  unnecessary  to.  retain  terms  in  the  describing 
function  whose  influence  is  not  felt  in  the  immediate  vicinity  of  gain 
crossover. 

Table  V  shows  the  application  of  the  adjustment  rules  to  the 
prediction  of  a  pilot's  describing  function  and  to  the  formulation  of 
the  crossover  model  for  several  simple  limiting-case  controlled  elements.* 


TABLE "V 

OPERATOR  DESCRIBING  FUNCTION  AND  CROSSOVER  MODEL 


WHEN  THE 

CONTROLLED  ELEMENT 
TRANSFER  FUNCTION  IS 

<Y0) 

THE  OPERATOR'S 
DESCRIBING  FUNCTION  IS 

(Yp) 

THE  CROSSOVER 
MODEL 

(cn  about  o\,) 

IS  YpYc 

^c_ 

J » 

Kpe-^e 

K^e-JaWe 

Kc 

Kps'3"'  , 

KpKc  e-J^e 

(Tj;Jtu+  1 )  i  Tl  <C  ^ 

Ti  JU3 

Kc 

Kpe“jajr*(TLJ<o+l)  ;  ~  «  at 

KpKcTLe“JtOT« 

(»2 

jo 

In  many  instances  more  complicated  controlled  elements  can  be  satisfac¬ 
torily  approximated  in  the  vicinity  of  gain  crossover  by  theBe  simple 
limiting  cases.  Thus,  for  example,  a  damped  second-order  controlled 
element,  Y0  -  Kc/[(jV<%)2+  (2§/%)  J<u+ l]  >  with  a  natural  frequency, 
<%«  Oq,  is  closely  approximated  near  by  r  double  integrator, 

Yc  -  Kco^/(jcn)2.  The  operator's  describing  function  would  then  include 
a  lead  equalization  (see  also  the  last  entry  in  Table  III). 


*Notr  that  the  neuromuscular  lag,  which  has  an  important  influence 
at  and  near  crossover  frequencies,  is  subsumed  in  a  larger  effective 
reaction  time  delay,  te  »  t  +  T$j  for  Kc/(jo>)  •  For  the  other  two  con¬ 
trolled  elements  the  lead  equalization  can  occur  at  high  frequencies, 
so  for  these  ■  T-!  Tto-Tl.  Actually,  te  is  a  catchall  term  which 
serves  to  incorporate  into  the  describing  function  the  phase  effects 
near  crossover  of  pure  transport  lags  and  poles  and  zeros  vitlj  break¬ 
points  much  greater  than 
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b.  oBJxort  Of  m  vamaam 

The  primary  purposes  of  the  experimental  series  are  the  ■validation 
of  the  existing  analytical-verbal  model  and  the  extension  of  this  model 
In  accuracy  and  detail. 

1.  Motel  VUJtetitt 


In  the  sense  in  which  It  Is  used  here,  "validation"  is  intended  to 
mean  the  confirmation  of  the  extrapolations  based  on  limited  data.  The 
delineation  of  crucial  experiments  designed  to  test  the  validity  of  the 
model  involve  the  following  preliminary  predictions: 


a. 


Describing  function  form  adjustment.  For  a  variety  of 
controlled  elements  the  adjustment  rules  should  be  used 


to  predict  the  operator's  describing  function  form.  These 
controlled  elements  should  be  such  as  to  require  a  complete 


range  of  operator  equalization  form  adaptation;  and  they  should 
be  "new"  in  that  they  have  not  previously  been  examined  experi¬ 


mentally  in  detail. 


b.  Mean-squared  error  minimization  and  variation.  Estimates 
of  mean-squared  error  due  to  system  forcing  function,  and 
its  variation  with  forcing  function  characteristics,  Bhould  be 
derived.  Special  control  situations  should  be  selected  which 
most  concretely  illuminate  the  experimental  consequences  of 
describing  funotion  adjustment  for  minimum  mean-squared  tracking 
error,  and  the  variation  of  this  quantity  with  farcing  funotion 
modifications . 


With  these  preliminary  analyses  and  predictions  in  hand,  an  experimental 
series  can  readily  be  delineated  to  provide  data  which  either  confirm  or 
deny  the  predictions  of  the  model. 

As  already  indicated  in  Table  III,  the  simple  controlled  elements 
Y0  -  Ko,  Ko/jm,  and  Kc/Qco)2  evoke  a  complete  range  of  equalizer  form 
adjustment,  so  these  become  prime  candidates  for  oritical  validation 
experiments.  Presuming  that  Elkind's  experiments  (Ref.  1?)  with  a 
variety  of  forcing  functions  for  the  simple  controlled  element  Y0  ■  Kc 
could  be  repeated,  in  part,  in  the  current  Beriei,  then  his  data  as  a 
whole  could  be  considered  in  direct  context  with  our  own.  The  crucial 
controlled  element  form  for  model  validation  would  then  be  Yc  ■  Kq/Jco 
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and  Kc/( Jo>)2.  Limited  teste  have  been  accomplished  for  the  first  of 
theee  (Refs.  24  and  47),  but  no  data  are  available  for  the  Becond. 

These  controlled  element  forms  are  further  attractive  in  that  they 
represent,  in  the  vicinity  of  crossover,  limiting  oases  of  a  great 
variety  of  other  possible  controlled  elements  of  far  greater  complexity* 

As  another  aspect  of  validation,  it  is  necessary  to  distinguish 
between  those  features  of  the  model  in  any  particular  situation  which 
must  be  evidenced,  and  which  therefore  would  tend  to  be  invariant  with 
different  operators,  and  those  features  which  might  be,  to  an  extent, 
a  matter  of  Indifference.  It  will  become  apparent  later  that  individual 
operators  may  display  what  can  only  be  described  aB  a  particular  tracking 
style.  The  style  of  one  operator  may  be  somewhat  different  from  the 
style  of  another  operator,  and  this  is  reflected  in  changes  in  the 
describing  function  model  in  those  regions  where  the  form  of  the  model 
is  not  critical  to  good  tracking  performance.  On  the  other  hand,  where 
the  describing  function  form  is  critical  to  good  tracking  performance, 
most  particularly  in  the  vicinity  of  the  crossover  frequency,  it  is  to 
be  expected  that  the  operator  will  not  exercise  a  choice  and  that  the 
tracking  performance  will  be  tightly  constrained.  The  three  controlled 
elements  mentioned  above  are  not  particularly  constraining  except  in  the 
region  of  crossover.  However,  certain  unstable  controlled  elements  tend 
to  require  a  more  uniform  behavior. 

8.  Model  Extension 

The  primary  limitations  to  the  existing  model  are  due  to  limitations 
in  the  experimental  data  on  which  the  model  is  based.  In  essence,  the 
model  is  about  aB  sophisticated  as  it  can  be  without  further  empirical 
knowledge.  So,  the  question  becomes  "Where  are  the  existing  data 
deficient?" 

a.  The  vast  majority  of  the  data  are  based  on  forcing 
function  bandwidths  Which  vere  very  low  In  frequency 
relative  to  ti^..  Hence,  most  of  the  important  crossover, 
features  of  the  model  are  based  on  extrapolations  and  the 
doctrine  of  compatibility  between  date  sources. 
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hr  Th»»*  lew  frequency  f  creins^f unction  bandwidth  data  iu 
not  disclose  any  differences  between  manipulators. 

c.  To  achieve  random-appearing  forcing  functions,  noise 
generators  or  their  equivalent  were  used.  (This  choice 

of  forcing  funatlon,  while  satisfying  the  random-appearing 
criterion,  Introduced  other  problems: 

(1)  Sampling  variability  (i.e.,  effects  of  finite  run 
length)  Introduced  several  uncertainties. 

(2)  Data  processing  equipment  was  complex  and  often 
notoriously  unreliable. 

(3)  The  remnant,  of  course,  came  out  as  a  continuous 
power  spectral  density.  Consequently  It  was  of 

such  form  that  it  was  not  possible  to  discriminate  among  the 
several  possible  remnant  sources  (i.e.,  remnant  can  be  due  to 
sampling,  nonllnearltles,  nonstationary  behavior,  actual  random 
noise  injection,  etc.,  but  the  distinctive  features  correspond¬ 
ing  to  each  of  these  tended  to  be  lost  In  a  general  snoothed- 
over  slush).  Also,  precious  little  remnant  data  were  available. 

d.  With  the  notable  exception  of  Elkind's  results,  the  effects 
of  variation  in  forcing  function  bandwidth  were  not  system¬ 
atically  explored.  Also,  the  consequences  of  changes  In  forcing 
function  amplitude  on  open-loop  describing  function  data  were  not 
known  in  detail. 

Many  other  data  deficiencies  exist,  in  both  scope  and  kind,  but  the 
four  major  items  listed  above  were  the  most  influential  in  the  program 
planning  phases  of  the  experiments  to  be  reported  here.  A  major  correc¬ 
tion  of  these  data  deficiencies  was  desired,  at  least  in  the  areas  most 
significant  for  man-machine  system  analysis.  Thus,  beyond  the  critical 
experiments  for  the  validation  of  the  analytical-verbal  model,  it  was 
expected  that  experimental  phases  of  the  program  would  yield  additional 
data  which  would  permit: 

a.  The  refinement  of  statements  concerning  the  adjustment 
rules  and  the  form  of  the  describing  function  in  the 
vicinity  of  the  crossover  frequency.  Such  data  might  Include: 

(1)  Phase  margin  data  (and  gain  margin  data  whore 
this  parameter  is  pertinent) 

(2)  Closed-loop  dominant  mode  characteristics,  suoh 
as  olosed-loop  damping  ratio, 
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(3)  Tighter  limits  on  vp'ues  of  t  and/or  t  +Tj^ 

(4)  Better  understanding  of  and  limits  on  the  u^, 
regression  condition 

(p)  Better  limits  on  the  maximum  value  of  lead 
time  constant,  T^  and  lag  time  constant,  Tj 

b.  An  inproved  understanding  of  the  general  effeots  of 
forcing  funotion  amplitude  and  bandwidth • 

c.  An,  in®) roved  understanding  of  the  possible  effects  on  the 
pilot's  describing  function  of  changes  in  the  manipulator. 

d.  An  extension  of  the  population  over  which  the  operator's 
describing  function  was  measured,  (it.  would,  of  course, 

enlarge  confidence  in  th«  model  if  it  could  be  shown  that  not 
Just  a  few  subjects,  but  many,  could  be  described  in  the  same 
way.) 

A  natural  result  cf  an  experimental  program  carried  out  for  these  purposes 
would  bo  a  set  of  definitive  data  measured  over  a  wide  range  of  frequencies 
for  elementary  and  limiting-case  ontrolled  elements. 

In  addition  to  all  the  purposes  described  above,  it  was  hoped  that 
the  experiment.  1  results  would  yield  the  information  which  would  permit 
a  variety  of  questions  concerning  human  operator  performance  in  tracking 
tasks  to  be  answered.  Some  of  the  questions  which  have  arisen  time  after 
time  in  connection  with  previous  investigations  are  suggested  by  the 
following  phrases: 

*  a.  Motor  response  models:  Pew,  if  any,  investigators  using 
random-input  describing  function  techniques  have  carried 
their  measurame:,  ts  to  frequencies  high  enough  to  reveal  any  but 
the  crudest  facts  concerning  the  motor  response  (neuromuscular 
lags)  of  the  human  operator. 

b.  Statlonarity  during  runs:  Data  rsduction  performed  in 
the  frequency  domain  requires  fairly  long  averaging 
times.  The  describing  function  which  is  measured  is  then  an 
average  describing  function  for  the  period  of  the  test.  It 
may  be  postulated  that  the  operator's  perfonnance  is  far  from 
invariant  over  a  period  of  minutes,  and  that,  for  example,  his 
describing  function  might  be  appreciably  different  at  the  end 
Of  a  long  run  than  at  the  beginning.  If  this  were  the  case, 
there  would  be  good  reason  to  doubt  the  validity  of  the  mathe¬ 
matical  description  under  any  circumstances  other  than  as  an 
expression  of  the  average  performance  in  a  particular  experimental 
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situation.  The  nature  of  oross-spectral  data  reduction  techniques 
and  the  forcing  function  frequencies  of  interest  preclude  the 
explicit  examination  of  variations  in  the  operator's  describing 
function  for  any  period  shorter  than  approximately  2  min,  but 
evidence  for  or  against  rapid  fluctuations  in  the  description  of 
the  operator  could  be  Indicated  by  considering  p  values  in  con¬ 
text  with  other  measurements.  Also,  the  question  of  whether  or 
not  the  description  would  be  appreciably  different  for  the  first 
2  min  and  the  last  2  min  of  the  10-min  run  could  quite  easily  be 
examined. 

*  c.  Pulsing  control  as  ar  expression  of  lead;  It  is  a  matter 

of  fairly  common  experience  that  in  control  tasks  which 
are  very  difficult,  either  because  of  potential  instability  of 
the  controlled  element  or  because  of  a  very  wide  bandwidth  of 
the  input  forcing  function,  the  operator  will  tend  to  exercise 
control  with  e ■  series  of  discrete  pulses.  The  reasons  why  he 
does  this  are  not  clear,  but  it  may  be  that  hr  can  in  thiB 
fashion  arrange  more  lead  or  phase  advance.  This  behavior, 
however,  has  some  nonlinear  properties  and  is  a  possible  source 
of  remnant. 

#  d.  Sampling  effects :  Some  authors  have  suggested  (Refs.  5, 

6,  and  66)  that  the  human  may  beh*  e  as  a  samp led- da t 
system  with  a  reasonably  constant  sampling  rate  for  otherwise 
stationary  conditions.  Physical  evidence  for  or  against  thiB 
hypothesis  is  very  sparse.'  If  the  human  operator  c.pproximu+es 
a  constant- sampling-period,  sampled-data  subsystem,  the  offset 
of  sampling  would  be  a  directly  testable  physical  source  of 
the  remnant. 

#  e.  Nonlinear  effects:  There  is  evidence  of  nonlinear 

behavior  in  both  the  perceptual  and  actuation 
meohanisms  of  the  operator.  For  example,  there  are  rate 
thresholds  below  which  the  perception  of  motion  is  not 
possible.  Further,  in  attempting  to  make  predictions  of 
future  motion  at  rates  above  the  threshold,  the  operator 
will  overestimate  slow  rates  and  underestimate  fast  ones. 

Some  such  nonlinear  effectB  might  well  be  a  suitable  source 
of  the  remnant. 

*  f.  Other  sources  of  remnant:  In  addition  to  Items  a-e 

aSove,  a  variety  o i  mechanisms  might  possibly  be  the 
source  of  the  remnant. 

It  waB  not  thought  desirable  to  attempt  to  resolve  all  of  these 
questions  within  the  framework  of  the  present  program,  but  certain 
experiments  were  planned,  to  be  oarrled  out  more  or  less  explicitly,  to 
resolve  some  of  them.  (Those  which  are  treated,  to  some  extent,  in  this 


report  are  narked  with  an  asterisk  in  the  margin.)  All  the  questions 
and  data  deficiencies  listed  were  carefully  considered  in  evolving  the 
final  experimental  plans,  and  where  at  all  possible  attempts  were 
Included  to  attaok  suoh  questions  within  the  main  framework  of  the 
planned  experiments. 
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Pragmatically,  the  most  compelling  justification  of  any  model  Is  Its 
capacity  to  subsume  past  experimental  results  and  to  predict  the  outcome 
of  future  experiments  especially  contrived  to  exercise  the  model  to  its 
limits.  The  existing  human  pilot  model,  discussed  In  the  last  chapter, 
was  constructed  on  the  basis  of  compatibility  with  past  results,  so  it 
performed  the  first  function  noted  prior  to  the  initiation  of  the 
current  program.  Its  next  test  was  in  forecast  of  extended  situations. 
These  predictions  were  the  basis  of  much  of  the  planning  of  the  experi¬ 
ments  and  choice  of  particular  experimental  situations. 

Three  kinds  of  predictions  are  summarized  here.  The  first  is, 
fundamentally,  application  of  the  adjustment  rules  and  rationale  to 
three  simple  controlled  elements.  Although  one  of  theBe,  Yc  ■  Kc,  has 
been  extensively  studied,  it  is  treated  again  here  with  the  model  to 
provide  the  basis  for  comparative  statements  between  this  and  the  other 
two  systems.  The  second  kind  of  forecast  is  also  concerned  with  describ¬ 
ing  function  adjustment,  but  in  a  rather  special  way.  The  intent  was  to 
find  epeoial  controlled  element  forms  which  would  tend  to  tightly  con¬ 
strain  the  operator's  choice  of  characteristics.  To  thiB  end  the  model 
is  used  to  explore  possible  control  situations  which  will  tend  to  con¬ 
firm  or  deny  conjectures  about  variability.  Finally,  the  third  type  of 
prediction  uses  the  model  to  make  mean-squared-error  estimates  for  sub¬ 
sequent  experimental  validation.  In  pursuing  this  objective,  unexpected 
results  vers  obtained  which  provide  the  raison  d'etre  for  two  adjustment 
rules  which  previously  lacked  a  theoretical  basis. 

A.  EQUALIZATION  ADJUSTMENT 

A  series  of  controlled  elements  which  require  the  pilot's  equalization 
selection  to  range  from  lag-lead  to  essentially  pure  .gain  to  lead-lag  have 
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already  been  diacuesed  in  the  last  chapter.  From  there  It  will  be  | 

recalled  that  these  pilot-adapted  forms  correspond  to  Yc(  jcc)  forms  of  | 

Kq,  Kq/JcDj  and  Kc/(j«})2,  reepeotlvely.  Detailed  analyses  of  man-machine  ! 

systems  involving  these  controlled  elements  will  be  sunsarlKed  below. 

1 .  y0  -  *e 

This  simplest  possiblo  controlled  element  is  also  the  most  extensively 
studied,  since  it  was  included  as  one  of  the  controlled  elements  in  three 
of  the  programs  listed  in  Table  11. '  The  most  definitive  investigation  was 
that  conducted  by  Elklnd  (Ref.  15).  Yet  some  characteristics  remained 
ill  defined  (e.g.,  phase  margins)  or  not  too  well  understood  (e.g., 
ufc  regression). 

The  operator-adapted  describing  function  fur  control  of  the  pure  gain  ’ 

controlled  element  will  be  5 


The  phase  margin  of  this  system  will  be 

<Pm  -  #  +  <  Vo(J«b) 

■  «  -  tan-1  TjOfc  -  tu^j  -  tan-1  TflOj,  +  tan-1  TlU\,  ( 1 2) 

The  lag  time  constant,  Tj,  will  be  much  greater  than  either  the  lead 
time  constant,  T&,  or  the  neuromuscular  lag,  to  achieve  good  low 
frequency  closed-loop  system  response.  In  fact,  the  lead  and  the 
neuromuscular  lag  will  occur  at  relatively  high  frequencies  and  they 
will  tend  to  cancel.  Therefore,  at  frequencies  in  the  crossover  region 
and  below,  the  primary  effect  of  these  two  terms  will  be  seen  in  their 
contribution  to  the  phase,  which  will  be  approximately 

tan"1  T ip  -  ban"'1  -  (Tj,  -  TrJoj  ,  <0  £  O(ccfc) 

where  0(  )  means  "of  the  order  of." 
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The  net  effect  is  thus  Indistinguishable  at  low  frequencies  from  a 
modified  "effective”  reaction  time  delay,  xe>  given  by 

Te  4  *+TN-TL  (13) 

Using  this  approximation,  the  open-loop  system  characteristics  become 
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(14) 


A  generic  Jm  Bode  plot  showing  the  open-loop  characteristics  of  Eq  14 
is  given  in  Fig.  4.  Also  illustrated  there  is  the  crossover  model 
approximation  obtained  by  noting  that  l/Tio*  «  1 ,  so  that  KpKQ/Tj  -  a*,. 
The  crossover  model  version  of  Eq  14  is  then 
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Because  oie  »  l/Ti  for  good  low  frequency  closed-loop  response,  the  phase 
margin  expression  (Eq  12)  au  modified  by  Eq  13  becomes  approximately 


to  4  rt“  (T-fffe)  -t«°* 

4  *§■  +  P  “  ''e^ 


06) 


where  p  ■  l/Tjo^.  The  crossover  model  phase  margin  will  be  the  same  as 
that  given  by  Eq  16  when  5  is  set  equal  to  aero. 

In  this  system  the  moBt  important  parameter  is  the  crossover  frequency. 
Unfortunately  is  also  one  of  the  more  difficult  things  to  estimate 
accurately  using  the  existing  model  because  of  the  absence  of  extensive 
phase  margin  data.  Consequently  an  estimate  of  a  likely  crossover 
frequency  range  is  usually  the  beat  that  can  be  done.  In  the  present 
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case  this  can  be  approached  as  follows.  Solving  for  the  crossover 
frequency  from  Eq  1 6, 

i  +  (17) 

Because  Tjj  and  Tl  approximately  cancel  for  this  case,  t0  will  approach  T. 
A  reasonable  value  for  re  is  then  about  0.20  sec,  which  corresponds  to 
the  higher  side  of  the  r  range  given  in  Chapter  II.  With  a  value  for 
selected,  j.  series  of  bounds  on  can  be  established  using  various 
values  of  3  and  qfy  in  Eq  17.  Any  value  of  3  less  than  0.1  is  compatible 
with  good  closed-loop  system  response,  so  a  3  ■  0.1  line  provides  one 
boundary.  The  system  will  be  unstable  if  the  phase  margin  is  less  than 
zero  degrees,  so  an  upper  bound  is  established  by  the  line 

%ax  *  (18) 

Finally,  an  upper  limit  for  the  phase  margin  will  be  taken  as  40°. 

(This  estimate  is  based  on  the  very  limited  Elkind  data  in  the  crossover 
region.)  The  line, 


established  a  lower  boundary  on  the  region  of  likely  crossover  frequencies. 

The  boundaries  developed  above  are  illustrated  in  Fig.  5.  Using  the 
middle  of  the  region,  an  estimated  crossover  frequency  will  be  about 
6.5  rad/sec. 

Boundaries  similar  to  those  shown  in  Fig.  5  can  be  developed  using  a 
constant  phase  margin  and  a  variable  ts  (Bee  Ref.  37)*  Uncertainties 
exist  in  both  r9  and  the  phase  margin  (which  are  ]io£  independent  quanti¬ 
ties — see  Adjustment  Rule  4b),  so  the  bounded  region  can  be  considered 
to  take  into  reasonable  account  nominal  variations  in  both  quantities. 
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Figure  5*  Crossover  Frequency  Estimation  for  Yy  ■  K0 


According  to  Adjustment  Rules  5a,  b ,  and  o,  the  crossover  frequency 
estimated  above  will  be  essentially  Independent  of  variations  in  Ky,  and 
will  also  be  constant  in  the  presence  of  increases  in  forcing  function 
bandwidth  until  <Cj_e  (Eq  5)  approaches  or  exceeds  6.5  rad/sec  or  so. 

When  exceeds  the  estimated  the  observed  crossover  frequency  will 
exhibit  a  substantial  decrease  (cc^.  regression). 
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a.  y0  -  Vi® 


This  single-integrator  controlled  element  is  an  idealised  rate  control, 
i.e.,  an  operator-applied  step  function  input  to  the  machine  results  in  an 
output  rate  proportional  to  the  step  input.  Because  of  the  free  s,  there 
is  no  requirement  for  lew  frequency  equalization,  so  the  adjustment  rules 
indicate  Tj  •  0.*  At  high  frequencies  the  operator  lead  can  be  used  to  com¬ 
pensate  for  the  neuromuscular  lag.  Thus  the  appropriate  open-loop  .system 
model  is 
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where  Te  ■  t  +  Tjj  -  T^.  This  complete  open-loop  system  model  is  identi¬ 
cal  to  the  crossover  model  (compare  with  Table  V) .  The  open-loop  system 
Bode  diagram  is  essentially  the  same  as  that  for  the  crossover  model  Bode 
diagram  for  the  Yc  ■  Kq  case  (dashed  curves  in  Fig.  4). 

Because  conditions  in  .the  region  of  crossover  are  so  similar  to  those 
for  the  Yc  ■  Kc  situation,  the  best  estimate  of  a*,  possible  using  the  circa 
i960  model  is  based  on  that  given  in  the  laBt  article  presuming  p  does  not 
appear  in  Eq  16.  Because  p  is  thus  effectively  zero,  the  estimated  mean 
might  be  dropped  to  about  6  rad/sec  or  so  by  using  the  average  of  the 
and  bounds  evaluated  for  zero  p.  [A  few  higher  phase  margins  have 

been  observed  (Ref.  24)  for  some  Yc  ■  Kc/jco  data,  so  the  lower  bound  should 
also,  perhaps,  be  reduced.]  As  far  as  the  numerical  values  are  concerned 
(6.5  and  6  rad/sec),  this  is  a  refinement  that  is  probably  not  warranted 
because  of  the  basic  uncertainties  present  in  both  estimates.  But  the 
relative  magnitudes  of  a*,  for  Yc  -  Kq  and  Kc/joa  are  expected  to  exhibit 
this  ordering. 

The  application  of  Adjustment  Rule  5  to  this  case  is  the  same  as  noted 
for  the  Yc  >  K0  situation. 


5 


Tfa  -  Xo/(J“)2 


The  double-integrator  controlled  element  provides  a  crucial  teat  for 
the  validation  of  the  analytioalrvetfbal  model  ainoe  experimental  data 
bad  never  been  obtained  parlor  to  the  current  program.  Thus  the  applica¬ 
tion  of  the  adjustment  rules  represents  pure  extrapolation  of  the 
adaptation  principles  believed  relevant  to  human  operator  description. 


The  operator  describing  function  adapted  for  the  control  of  the 
Ke/(jto)2  controlled  element  needs  no  low  frequency  equalization  (Tj  -  0), 
but  does  require  the  use  of  some  lead  equalization  to  meet  the  simple  and 
primary  requirement  of  stability.  Thus,  in  the  complete  open-loop 
describing  function 


?  ~  (Jus)2  (%ja>  +  1) 


(21) 


it  is  immediately  apparent  that  Tj,  >  t  +  Tj^  is  a  necessary  condition  for 
stability.  A  low  frequency  approximation  to  Eq  21  is 


«.  * 
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(22) 


where  Te  «  t  +  Tjj.  This  approximate  expression  is  suitable  for  the  low 
frequency  range  up  to  and  somewhat  beyond  crossover.  Note  especially 
that,  because  the  equalizing  lead  is  needed  at  much  lower  frequencies 
to  obtain  stability,  t9  here  does  not  contain  the  compensating  effect 
of  a  high  frequency  lead  tending  to  offset  the  neuromuscular  lag. 
Consequently,  re  for  the  Yc  »  Kc/(jcu)2  case  will  be  considerably  larger 
than  those  values  estimated  for  the  Y0  ■  Kq  and  Y0  -  Kq/Jco  systems. 

Note  also  that  Eq  22  approaches  the  crossover  model  (Table  V)  only  when 

Tl%  »  1  • 

To  gain  a  better  appreciation  for  the  likely  position  of  l/T^, 
consider  the  detailed  unified  servo  analysis  (Ref.  59)  presentation  shown 
in  Pig.  6.  Here  l/T^  is  located  at  a  frequency  which  is  considerably 
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leas  than  the  crossover  frequency.  The  closed-loop  system  in  this  case 
vill  save  a  law  frequency  characteristic  given  by 


(23) 

The  approximate  cancellation  of  the  low  frequency  closed-loop  lag,  T[, 
by  the  lead  equalization  time  constant,  T^,  indicated  in  Eq  24  results 
in  an  excellent  low  frequency  closed-loop  response  characteristic. 

Now,  by  way  of  contrast,  Imagine  that  the  lend  break  point  is  much 
closer  to  l/te  and  the  crossover  frequency.  The  dipole  pair  in  Eq  23 
will  then  no  longer  nearly  cancel.  Also,  for  separations  between  l/Tj, 
and  i/t0  of  about  6  (tlAb  -  3  +  VS  -  5  -83,  obtained  from  the  require¬ 
ment  that  two  real  values  of  d[YpYc(-a)]/da  -  o),  the  two  poles  starting  at 
the  origin  can  never  reach  the  real  axis.  Further,  the  closed- loop 
characteristics  become  increasingly  sensitive  to  variatione  in  Tl  and 
as  their  separation  becomes  less.  Consequently,  it  is  clear  that  1/Tl 
should,  ideally,  be  located  at  values  much  less  than  ctfe.  This  completes 
the  prediction  of  the  detailed  form  of  Yp  for  this  case. 

The  next  problem  is  to  estimate  the  approximate  crossover  frequency. 
This  will  be  accomplished  in  a  fashion  similar  to  that  used  for  Y0  -  Kc» 
The  phase  margin  corresponding  to  Eq  22  is 

cpM  .  tan-1  Tlc^  -  Teofc  (25) 

which,  if  1  /\%  «  1  >  becomes  approximately 

«PM  ;  T  " 

i  Y  -  p  -  vfc  (26) 
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where  p  ■  l/T^c%.  Rearranging  Eq  26,  an  expression  for  the  crossover 

f rtn^uouujr  1b  ubuniiiou; 


% 


_  1  IJK  - 

-  t;It  - p  “ 


(27) 


This  equation  is  similar  to  Eq  17  for  the  Y0  ■  Kq  controlled  element, 
and  the  likely  range  for  cuq  values  will  be  determined  In  the  same  general 
way.  The  pertinent  boundaries  in  this  case  are: 


p 

<  0.1 

,  for  good  low  frequency  response 

%ax 

■  ii-r-e) 

,  for  the  -  0°  bound 

%0 

•  4GH 

,  for  the  <py  ■  40°  bound 

Before  these  boundaries  can  actually  be  drawn,  a  value  for  re  must  be 
estimated.  A  well-established  and  representative  number  for  t  is  0.2  sec 
to  which  must  be  added  an  estimate  for  Tjj.  As  noted  in  Chapter  II,  Tjf  1b 
variable  over  a  fairly  wide  range,  but  a  typical  value  of  0.1  is  often 
used  as  representative  of  the  order  of  magnitude.  This  will  be  the  pro¬ 
cedure  used  here.  With  xe  A  0.5  sec,  the  boundaries  appear  as  shown  in 
Fig.  7.  The  center  of  the  probable  crossover  region  is  oe^  ■  4  rad/seo. 
Just  as  in  the  Y0  ■  Kq  and  Y0  *  Kq/Jco  cases,  one  should  not  make  too 
much  of  the  exaot  numerical  value.  Instead,  the  important  estimates  are 
the  general  range  of  values,  l.e.,  2.5  £  £  5  rad/Bec,  and  the  order 

of  the  ccfc  values  for  the  various  controlled  elements,  i.e., 

>  ^Kq/J*  >  Hkc/(J<s)2*  Finally'  here  again  it  can  be  antici¬ 
pated  that  reductions  in  controlled  element  gain  will  be  offset  by 
increases  in  the  operator's  gain  and  vice  versa  1  and  that  a*,  regression 
will  occur  for  forcing  function  bandvidths  of  about  4  rad/sec. 
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Figura  7.  Crossover  Frequency  Estimation  for  Y0  ■  K0/(Joj)‘5 
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The  three  controlled  elements  considered  above  can  all  be  stabilised 
easily  by  the  human  without  any  fundamental  conflict  between  stability 
and  low  frequency  response.  The  operator  is  under  some  constraint  at 
frequencies  in  the  region  of  crossover,  but  elsewhere  his  transfer 
properties  can  conceivably  vary  a  great  deal  without  causing  difficulty. 

For  Instance,  at  very  low  frequencies  the  open-loop  describing  function  YpYc 
is  large  relative  to  one  for  all  three  controlled  elements,  so  over  the  low 
frequency  range  the  resulting  closed-loop  systems  all  exhibit  nearly  unity 
input— output  relationships;  the  effect  on  this  property  of  minor  reductions 
or  increases  in  the  low  frequency  portion  of  the  operator's  describing 
function  will  be  Insignificant  practically.  In  other  words,  the  systems 
treated  thus  far  have  not  been  equally  sensitive  over  the  entire  frequency 
range  of  interest  to  changes  in  the  operator's  characteristics,  but  have, 
Instead,  been  far  more  critically  affected  by  Yp  changes  In  some  ranges 
than  In  others.  One  could  conjecture  that  the  operator  might  be  able  to 
take  some  advantage  of  this  situation  by  being  relatively  haphazard  over 
frequency  regions  where  such  behavior  makes  little  difference,  while  being 
far  morn  precise  where  it  really  counts. 

There  are  controlled  elements  of  great  practical  importance  which 
tend  to  be  far  less  forgiving  of  variations  in  the  controller  character¬ 
istics.  The  most  common  of  these  are  vehicles  which  have  various  kinds 
of  instabilities  over  which  control  must  be  exerted  to  create  a  stable 
man-machine  system.  The  resulting  system  is  conditionally  stable,  because 
either  too  high  or  too  low  a  gain  will  result  in  Instability. 

|  Idealizations  of  two  controlled  elements  which  result  in  conditionally 

|  stable  olosed-loop  systems  have  been  selected  for  further  examination. 

I  These  are i 


*c  - 


The  first  can  be  thought  of  as  a  generalization  of  the  Yc  -  Kc/jm 
controlled  element  (-wherein  1/T  ■  o),  and  the  second  as  a  similar 
generalization  of  the  Y0  ■  Kq/ ( j<o)®  configuration.  Just  as  in  their 
special  l/T«*  0  cases,  the  major  difference  between  these  two  controlled 
elements  derives  from  the  necessity  of  using  low  frequency  lead  in  the 
second  to  attain  stability,  whereas  lead  is  not  required  for  stability 
in  the  first,  so  any  lead  generated  can  be  used  to  improve  the  high 
frequency  response. 


This  simple  divergence  controlled  element  has  been  studied  analytically 
previously  in  Ref.  1 ,  and  certain  facets  have  been  examined  experimentally 
in  Ref.  11 . 

Figure  8  presents  a  detailed  analysis  of  a  system  composed  of  the 
unstable  divergence  and  a  pure  time  delay,  t#.  As  might  be  expected, 
when  l/T  is  very  small  relative  to  1/t6  the  system  characteristics  appear 
very  Bimilar  to  those  for  Yc  ■  Kc/jm.  For  this  limiting  case  it  will  be 
recalled  that  the  lead  equalization  was  used  primarily  to  offset  high 
frequency  neuromuscular  syBtem  effects.  Consequently,  for  small  values 
of  l/T,  at  least,  the  operator  describing  function  for  control  of  this 
controlled  element  will  be 

Yp  i  Kpe-Ju/r«  (28) 

where  xe  ■  t  +  Tjj  -  Tl»  This  operator  describing  function  will  probably 
be  retained  as  l/T  increases  toward  1/t0  since  neither  lead  nor  lag 
modifications  are  oapable  of  improving  the  system  significantly.  As 
can  be  appreciated  from  Fig.  8,  the  value  of  1/t6  provides  an  upper 
limit  on  the  operator's  control  capacity,  so  Tp,  will  play  an  important 
role  indeed  in  reducing  *te  to  its  lowest  effective  value.  An  experi¬ 
mental  issue  of  great  interest  here  is  the  extent  to  which  the  operator 
may  possibly  be  able  to  reduce  the  neuromuscular  lag  and  reaction  time 
delay  when  the  pressure  to  do  so  is  extreme. 
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The  conditionally  stable  nature  of  this  system  Is  revealed  by  the 
low  and  high  gain  stability  limit  zero-db  lines  shown  in  Fig.  80.  The 
operator’s  gain  characteristics  must  be  such  as  to  establish  and  main¬ 
tain  the  actual  system  zero-db  line  somewhere  between  these  two  extremes. 
Deviation  from  this  policy  at  low  frequencies  results  in  an  aperiodic 
divergence,  whereas  at  high  frequencies  an  oscillation  is  the  penalty 
incurred  for  exceeding  the  stability  limit. 

The  tolerances  on  operator  gain  adjustments  can  be  made  extremely 
small  by  proper  adjustment  of  l/T.  As  the  divernence  time  constant 
decreases,  the  spread  between  l/T  and  l/re  also  decreases.  Maximum 
phase  margin  becomes  smaller  and,  most  important  of  all,  the  range  of 
stable  gain  levels  becomes  smaller.  This  range  is  theoretically  reduced 
to  a  single  value  of  gain  when  the  real  axis  breakaway  (o^)  coincides 
with  the  origin  on  the  root  locus.  The  position  of  the  breakaway  is 
easily  determined  since  it  coincides  with  the  point  of  zero  slope  on 
the  G(-cr)  amplitude  plot.  Thus,  the  breakaway  point  is  that  value  of  a 
for  (rhich  dO/do  is  zero.  Since  0(*-<j)  is 


G(-o) 


~(KpKcT) 
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To  +  1 


(29) 


then  the  slope  of  the  Siggy  plot  will  be 


dO(-a)  .  TTeQ(-°) 
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which  is  zero  when 
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The  breakaway  will  coincide  with  the  origin  when  the  divergence  time 
constant,  T,  is  equal  to  Ta.  Since  the  low  frequency  effects  of  the 
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neuromuscular  system  can  be  largely  canceled  by  the  lead  equalization, 
t6  Is  essentially  t.  Hence,  operator  control  of  divergences  having 
time  constants  of  the  order  of  0.15  sec  are  theoretically  possible. 

Early  experimental  evidence  on  the  control  of  such  divergences  tends 
to  support  this  conclusion  (see  Ref.  1,  p.  4j),  and  subsequent  experi¬ 
ments  (Ref.  11)  amount  to  conclusive  demonstration. 

For  the  experimental  program  oontenplated  here,  the  main  objective 
is  to  explore  the  effects  of  system  constraints  on  operator  variability. 
To  obtain  the  precision  desired  in  describing  function  measurements, 
relatively  long  tracking  runs  are  required,  some  with  large  forcing 
function  bandwidths.  Under  these  conditions,  divergence  time  con¬ 
stants  near  the  stability  limit  are  likely  to  be  too  extreme,  tending 
to  induce  fatigue  and  frustration  in  both  operator  and  experimenter. 
Much  smaller  values  of  l/T,  say  1  to  2  rad/sec,  impose  significant 
constraints  on  operator  variability  while  being  compatible  with  a  more 
relaxed  experimental  procedure.  For  example,  with  Te  -0.2  sec  and 
T  -  0.5  sec,  the  breakaway  frequency,  Is  3  rad/sec  and  the  total 

stable  gain  region  is  10  db. 

S.  Yq  ■  — 7— - — rr 

*(*  "  T) 

JuBt  as  for  its  double-integrator  special  case  (l/T  -  0)  this 
controlled  element  requires  low  frequency  lead  equalization  for 
stability.  Because  of  ‘the  divergence,  some  minimum  gain  level  is 
required,  thereby  leading  to  the  conditionally  stable  nature  of  this 
system.  The  appropriate  operator  describing  function  is 

Yp  i  K^e-J<J,T®  (TlJo>  +  1 )  (52) 

where  te  here  is  approximately  t  +  Tjj,  Just  as  for  Yc  -  Ko/(jci)2. 

Figure  9  presents  a  detailed  unified  servo  analysis  of  this  system. 
There  the  l/Tj,  lead  break  point  occurs  before  the  l/T  divergence. 

This  relative  position  is  appropriste  for  the  moderate  to  large  values 
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ft*1  l/T  vhioh  are  of  primary  interest  heirs ,  for  very  tumll  values  of  1  /T, 
however,  the  relative  position  of  these  break  points  will  be  reversed. 

That  this  relative  order  of  l/Tp,  and  l/T  is  pertinent  can  be  readily 
appreciated  by  Considering  the  alternative. 

In  the  limit,  as  l/Tj,  approaches  zero,  this  man-maohlne  system  tends 
to  approach  in  form  that  for  the  simpler  first-order  divergence  considered 
above.  In  such  a  limiting  circumstance  the  controllability  limit  would 
occur  when  T  «*  Te.  Because  the  lead  equalization  is  not,  in  thlB  case, 
tending  to  cancel  neuromuscular  lags,  Te  is  substantially  larger  than  that 
for  the  simple  divergence. 

Although  the  root  loci  for  the  simple  divergence  of  Fig.  8  differs 
drastically  from  that  Bhovm  in  Fig.  9,  the  essential  character  of  these 
systems  in  the  region  of  crossover  is  similar.  This  is  most  easily 
seen  by  comparing  Figs.  8c  and  9c.  The  differences  so  prominent  on  the 
root  loci  occur  at  gains  which  are  lower  than  those  of  direct  concern. 

The  key  similarity  in  both  circumstances  is  that  as  1 /T  is  increased, 
the  tolerances  on  available  gain  variations  decrease.  The  central  dif¬ 
ferences  between  the  two  constraining  situations,  already  mentioned  above, 
are  the  necessity  for  low  frequency  lead  generation  in  one  instance  and 
the  differences  in  Te. 

a.  vmcemjm  mbabunb  ahd  mxnxmxzaiior  mmsmm 

The  mean- squared  error  generated  in  a  man-machine  system  with 
stationary  forcing  functions  has  significance  from  at  least  three  points 
of  view.  First,  mean- squared  error  or  any  closely  related  error  func¬ 
tional,  such  as  fej ,  is  often  used  for  the  relative  assessment  of  com¬ 
peting  man-oaohine  systems*  Second,  the  stationarity  of  mean-squared 
error  values  gives  a  direct  and  easily  obtained  indication  of  over-all 
system  stationarity.  Learning  curves,  for  example,  usually  consist  of 
a  record  of  mean- squared  error  or  similar  quantity  versus  time,  with 
stationary  conditions  being  presumed  for  the  system  as  well  as  the 
performance  measure  when  it  has  achieved  what  appears  to  be  a  constant 
minimum  value.  Third,  there  is  some  reason  to  believe  that  the  operator's 
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describing  function  parameters  are  finally  adjusted  so  that  closed-loop 
low  frequency  performance  in  operating  on  the  forcing  function  minimizes 
or  nearly  minimizes  tbe  mean-squared  tracking  error.  This  is,  In  fact, 
Adjustment  Buie  4a. 

Because  of  the  importance  of  mean-squared  error,  it  was  desirable  to 
develop  some  simplified  relationships,  using  the  pilot  model,  that  can 
be  used  to  make  predictions  which  can  subsequently  bo  tested  by  experi¬ 
ment.  unfortunately,  not  enough  was  known  about  remnant  to  include  it 
in  any  mean-squared  error  formulations.  However,’  it  iB  not  difficult 
to  separate  the  mean-squared  error  into  two  components  —  that  due  to 
remnant  and  that  due  to  forcing  function.  Attention  here  will  be  devoted 
to  the  mean-squared  error  component  due  to  forcing  function,  i.e,,  the 
first  term  in  Eq  4. 


by 


The  mean-squared  error  component  due  to  forcing  function  is  given 


(35) 


where  «eftl  is  the  error  power  spectral  density  due  to  An  extremely 

simple  yet  remarkably  universal  open-loop  system  describing  function 
model  is  the  crossover  model  of  Eq  10,  which  is  repeated  below. 
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This  model  has  a  high  degree  of  validity  for  the  Y0  ■  Kq/Jcu  system  and 
is  also  quite  good  for  most  Y0  ■  Kq  systems.  It  can  have  some  validity 
for  special  cases  in  other  systems  also,  such  as  the  Y0  ■  Kc/( Joi  -  l/T) 
system  when  1 /T  is  small.  Although  this  open-loop  describing  function 
also  has  merit  as  an  approximation  to  describing  functions  for  Y0  ■  Ko/(jcu) 
in  the  crossover  region,  its  description  of  the  phase  at  frequencies  not 


too  far  below  Is  Inaccurate.  Still,  mean-squared  error  oalculatlons 
baaed  on  this  simple  form  will  nave  a  measure  of  applicability  to  many 
of  the  situations  to  be  examined  experimentally* 

The  simplified  systems  for  which  mean-squared  error  calculations 
will  be  made  are  shown  In  Fig.  10.  The  crossover  model  system,  Fig.  10a 
leads  to  Integrals  which  cannot  be  evaluated  without  making  simplifying 
approximations,  although  it  can  easily  be  evaluated  numerically.  Both 
of  these  approaches  will  be  used  below.  For  all  the  calculations  the 
forcing  function  spectral  density  will  be  taken  as  the  rectangular 
spectrum  given  by  Fig.  10b. 
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b)  Forcing  Function  Powor  -  Spectra!  Density 


Figure  10.  Elements  of  ef  Calculations 


The  error- to- input  describing  function  for  the  crossover  model  is 

given  by 
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When  Eq  34  ie  inserted  into  Eq  33 
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This  integral  is  not  readily  evaluated.  However,  if  sin  cwe  is  replaced 
by  its  argument  (a  good  approximation  for  Bystems  wherein  is  such  that 
«  l)  the  integral  is  easily  found.  The  result  is 


if  the  hyperbolic  arctangent  is  carried  only  to  the  second  term  in  its 
expansion.  The  exceptionally  simple  result  (Ref.  37)  given  by  Eq  37  is 
often  referred  to  as  the  "one-third  law."  It  will  later  be  seen  to  hold 
remarkably  well  as  long  as  o^/’o^  1b  somewhat  less  than  1 . 

Equation  33  can  also  be  integrated  numerically.  The  result  of  such 
an  operation  is  shown  in  Fig.  11.  In  this  plot  Te  serves  as  a  norna.1- 
izing  parameter,  so  the  figure  ie  essentially  relative  mean-squared 
error  versus  crossover  frequency,  with  forcing  function  bandwidth  as  a 
parameter.  These  results  are  of  very  great  Interest  indeed,  for  they 
imply  much  about  operator  adjustment!  It  will  be  noted  that  the  minimum 
values  of  mean-Bquared  error  for  the  smaller  cm  conditions  occur  at  the 
stability  limit.  This  implies  that  crossover  frequency  should  be  adjusted 
at  or  near  this  value  and  be  kept  there,  l.e.,  ■  'constant.  This  is  a 

highly  enlightening  analytical  Justification  for  Adjustment  Rules  3a 
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and  5b,  which  were  based  on  extrapolations  of  experimental  observations. 
Bui  Fig.  11  contains  still  wore  inform lion.  Notice  that  at  a  value  of 
lying  between  1 .2  and  1 .4  the  minimum  mean- squared  error  will  no 
longer  be  along  the  stability  limit.  Instead,  smaller  values  of  mean- 
squared  error  can  be  obtained  by  a  drastic  reduction  in  the  crossover 
frequency.  This  phenomenon  is,  of  course,  the  regression  noted  in 
Adjustment  Rule  5c.  Judging  from  Fig.  11,  such  regressions  should  occur 
when  cng./a^j  becomes  greater  than  about  0.8,  which  is  completely  compatible 
with  the  experimental  data  on  which  the  order  of  magnitude  condition, 
ajj/ocfc  ■  1 ,  was  based. 
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Figure  1 1 .  Mean-Squared  Error  Based  on  Crossover  Model 
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The  experimental  portions  of  this  study  were  Intended  to  make  explicit 
and  to  quantify  those  aspects  of  human  adaptive  behavior  which  can  be 
demonstrated  in  essentially  stationary  experimental  situations.  Great 
emphasis  was  placed  on  the  isolation  and,  wherever  feasible,  control  of 
sources  of  variability  so  that  the  experimental  findings  could  be  pre¬ 
sented,  insofar  as  possible,  in  a  deterministic  rather  than  statistical 
fashion.  This  was  accomplished  by  precise  measurements,  sometimes  on  two 
independent  analysis  machines j  by  using  specially1  oontrived  experimental 
and  measurement  procedures;  and  by  paying  meticulous  attention  to  detallB. 
These  means  have  been  successful  enough  to  allow  definitive  results  to  be 
obtained  using  a  restricted  number  of  experimental  runs. 

The  general  measurements  and  task  variables  involved  in  the 
experiments  are  summarized  In  Fig.  12.  The  task  variables  under  the 
experimenter's  control  are  enclosed  In  dashed  boxes;  those  actually 
controlled  In  this  experimental  series  were  the  forcing  function  and 
controlled  element.  For  the  data  to  be  discussed  here,  Yc  was  set 
equal  to  Kq,  Ke/joi,  Kq/ ( Jos ~1  /T) ,  Kc/(jco)2,  and  Ko/ju>(ja>“1 /T) .  The 
basis  for  the  selection, of  all  these  controlled  elements  has  been  laid 
in  the  previous  ohapter.  The  selection  of  i(t)  will  be  discussed 
subsequently  in  this  chapter. 

The  measurable  signals  in  the  control  loop  are  characterized  by 
power  speotra  such  as  ®aa,  mean-squared  values  such  as  e^,  and  ampli¬ 
tude  probability  distributions  such  as  pa.  The  operator  is  character¬ 
ized  by  the  describing  function  Yp  and  the  remnant  4^.  These  quantities 
are  sufficient  to  completely  define  the  man-machine  system  in  spectral, 
average,  and  amplitude  distribution  terms.  The  key  analyses  necessary 
are  those  pertinent  to  determination  of  power  spectra  and  cross  speotra. 


.General  Measurements  and  Task  Variables 


The  analytical  details  Involved  in  proceeding  from  the  spectral  quantities 
to  describing  function,  remnant,  etc.,  are  so  veil  known  and  documented 
(e.g.,  Refs.  13,  21,  24,  29,  34,  35,  50,  51 >  etc.)  as  to  make  detailed 
repetition  here  unnecessary.  However,  it  is  worth  mentioning  again  that 
the  correlation  coefficient,  p,  as  used  in  these  experiments  has  been 
made  especially  sensitive  to  operator  time  variations  by  virtue  of  the 
forcing  function  type  selected  and  the  analyzer  used.  Also,  another 
dimensionless  quantity,  has  been  defined  to  indicate  the  relative 
amount  of  remnant  in  the  total  output.  Finally,  because  both  periodic 
and  random  time  functions  are  present  in  the  experimental  systems,  the 
power  spectra  may  have  two  different  natures,  line  and  continuous. 
Terminology  and  definitions  which  are  consistent  with  both  possibilities 
are,  therefore,  needed.  These  are  given  below. 

In  general  the  symbol  ^(to)  will  be  used  to  denote  the  power  spectral 
density  of  a  time  signal  x(t) .  In  more  compact  style  the  power  spectral 
density  is  also  referred  to  as  the  power  spectrum  or,  even  more  simply, 
as  the  spectrum.  The  autocorrelation  function,  Rxx(t) ,  of  the  same  signal 
is  defined  to  be 


/•T 

(t)  -  Lim  «  f  x(t)x(t  +  x)dt 


-  x(t)x(t  +  x)  (38) 

The  power  speotral  density  is  proportional  to  the  Fourier  cosine  trans¬ 

form  of  the  autocorrelation  function,  R^t), 

■  **■  f  KxxC'O  C0B  CtfrdT  (39) 

Jo 

The  mean- squared  value  of  the  time  signal  will  be 

_  .  /»« 

x2  ■  ®xx(w) d* 

.  R^O)  (40) 


If  x( "t)  is  si  periodic  func wicn f 


E  <Px(%)  Bln  (%t  *  tB) 

m»1 


then  the  autocorrelation  function  is 

_  ,  ,  &  *4(®ta) 

Rxx(t)  -  E  — s —  c°3  %T 

m»1 


and  the  power  spectral  density,  considering  only  positive  frequencies, 
will  he 


fc^co)  ■  *  £  (^(<t^)6(o  -  cj^) 
m«1 


The  magnitudes  of  the  delta  function  spectral  lines  in  are  it  timeB 
the  square  of  the  peak  amplitudes  of  the  sinusoidal  components  in  x(t) . 
When  plots  r  presenting  the  line  spectra  are  made,  the  actual  quantities 
used  are  <p^(%),  expressed  in  power  db,  i.e.,  10  log10  4^-  1516  mean 

square  of  x(t),  following  Eq  40,  is 


if*  M  - 

7Z  I  *  E  -  %)  to 

wq  m«1 


* 1  *SK> 

m-1 


When  x(t)  Is  a  stationary  random  function  the  autocorrelation  func¬ 
tion  definition  resales  unchanged,  and  the  power  spectral  density 
definition  based  on  the  autocorrelation,  Eq  59,  and  the  mean-squared 
value  expression,  Eq  40,  are  also  unmodified.  The  form  of  the  power 
spectral  density  will,  however,  be  entirely  different  from  the  sum  of 
delta  functions  occurring  for  the  periodic  case,  i.e.,  the  power  spectral 
density  will  be  continuous.  If,  for  one  reason  or  another,  the  component 
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of  the  power  spectral  density  which  is  due  to  a  random  component  in  x(t) 
is  required  to  ha  specifically  identified,  then  the  notation  tp^  is  used. 
Thus,  if  both  periodic  and  random  components  appear  in  x(t),  the  power 
spectral  density  will  be 

M 

*xtf(“)  ■  +  *  E  q^(a^n)0<as  -  %)  (45) 

m-1 

Power  db  are  also  used  for  plots  of  cp^cc),  !•«•>  10  logiQ  ^(cn). 

The  experimental  measurement  of  the  quantities  summarized  above  is 
critically  dependent  on  the  manner  in  which  the  measurements  are  executed 
and  on  the  apparatus  associated  with  the  control  tasks  and  signal  analysis. 
This  chapter  Is  therefore  devoted  to  descriptions  of  the  physical  layout 
and  equipment,  the  task  variables  (forcing  function  and  controlled  element), 
the  measurement  situations  and  data  presentation  conventions,  and  other 
empirical  aspects.  Particular  attention  is  given  to  the  description  of 
special  procedures  evolved  in  the  course  of  the  program  to  avoid  pitfalls 
or  to  Increase  the  reliability  of  the  data. 

A.  mfZOAL  XAYOlff  AMP  IQKIMI1R 

The  experiments  were  performed  in  a  laboratory  area  aoneiBting  of 
two  connected  rooms.  The  larger  of  the  two  rooms  contains  all  of  the 
electronic  equipment  for  performing  and  analysing  the  experiments.  The 
smaller  room  contains  the  manipulator  and  display.  In  this  way,  the 
operator  is  isolated  from  the  measuring  equipment  and  othar  disturbances. 

An  intercom  is  used  for  communication  between  operator  and  experimenter. 

The  experimental  equipment  consists  mainly  of  three  basic  systems. 

The  first  is  the  watt-hour  meter  analyser  shown  in  Figs.  15  and  14.  This  is 
an  analog  device  capable  of  computing  the  cross  spectra  between  the  forc¬ 
ing  function  and  two  other  functions  at  eaoh  of  the  ten  frequencies 
contained  in  the  input,  in  real  time..  This  device  is  more  extensive 
than,  but  is  otherwise  similar  to,  apparatus  originally  used  by  Russell 
(Ref.  47)  to  obtain  human  response  measurements.  It  contains  ten  sinu¬ 
soidal  function  generators,  Fig.  15,  from  which  the  forolng  function  is 
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synthesized.  Since  a  separate  signal  is  available  at  each  of  these 
frequencies  (both  sine  and  ooslne  components) ,  the  calculation  of  the 
cross  spectrum  between  the  forcing  funotlon  and  another  signal  is  rela¬ 
tively  single  (see  Refs.  44  and  21).  At  each  of  the  input  frequencies, 
a  multiplier  and  an  integrator  are  required  for  the  real  and  imaginary 
parts.  These  functions  are  provided  by  watt-hour  meters,  Fig.  l4.  The 
use  of  these  devices  makes  possible  simultaneous  computation  of  cross 
spectra  at  all  input  frequencies.  The  limitations  of  this  equipment 
are  that  measurements  are  restricted  to  input  frequencies  and  it  is  not 
possible  to  measure  power  spectra.  The  major  use  of  this  equipment  is 
to  calculate  the  describing  function  during  the  actual  experiment,  so 
that  the  results  are  immediately  available.  Mean-squared  error,  e^, 
and  operator  output,  c2,  are  also  measured  by  other  vatt-hour  meters  in 
the  same  gear. 

The  seaond  system  Is  the  magnetic  tape  recording  equipment,  Fig.  12. 
During  am  experiment  the  forcing  function,  i(t) j  error  signal  (input  to 
operator),  e(t);  operator's  output,  c(t);  and  controlled  element  output 
(system  output),  m(t),  are  reaorded  on  tape.  These  signals  are  later 
transcribed  onto  a  continuous  tape  loop  which  provides  repetitive  data 
for  further  processing. 

The  third  basia  component  of  the  data  reduction  equipment  1b  the 
speotral  and  aross-epeatral  analyzer  (Ref.  27) •  Besides  providing  an 
alternative  means  for  cross-spectrum  measurement,  this  equipment  (Fig.  1 6) 
fills  the  need  for  the  additional  measurements  of  power  spectra  needed  ( 
for  remnant  calculations.  It  is  capable  of  computing  the  cross  spectrum 
between  any  two  signals  without  necessarily  requiring  the  forcing  func¬ 
tion  to  be  mads  up  of  sinusoids.  However,  because  of  its  sequential 
operation,  the  input  data  must  be  repeated  once  for  each  frequency  for 
which  a  value  is  desired,  so  that  on-line  operation  is  not  possible. 

Input  data  for  this  analyzer  are  obtained  from  the  magnetite  tape 
recordings.  From  these  tapes  are  obtained  the  correlation  coefficient, 
the  describing  funotlon  calculations  for  corroboration  of  the  vatt-hour 
meter  analyzer  results,  and  remnant  spectra  as  print-outs  on  a  paper  tape. 


FIGURE  15-  TAPE  RECORDING  EQUIPMENT  FIGURE  16.  CROSS  SPECTRAL  ANALY2ER 


The  controlled  element  is  generated  with  an  analog  circuit.  This 
circuitry  Is  contained  in  a  separate  rack  (Fig.  17)*  By  switching 
feedback  components,  a  variety  of  response  characteristics  can  be 
obtained.  These  include  the  Halting  oases  used  for  the  majority  of 
’  experiments,  Kq»  Ka/jtOj  and  Xa/(J«)a>  and  the  unstable  oases  which 
approach  Kq/Joj  and  Ko/(J«)2  as  l/T  approaches'  zero,  ^/(Jco  -  l/T)  and 
Kq/JcjCJo-  l/T). 

The  manipulator  used  for  these  experiments  (Fig.  1 8)  was  designed 
to  exhibit  the  minimum  possible  Inertia  and  damping,  so  that  it  behaved 
essentially  as  a  spring  restraint.  Only  the  lateral  degree  of  freedom 
(for  transverse  or  roll  control)  was  used  in  this  experimental  series. 
The  stick  was  manipulated  to  minimize  the  error,  displayed  on  a  6" 
oscilloscope  face  as  the  horizontal  distance  of  a  spot  from  the  center. 

In  order  to  estimate  the  accuracy  with  which  describing  functions 
can  be  measured  using  this  equipment,  an  extensive  aeries  of  teste  were 
performed  (Refs.  44  -  46) .  The  results  of  these  teste  indicate  an 
accuracy  for  both  the  watt-hour  meter  analyzer  and  the  spectral  and 
eross-apeotral  analyser  of  about  O.J  db  in  amplitude  and  4°  in  phase 
angle  over  a  dynamic  range  of  40  db  in  amplitude  and  a  frequency  range 
of  2  decades  (0.14  to  l4  rad/sac) .  Typical  measurements  are  shown  in 
Fig.  19.  For  these,  the  tracking  loop  contained  only  the  controlled 
element,  whloh  was  set  to  give  the  indicated  transfer  functions. 
Excellent  agreement  was  obtainad  between  the  two  analyzers  and  the 
theoretloel  response  for  the  known  controlled  elements. 

b.  nscouKXOM  or  m  maaa  nvoixov 

Since  the  watt-hour  meter  analyzer  can  measure  only  at  the  forolng 
function  frequencies,  the  selection  of  these  frequencies  and  the  forcing 
function  spectral  shape  is  of  great  importance  in  minimizing  measurement 
variability.  In  Refs.  21,  40,  4l,  and  54  conditions  pertinent  to  the 
use  of  deterministic  signals  to  simulate  random-appearing  signals  having 
Gaussian  amplitude  distributions  are  presented.  Requirements  for 
independence  of  sl&e  waves  and  appropriate  choice  of  run  lengths  for 
minimizing  the  measurement  error  are  diaoussed  in  detail.  In  Ref.  54 
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Figure  19.  Aaouracy  and  Compatibility  of  Analyzers 


it  is  shown  that  a  close  approximation  to  a  Gaussian  amplitude  distri¬ 
bution  oan  be  achieved  when  as  few  as  five  properly  chosen  independent 
sine  waves  are  summed •  Therefore  the  ten  Independent  signals  available 
are  more  than  adequate ,  and  even  allow  some  flexibility. 

The  desirable  characteristics  of  a  forcing  function  composed  of 
independent  sine  waves,  and  to  be  uaed  for  human  response  studies,  can 
be  specified  as< 

1.  Random-appearing,  so  that  the  operator  cannot  detect 
any  internal  ooherence  in  the  forcing  function  and 
thereby  adopt  a  higher  level  of  behavior 
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2.  Having  a  Gaussian  amplitude  distribution,  bo  that 
Oaussian-lnput  describing  function  theory  may  be 

used 

J.  Extending  over  about  2  decades,  including  loir 
frequencies  below  0.5  rad/sco 

4.  Having  a  sinply  definable  cutoff  frequency 

5.  Providing  energy  In  the  crossover  mage  without 
disturbing  the  operator's  low  frequency  performance 

6.  Being  composed  of  low  frequency  sine  waves  which 
are  Integral  multiples  of  the  run  length  to  minimize 
averaging  error 

The  forcing  function,  l(t),  selected  for  the  set  of  experiments 
described  in  this  report  —  the  so-called  augmented  rectangular  input 
spectrum — conforms  to  these  requirements.  This  input  is  designated 
a>L,  cj.  where  ajj.  is  the  cutoff  frequency  in  radlans/second  and  id 
the  rms  amplitude  of  this  input  as  measured  on  the  pilot's  display. 
Ordinarily,  cj.  was  1/2",  so  is  specified  in  our  notation  only  when 
it  differs  from  l/2".  The  frequency  setting,  <%,  and  the  number  of 
periods  for  each  component  in  the  fixed  240-sec  run  length,  Tr,  are 
as  follows  t 


% 

n  -  Tro* 

0.157 

6 

0.262 

10 

0.393 

15 

0.602 

23 

O.969 

37 

1.49 

37 

2.34 

97 

4.03 

154 

7.57 

289 

13.8 

527 

Three  cu&  values  were  used — 1 .5,  2.5,  and  4.0  rad/sec.  For  the  mj., 

1/2"  forcing  funotlon,  alone  characterizes  i(t).  To  define  these 
three  inputs,  the  amplitudes  at  the  lowest  six,  seven,  or  sight  frequencies 
were  set  equal,  for  cutoff  frequencies  of  1  »5j  2.5*  or  4.0  rad/sea, 
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respectively.  The  amplitudes  of  the  remaining  frequencies  were  set  to 
one  tenth  of  the  low  frequency  amplitudes  (20  db  attentuatlon) .  The 
logarithmic  spacing  (which  facilitates  fitting  theoretical  curves  to 
the  measured  values)  was  selected  to  insure  the  effective  independence 
of  the  sine  wave  components  over  the  run  length  of  interest. 

Figure  2o  shows  measurements  of  the  three  forcing  functions  which 
Illustrate  the  augmentation  of  the  well-defined  rectangular  cutoff 
frequency  spectral  shape  by  a  series  of  very  low  energy  signals  extend¬ 
ing  up  to  13.8  rad/sec.  Specific  tests  were  performed  to  assure  that 


0.1  1.0  10.0 


w(  rad /sec) 

Figure  20-  Measured  Input  Power  Spectra  Magnitudes 

these  forcing  function  spectra  were  indeed  nearly  Gaussian.  Figure  21, 
vhloh  is  representative  of  the  data  for  o&  ■  1>3  and  ay,  m  ^»0  as  well, 
la  a  plot  on  probability. paper  of  the  cumulative  probability  distribu¬ 
tions  for  seven  a*  ■  2.3  inputs.  A  straight  line  is  produced  by  a 
Gaussian  distribution.  Amplitude  distribution  analyses  were  also  made 
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Probability  Distribution  (for  w,  =  2.5 i 


for  these  three  inputs  using  a  chi-squared  test  (e.g.,  Ref.  23). 
These  indicate  that  the  inputs  are  all  Gaussian  at  the  0.05  level. 


a 


Figure  21  .  Measured  Cumulative  Probability  Distribution 
of  Forcing  Function  Amplitude 

For  each  of  the  forcing  functions,  the  chi-squared  test  for  goodness  of 
fit  to  a  Gaussian  distribution  used,  as  a  null  hypothesis,  the  hypothesis 
that  any  difference  between  the  observed  distributions  and  the  Gaussian 
distribution  was  the  result  of  ohanoe  forces.  For  this  test, 
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where  0*  represent*  the  observed  measurements  in  a  given  category,  which 
is  a  decile  here,  and  Ai  1b  the  number  of  observations  expected  In  a 
decile  if  the  distribution  Is  Gaussian  with  the  measured  mean  and  vari¬ 
ance  .  Because  the  categories  used  were  deciles,  k  »  10-  Slnoe  the 
Gaussian  distribution  has  two  parameters  and  these  parameters  are 
unknowns  in  thiB  case,  there  are  seven  degrees  of  freedom.  Using  this 
hypothesis  as  the  test  hypothesis,  any  difference  between  the  observed 
distributions  and  the  Gaussian  distribution  which  turned  out  to  be 
significant,  l.e.,  not  Gaussian  at,  say,  the  0.05  level  of  significance, 
could  be  expected  to  occur  only  one  time  in  twenty  (5  percent)  or  1>jbb 
if  solely  chance  forces  caused  that  difference.  All  of  the  chi-squared 
tests  for  individual  inputs,  i.e.,  five  1 .5,  seven  2.5,  and  three  4.0 
forcing  function  samples,  failed  to  reject  the  null  hypothesis.  The 
chi-squared  values  were  then  averaged,  and  theBe  averages  are  plotted, 
together  with  the  ranges,  in  Fig.  22. 

The  augmented  rectangular  forcing  function  shape  was,  of  course, 
Intended  to  overcome  the  past  deficiencies  In  human  response  data  due 
to  their  limited  frequency  range  and  to  the  lack  of  precise  measure¬ 
ments  in  the  critical  crossover  region 1  Hopefully,  the  tiny  high 
frequency  forcing  function  components  greater  than  cuj.  do  little  to 
disturb  the  operator's  low  frequency  performance.  Ideally,  the  low 
frequency  characteristics  should  be  substantially  the  same  as  those 
measured  with  a  pure  rectangular  forcing  function  spectra,  l,e, ,  without 
the  augmentation  frequencies. 

In  order  to  demonstrate  the  relative  lack  of  influence  of  the  low 
energy  higher  frequency  Blgnals  on  the  low  frequency  describing  func¬ 
tion  behavior,  we  refer  to  Elklnd's  definitive  set  of  data  for  Y0  -  1 
(Ref.  IJ).  Figure  2J  shows  Yp  for  Elklnd's  R,40  forcing  function — 
t.  pure  rectangular  forcing  function  spectra  consisting  of  40  equally 
spaced  equal-amplitude  sine  waves  up  to  -  2,5  rad/sec  and  of 
1.0  Inches  rms  magnitude  —  and  Yp  for  Elklnd's  B-6  forcing  function 
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Figure  22.  Average  Value#  for  Augmented  1 .5,  2,5,  and  4.0 


which  was  an  augmented  rectangular  forcing  function  spectra  cor-r 
144  sinueoida  and  eq.uivale.nt  to  q  ■  J.O,  1"  in  our  notation. 
be  noted  that  the  B-6  data  provide  a  convincing  extrapolate  u-;> 

R.40  data  to  the  croaaover  frequency.  Subsequently  we  shall  i 

that  the  Y0  ■  Kq  data  developed  in  the  current  program  are  con^: 

with  Bikind 1 s  data.  It  therefore  follows  that  the  augmented  r<  >  jolar 

forcing  functlone  used  here  provide  a  reasonable  extrapolation  low 

frequency  measurements  to  and  beyond  the  crossover  frequency  and  that 

the  high  frequency  shelf  does  not  materially  affect  the  operator's 

behavior. 
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Figure  23.  Comparison  of  Elkind '  s  Data 
for  Rectangular  (B.4o)  aM.Augmented  (U-6)  Forolng  Functions 
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The  measurement  situation  was  as  follows.  The  pllpt  manipulated  a 
stiok  (Pig.  24)  which  produced,  as  an  electrical  output  signal,  the  pilot's 
output,  o(t),  which  was  fad  to  the  controlled  element,  Y0  (refer  to  Fig.  12), 
The  damping  and  Inertia  of  this  stiok  were  made  as  low  as  possible . 

The  circuit  gains  were  set  to  give  the  following  deflection  sensitiv¬ 
ities  on  the  display  for  a  pure  gain  controlled  element  of  unity: 

Kb  ■  1  inch  (display) /6  deg  (stick) 

Fe  m  2.21  oz/deg  stick  (applied  with  a  4  inch  moment  arm) 


With  the  4  inch  moment  arm,  lateral  motion  of  the  operator's  hand  amounts 
to  about  0.07  inohes  (stick)  per  degree  of  stick  rotation.  Accordingly, 
the  sensitivity  can  be  expressed  in  terms  of  the  linear  motion  of  the 
operator's  hand  by  dividing  the  angular  sensitivity  (in  inches  per  radian) 
by  the  moment  arm,  i.e., 

Ks  ■  1  lnoh/6  deg  x  57*3  deg/ rad  x  0.25  inoh""1  -  2.50  inches/inch 

For  a  pure  gain  controlled  element  with  gains  other  tban  unity,  the  actual 
displacement  on  the  display  as  a  result  of  a  stiok  displacement  was  the 
value  given  above,  Kg,  multiplied  by  the  controlled  element  gain,  Kq.  In 
summary,  for  pure  gain  controlled  elements: 


Kg  «  1 t  K, 
or  Kb 

Kf 


0.167  inches/deg 

2.56  inohee/inoh  (4  inoh  moment  am) 
0.0754  inches/oz 


Kc  1 «  Ka 
or  Ka 


■  0.l67Ko  inohes/deg 

■  2.58  K0  inches/ inoh 


Kf  ■  0.0754  K0  inches/oz 


For  frequency-dependent  controlled  elements  the  dynamioe  of  the  con¬ 
trolled  element  intervene  between  the  stick  output  and  the  display.  Tha 
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deflection  sensitivities,  when  referred  to  the  display  as  a  commpon  point, 
then  become  dynamic  entities.  So,  for  a  unit  step  displacement  of  1  deg 
allied  to  the  Stic*  tne  deriection  on  the  scope  will  be  0.167  oC  '  [X0/‘J 
if  this  convention  Is  used. 

For  convenience  and  simplicity  it  is  much  easier  to  think  of  the 
deflection  sensitivities  in  terms  of  the  steady -state  characteristics 
which  are  ultimately  approached  after  a  unit  stick  deflection,  for  this 
type  of  description  the  appropriate  expressions  are 

Ks  ■  O.I67  [s^Xc(8)]  B*o  inc4®B/aeS 

Kb  ■  2.33  b\(s)]Bm0  inches/inch 

Kf  ■  0.0751*  fsl%(s)]B»o  inehes/oz 

where  k  is  the  order  of  the  free  s  term  in  Y0 

For  example,  If  Yc  ■  Kq/b  the  deflection  sensitivity  will  be  0.1 67  K,, 
iaohes/see  (display) /deg  (stick).  The  expressions  above  are  adequate  to 
characterize  system  gains  even  when  the  controlled  element  is  unstable, 
for  although  a  "steady- state"  never  ooours  there  la  still  a  definable 
"steady-state"  component  in  the  response. 

Further  conversion  factors  can  bs  defined  based  on  the  operator's  eye 
to  scope  distance  of  about  29  inches. 

Roll  axis  control  using  the  stiak  was  accomplished  by  grasping  the 
control  handle  and  exerting  lateral  foroe  on  it,  as  illustrated  in  Fig.  24. 

The  output  of  the  oontrolled  element  (system  output),  m,  was  subtracted 
from  the  input  forolng  function  to  produce  the  error  signal,  e.  This 
function  was  displayed  to  the  operator  on  the  oscilloscope*  Because  the 
watt-hour  meters  must  opsrate  at  60  ops,  part  of  the  oontrol  loop  used 
a  60  ops  carrier.  The  signals  from  the  d-o  portion  of  the  loop  were  fed 
to  the  tape  recorder  for  future  processing*  The  error  and  pilot's  out¬ 
put  signals  were  passed  through  pre-emphasis  networks  before  recording 
in  order  to  increase  the  slgnal-to-nolse  ratio  on  the  tape  at  high 
frequencies. 

The  computation  time  for  the  experiment  was  accurately  oontrolled, 
equal  to  an  Integral  number  of  period*  of  the  input  frequencies.  Prior 


77 


Pott  A  tit  — 
Pattering 
Foret  Spring 


intercom 


Vi  scat  Display 


-  Extreme  Stick 
Positions 


Pivot 

Point 


Pitch  Axis 
damp 


!  w 

.  i"f  mL. 


Pitch  Axis 
Pivot 


1/  ■■■<*! 


Peit  Axis 
Pivot  Point 


Pitch  Axis 
Locking  Bar 


x  Arm  Pest 
Pivot  Point 


Pitch  Axis 
Restoring  Spring 


Not*: 

Pitch  OKI*  locked  during 
roll  exit  control  trlolo 


FIGURE  24.  STICK  MANIPULATOR 


to  the  beginning  of  the  computation,  10  to  15  sea  wers  allowed  for  the 
operator  to  reach  a  eta Die  tracking  condition.  Before  and  after  each 
experiment.,  the  readings  of  the  watt -hour  matara  w»r»  recorded.  Proa 
these,  the  operator's  describing  function  was  calculated  at  each  of  the 
input  frequencies. 

The  Yp  data  have  generally  been  presented  as  dimensionless  quantities. 
To  do  this,  the  describing  function  has  bsen  defined  in  terms  of  the 
pilot 'e  output  and  error  signals,  referred  to  deflection  on  the  display 
(in  inches) .  Similarly,  the  power  spectral  densities  are  referred  to 
the  display  and  are  expressed  in  units  of  square  inches  per  radian  per 
second. 


In  cases  where  the  power  spectral  density  has  both  a  continuous 
portion  and  discrete  lines,  the  line  components  are  plotted  ae  If  the 
total  power  in  the  line  were  distributed  over  the  bandwidth  of  the  meas¬ 
uring  equipment.  In  this  way  the  integrated  speotrum  still  gives  the 
total  power,  but  the  peak  amplitude  shown  for  the  line  has  no  significance. 

s.  mmmm  nomw 
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In  processing  the  data,  cross  spectra  have  bean  cosputed  between  the 
input  function  and  the  signal  to  be  analysed.  Prom  these  it  Is  possible 
to  oanpute  the  describing  function t 
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Using  this  technique,  the  ef feats  of  signal  coaponsnts  at  othsr  frequen¬ 
cies  or  of  a  random  nature  are  minimised,  since  these  are  unoorrelated 
with  the  Input  and  tend  to  average  out.  It  is  possible  that  the  describ¬ 
ing  function  may  not  be  stationary  but  may  exhibit  a  low  frequency  time 
dependenos.  In  these  oases  the  computed  value  of  Yp  is  an  average  value 
for  the  length  of  the  experiment. 

A  slight  modification  has  bean  mads  to  the  Yp  calculation  for  the 
results  from  the  spectral  and  oross-spectral  analyser.  Since  the  fre¬ 
quency  at  which  the  measurement  is  mads  is  set  independently  of  the  input 
funotion  generator,  a  slight  error  In  the  setting  will  affeot  the  results. 
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To  eliminate  this,  the  ratio  of  the  cross  spectrum  to  the  input  power 
spectrum  vac  calculated  for  each  signal,  and  these  r»tiu»  ««n  then  used 
to  calculate  Yps 


v  *ia(<%)/*li(<»fa) 
p  "  *ieK>/*liW 


(W) 


The  bulk  of  the  describing  function  data  presented  in  this  report  has 
been  calculated  from  the  watt-hour  meter  results.  In  certain  cases  the 
spectral  and  cross-spectral  analyzer  results  have  been  used;  mainly  where 
the  watt-hour  meter  results  were  inaccurate  due  to  small  signal  levels. 
Generally,  the  describing  function  was  calculated  from  watt-hour  meter 
measurements  of  the  error  signal,  e,  and  pilot's  output,  c.  However, 
for  many  of  the  controlled  elements  that  were  used,  this  results  in  rela¬ 
tively  inaccurate  data  at  the  lowest  frequencies.  For  Yc  ■  Kq/Jo  or  K0/(Jcu)2 
the  controlled  element  gain  is  very  large  at  low  frequencies  and  very 
small  at  high  frequencies.  As  a  result,  the  operator  is  required  to 
produoe  only  email  aontrol  dofleotions  at  low  frequencies  to  obtain 
large  system  outputs,  m(t).  However,  at  high  frequencies  he  must  pro- 
duoe  large  output  deflections  in  order  to  obtain  reasonable  system 
outputs.  The  power  spectrum  of  the  pilot's  output,  «co,  then  increases 
with  frequenoy,  with  the  major  portion  of  hie  total  output  appearing  near 
the  input  cutoff  frequenoy.  Conversely,  the  syetem  output,  will 
be  relatively  flat  at  low  frequencies  and  will  decrease  rapidly  above 
crossover.  Since  it  is  neoeasary  to  adjust  the  analyser  gains  to  avoid 
saturation  due  to  the  large  frequency  components  present,  the  very  small 
components  at  low  frtqutncles  in  ths  pilot's  output  do  not  product 
sufficient  deflections  on  the  watt-hour  meters  to  obtain  acourat# 
results.  However,  if  ths  system  output  is  used,  correspondingly  poor 
data  will  be  obtained  at  high  frequencies.  For  most  of  ths  measurements 
it  was  desired  to  obtain  the  best  possible  results  near  oroseover,  and 
for  this  reason  ths  pilot's  output  was  used  in  the  analyslo. 

In  oertain  oases,  however,  to  dear  up  questions  wbieh  arose  about 
the  low  frequency  behavior,  the  analysia  was  based  on  the  system  output. 

In  particular,  there  was  a  strong  indication  that  large  phase  delays 
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exist  at  lew  frequencies.  These  have  been  noted  before  and  generally 
ignored  (Ref.  34,  Fig.  4-24C,  for  example).  Measurement  of  resulted 
In  too  loir  an  accuracy  to  determine  If  these  delays  actually  existed  or 
if  they  were  due  to  me*mir»m»nt  errors.  The  use  of  4^,  on  the  other 
hand,  provided  sufficient  accuracy  to  verify  their  existence. 

A  similar  problem  was  encountered  when  measuring  the  remnant  with 
the  spectral  and  cross-spectral  analyser.  In  this  case,  the  existence  of 
tape  background  noise  and  extraneous  noise  (such  as  that  from  the  play¬ 
back  amplifiers)  was  sufficiently  large  to  causa  large  errors  at  low 
frequencies  for  4ao  and  at  high  frequencies  for  9^.  Therefore,  the 
measurements  of  4^  and  p  were  based  on  4^  at  low  frequencies  and  4ec 
at  high  frequencies. 

The  Inputs  used  made  it  necessary  to  include  pre-emphasis  of  the 
error  and  pilot's  output  Bignals  before  recording  to  Increase  the 
amplitude  of  the  high  frequency  components  (above  4-6  rad/sea)  by  20  db, 
Without  it,  the  tape  noise  made  it  Impossible  to  obtain  good  acoumoy 
at  Input  frequencies  above  cutoff,  due  to  the  20-db  drop  In  input  aspli- 
tudes  above  this  frequency.  In  addition,  the  pre-aqphaols  also  Increased 
the  aooumoy  of  the  remnant  measurements  at  high  frequencies,  since  it 
effectively  increased  the  signal-to-nolee  ratio  at  the  Input  to  the 
spectral  and  Qross-spectral  analyser. 

A  large  portion  of  the  measurement  program  was  based  on  calculations  of 
the  correlation  coefficient,  p,  whioh  can  be  measured  with  the  speatral  and 
aross-spaatral  analyser.  This  value  provides  a  means  for  estimating  the 
degree  of  time  variation  in  the  operator* ■  transfer  funotion.  However,  the 
bandwidth  of  the  measuring  equipment  must  be  considered  in  determining  the 
prtolee  meaning  of  a  particular  p  value.  For  the  spectral  and  cross- spectral 
analyser  two  different  bandwidths  were  used,  depending  on  the  spacing  of  the 
input  frequencies  and  the  stability  of  the  equipment.  Generally,  for  low 
frequencies  the  narrow  bandwidth  (&pj  *  0.1 4l  rad/seo  between  3-db  points)  was 
used,  while  for  high  frequencies  the  wide  bandwidth  (dug  ■  0.626  rad/sec)  was 
used.  Since  the  pilot's  or  system's  output  is  approxlmatsly  the  sum  of  a 
line  spectrum  (linearly  correlated  with  the  input)  and  a  random  continuous 
speetrum,  the  correlated  power  measured  by  the  speatral  and  cross-spectral 
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analyser  Is  Independent  of  bandwidth,  while  the  uncorrelatad  power  is 
roughly  proportional  to  bandwidth*  Thus,  the  measured  value  of  p  for 
high  frequencies  (where  the  wide  bandwidth  is  used)  would  appear  to  be 
smaller  than  those  for  low  frequencies  (where  the  narrow  bandwidth  la 
used) .  To  avoid  this  discrepancy  the  high  frequency  values  were  cor¬ 
rected  to  be  equal  to  the  values  that  would  have  been  measured  if  the 
narrow  bandwidth  had  been  used,  even  though  they  were  actually  measured 
with  the  vide  bandwidth.  To  make  this  correction,  the  signal-to-nolse 
ratio  (correlated  to  uncorrelatad  power)  at  a  frequency  c%  can  be  expressed 


(1  “  P?) 


0  -  p|) 


*<PoK) 


rrr  for  the  narrow  bandwidth 


nq|(%) 

9nnK)4“e 


for  the  wide  bandwidth 


However,  the  ratio  <p§Afon  is  independent  of  the  measuring  equipment, 
and  oan  be  eliminated  between  the  two  equations .  Ibis  gives 

<nr^>  ■  trr$  <- 

or  by  rearranging 
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This  equation  oan  be  used  to  calculate  the  value  of  p  that  would  be 
measured  with  the  narrow  bandwidth,  Aoj ,  from  the  value,  pg,  measured 
with  the  wide  bandwidth,  As  a  check,  several  values  were  aotually 

measured  with  both  bandwldths.  The  results  wsre  as  follows; 


This  second  method  is  not  convenient  for  calculating  Yp,  but  it  serves 
as  a  computational  chest.  The  values  of  Yp  and  YpYu  are  found  from  t-h* 
first  method  of  computing  Yp  and  the  theoretical  value  of  Y,j.  Using 
confuted  values  of  p  and  YpYc  and  Eq  ?6,  the  ratio  ®#e/$ee  0S,n  be  calcu¬ 
lated  from 


This  ratio  can  also  be  measured ,  and  should  agree  with  the  calculated 
valu<' ,  Extensive  measurements  of  thiB  type  were  sale  for  one  of  the 
runs  using  -  1 .5,  l/4"  with  a  Ye  ■  Kc/(j<s)2.  The  results  are  shewn 
in  Figs.  25  and  26. 
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Figure  2%  Internal  Consistency  Check,  Yc  Auylltude  and  Phase 
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This  chapter  presents  the  experimental  describing  function  data 
obtained  to  serve  as  a  data  base  for  the  validation  and  elaboration  of 
the  analytical-verbal  model.  The  data  will  be  presented  in  various 
aggregations  designed  to  Illustrate  and  illuminate  the  findings.  Each 
set  of  data  taken  mo  an  element  in  a  grand  design  conceived  to  fulfill 
the  experimental  objectives  discussed  in  Chapter  II.  These  objectives 
were  sharpened  and  made  concrete  as  specific  experimental  conditions  by 
the  pre-experimental  analyses  of  Chapter  III.  The  resulting  experimental 
matrix  ia  shown  in  Fig.  27a. 

In  the  experiment  a  total  of  nine  subjects  were  used.  It  is  obvious 
that  filling  the  boxes  in  Fig.  27a  for  Kq  variation,  01  variation, 

amplitude  variation,  and  degree  of  training  for  eaoh  of  nine  pilots 
is  clearly  both  overly  ambitious  and  not  particularly  useful.  Instead, 
the  matrix,  as  indicated  in  Fig.  27b,  has  been  selectively  filled  in 
those  particular  blocks  where  the  findings  would  shed  the  most  light 
on  the  over-all  model  for  the  pilot.  The  numbers  in  the  blooka  indicate 
rune  by  that  pilot  for  the  condition  noted. 

Besides  using  Judicious  selection  to  decrease  the  magnitude  of  the 
experimental  task,  major  reductions  in  the  amount  of  data  required  can 
be  accomplished  without  decreasing  the  total  scope  by  making  maximum 
use  of  available  data.  Since  an  extensive  amount  of  data  Is  available 
for  Ye  ■  Kq  (Hef.  1 3),  it  was  important  to  establish  oonslstenoy  with 
these  data.  Then  the  sets  of  data  generated  in  the  program  would  be 
closely  enough  related  to  past  work  to  avoid  expansive  and  time  con¬ 
suming  repetition.  Similarly,  since  the  range  of  realistic  input  ampli¬ 
tudes  was  small  and  past  evidence  indicated  amplitude  independence  of 
the  measured  describing  function  over  a  faotor  of  ten  in  0 i,  a  linearity 
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check  was  sale  and  i*  input  amplitudes  were  standardized 

at  a  magnitude  that  was  both  realistic  and  treatable  for  most  condi¬ 
tions—  cj,  m  1  /2"i 

The  subjects  wera  two  highly  experienced  civilian  engineering  test 
pilots,  six  naval  test  pilots,  and  one  light-alreraft-qualified  civil¬ 
ian  pilot  subjected  to  extensive  training.  Thus,  by  dint  of  profes¬ 
sional  experience  and  special  training  the  subjects  were  select  members 
of  an  unusually  homogeneous  group  which  exhibited  high  grade  skills  in 
tasks  similar  to  those  performed  in  the  experiments.  The  subjects  were 
briefed  on  the.  purposes  of  the  experiments,  including  the  anticipated 
ultimate  use  of  the  data  obtained  in  vehicle  handling  qualities  studies. 
Their  level  of  Interest  and  cooperation  was  extremely  high.  Conse¬ 
quently,  by  subject  selection  and  thorough  briefing,  variability  and 
performance  differences  due  to  population  lnhomogeneities  or  laok  of 
motivation  were  minimised. 

To  minimize  variability  due  to  practice  effects,  each  pilot  was 
trained  to  a  stable  level  of  performance  as  measured  by  e^/o£.  In 
general,  this  means  that  at  least  ten  and  often  twenty  trial  runs  of 
two  minutea  duration  with  each  different  controlled  element  were 
carried  out  before  reoorded  runs  were  made  with  a  given  controlled 
element.  Figure  28  is  a  representative  proficiency  curve. 

The  remaining  sections  of  this  chapter  are  devoted  to  the  describing 
function  data  presentations.  The  next  aeatlon  covers  the  establishment 
of  connections  with  existing  data,  and  considers  the  effect  of  forcing 
function  amplitude  on  the  describing  funatlon  for  a  orltio&l  condition. 
The  third  section  treats  the  variability  of  describing  funotion  data  by 
examining  selected  data  sets  which  Illustrate  intra-  and  lntersubject 
effects.  Finally,  the  fourth  section  of  the  chapter  presents  grand- 
average  describing  functions  which  exhibit  the  central  variations 
examined  in  this  study —  those  caused  by  controlled  element  and  forcing 
funotion  changes.  This  section  also  subjects  some  important  qualitative 
conclusions  to  statistical  examination. 
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In  order  to  tis  In  with  the  body  of  date  generated  from  ¥ e  ■  Kg  *  1 
In  Ref.  1 J,  a  tracking  situation  as  similar  as  feasible  to  Elklnd's  had 
to  be  considered.  Ideally,  the  tle-ln  experiments  should  be  conducted 
with  forcing  functions  and  manipulator  similar  to  those  to  be  used  In 
our  other  experiments,  yet  also  similar  enough  to  Elklnd's  to  effeot  a 
reasonable  connection.  Fortunately  Elklnd's  B6  forcing  function  amounts. 
In  our  notation,  to  q  ■  3*0,  1",  so  the  cc*.  ■  2.5,  1"  forcing  function 
was  thought  to  provide  reasonably  close  approximation.  The  lightly 
restrained.  Btick  manipulator  used  in  our  series  differs  substantially 
in  form  from  Elklnd's  freely  moving  pencil-like  pip  tracker,  although 
the  movements  in  both  oases  were  generally  lateral  (with  more  rotation 
involved  in  our  series) .  Yet,  in  our  past  work  ve  had  been  able  to 
show  reasonable  oanneotlons  with  Elklnd's  data  even  using  am  aircraft 
oenter  stlok  (Ref.  34) ,  so  any  differences  due  to  the  manlpulatora  were 
expected  to  be  slight.  Consequently  a  aeries  based  on  the  use  of  an 
u>l  ■  2.3,  1"  forcing  function  was  planned. 

Three  highly  trained  pilots,  Nos.  2,  4,  and  6,  tracked  two  runs 
eaah  for  Ya  ■  Kq  ■  1 ,  2,  and  3,  respectively .  The  differences  between 
Yp  measurements  for  successive  rune  for  eaob  pilot  were  very  slight, 
and  the  two  rune  were  averaged.  These  data  were,  in  turn,  averaged  in 
the  YpY0  form  to  make  them  comparable  across  Y0  values.  These  averages 
are  shown  in  Fig.  29,  with  the  hash  mark*  indicating  the  range  of  the 
data,  i.e.,  the  looatlon  of  the  high  and  low  pilot  averages  used  in  the 
three  pilot  grand  averages.  Elklnd's  oompmrable  data  for  an  ■  3*0,  1" 
are  shown  for  oonparison  in  Fig.  29.  These  data  are  averages  of  four 
four-minute  rune,  two  from  one  sub J sot  and  one  run  from  each  of  two 
other  subjects.  It  is  dear  from  Fig.  29  that  the  STI-FIL  results  are 
remarkably  compatible  with  the  Elklnd  data.  In  fact,  this  extremely 
close  correspondence  between  data  taken  years  apart,  by  different 
experimenters  at  different  looatione,  with  different  subjects,  different 
analysis  apparatus,  and  slightly  different  forcing  functions  and  manlpu¬ 
latora,  etc.,  is  very  satisfying.  On  the  basis  of  this  most  compelling 
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evidence  it  wan  concluded  that  extensive  measurements  with  Yc  a  Kq  were 
unnecessary  and  that  Elkind's  data  could  and  should  reasonably  be  con* 
sidered  a  subset  of  ours  to  the  extent  required  in  model  building  and 
verification. 

Another  important  item  preliminary  to  the  major  experimental  eeries 
relates  to  the  effect  of  forcing  function  amplitude.  Elkind  had  demon* 
stnated  (Hef.  1 3),  by  means  of  comparleons  at  two  low  frequencies  (less 
than  <ife)  of  the  closed-loop  describing  function,  Yp/(1  +Yp),  that  theee 
closed-loop  measurements  are  independent  of  forcing  function  amplitude 
over  the  range  from  0.1  to  1  inch  rms.  This  result  is,  of  course,  totally 
in  keeping  with  what  would  be  expected  In  a  good  closed-loop  system.  Of 
much  greater  interest  is  the  general  effect  of  forcing  function  ampli¬ 
tude  on  Yp.  Tor  reasons  of  experimental  economy  one  might  hope  that  any 
euoh  effeot  would  be  small  over  a  reasonable  range,  although  on  other 
grounds  (e.g.,  see  Table  Z)  it  is  certain  that  some  amplitude-sensitive 
effects  are  present.  Consequently,  in  order  to  examine  the  appropriate¬ 
ness  of  the  restrloted  forcing  function  amplitude  values  in  Tig.  27b, 
a  check  of  linearity  was  mads  for  a  controlled  element,  Y0  ■  ^/(Joa)2, 
whlab  is  difficult  to  aontrol  end  which,  on  an  k  priori  basis,  may  possibly 
result  in  significant  nonlinear  controller  aotlon.  A  typical  comparison 
for  one  pilot  controlling  Y0  -  10/(jou)2  is  prssentsd  in  Tig.  JO,  in  which 
the  avenges  for  three  runs  for  *  1.5,  lA"  end  cu*  ■  1.2,  1  /£"  ere 
presented  for  Pilot  8.  The  comparison  was  mads  with  an  individual  pilot 
rather  than  an  avenge  over  pilots  so  as  to  aacentuats  any  possible  dif¬ 
ferences  rather  than  obscure  them  by  averaging  additional  pllota.  An 
examination  of  these  data  in  Tig.  30  makes  it  dear  that  there  ia  no 
evidence  of  nonlinear  behavior  in  the  aanae  of  a  describing  funotion 
variation  with  cj_.  This  point  will  be  examined  further  subsequently 
when  the  output  spectrum,  *oe(cn),  is  examined  in  fine  detail. 

o.  vajrxaixlxtt  or  moxoan  rumor  data 

The  condition  under  whloh  the  measurements  of  the  pilot's  desorib ' *<« 
funotion  and  associated  characteristics  of  human  control  behavior  were 
made  was  that  of  a  highly  skilled  and  highly  motivated  operator  performing 


94 


«•»»)  V  * 


Figure  50.  Comparison  of  Yp  for  Different  Forcing  Function  Amplitude*; 

Y0  ■  Xc /(d«)  »  “i  ■  1 •? 
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a  carefully  *nd  tightly  defined  task.  The  pilots  were  told  Ilia  l 
although  the  control  and  display  were  only  Idealisations  of  flying 
situations,  there  were  sisd. lari ties  in  that  the  controlled  elements 
could  be  thought  of  as  being  representative  of  certain  aircraft-like 
dynamics  in  the  case  of  Ye  ■  Kc/jco,  of  space-vehicle-like  dynamics  in 
the  case  of  Y0  ■  Kc/(jcu)2,  and  of  an  unstable  vehicle  In  the  cases  of 
Yc  ■  Kq/ (Jus  —  l/T)  and  Y0  ■  Kq/Jc^Juj  -  l/T).  The  instructions  were  to 
minimize  the  error  and  to  attempt  to  do  this  in  the  context  of  their 
flying  experience. 

By  confining  the  tests  to  this  highly  motivated,  highly  skilled, 
carefully  selected  sample  of  the  population  the  intent  was  to  reduce 
the  effects  of  population  inhomogeneities  on  the  data.  Still,  sources 
of  both  intra-  and  inter subject  variability  exist  over  and  above  the 
expected  variations  due  to  changes  in  controlled  element  dynamics  and 
forcing  function  characteristics.  The  nature  of  the  subject-centered 
variations,  for  describing  functions,  will  be  examined  below. 

1 .  Intrapilot  Variability—-  Run -to -Run  and  Differential  Xq  Iffeote 


a.  Run-to-run  variability.  The  first  intrapilot,  variability  of 
interest  is  of  a  run-to-run  nature  —  a  pilot  compared  with  himself 
when  he  tracks  the  same  input  successively. 


The  representative  individual  run  data  and  the  four  controlled 
elements  used  to  examine  the  repeatability  of  successive  runs  by  the 
same  pilot  may  be  characterized  as  follows: 


Fig.  No.  Pilot 


Controlled  Element 


31 

32 

33 

34 


8  k/Jgu:  Easiest  to  control 

6  K/(J<b-2):  Unstable,  but  is  easy 

to  control  if  the  pilot  is  atten¬ 
tive 

8  K/(jo>)2:  Difficult  to  control  and 

requires  considerable  practice  to 
do  well 

8  K/ju>( ju>- 1 .5) :  Very  difficult  to 

control;  momentary  lapses  of  atten¬ 
tion  will  cause  a  loss  of  control 
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On  examining  Figs.  JM  —  3h  we  observe  t.bet.  ail  four  controlled 
elements  indicate  considerable  run-to-run  variation  at  low  frequencies. 

There-  is  much  less  variation  In  the  mid-band,  which  is  in  the  region  or 
crossover,  and  then  a  bit  more  at  the  higher  frequencies  beyond  cross¬ 
over.  TAils  hind  of  variation  1b  particularly  well  displayed  In  Figs.  31 
and  33.  As  discussed  In  Chapter  III,  this  behavior  Is  consistent  with 
the  demands  of  the  closed-loop  system.  The  influence  of  Yp  variations 
on  closed-loop  system  behavior  is  relatively  small  at  low  frequencies, 
but  the  behavior  of  Yp  in  the  region  of  crossover  has  a  large  influence 
on  closed-loop  performance,  which  is,  after  all,  what  the  pilot  observed. 

For  the  unstable  elements,  Kc/(jco— 2)  and  Kc/jo(joJ  —  1 .5) ,  there  is 
evidence  of  constrained  behavior  throughout  the  entire  measurement  range. 

This  Is  particularly  apparent  for  both  amplitude  ratio  and  phase  for 
Kq/( Jcn-2),  Fig.  32.  With  the  very  difficult  task  of  Fig.  success¬ 
ful  tracking  forced  a  restricted  form  of  gain  variation  behavior  on  the 
pilot,  although  the  phase  variations  appear  quite  large.  These  occurred 
beoause  the  task  was  difficult  enough  to  Induce  brief  periods  of  loss  of 
control  due  to  Inappropriate  variations  in  temporal  action —  literally 
time-varying  phase  shifts.  These  are  reflected  by  phase  variability  in 
the  measurements. 

In  summary,  pilots  are  capable  of  extremely  high  repeatability  where 
it  Is  necessary i  e.g.,  at  crossover  or  in  the  control  of  an  unstable  Yc 
where  Yp  is  rigidly  constrained.  The  extent  to  which  this  repeatability 
can  be  maintained  is  a  function  of  frequency  (relative  to  o^)  determined 
by  the  form  of  the  controlled  element  and,  near  the  limits  of  control 
[e.g.,  Kq/Jo3(Jo3-1  .3)] ,  a  function  of  the  task  difficulty. 

b.  Effect  of  Xp  variation.  We  have  indicated  previously  that  the 
form  of  the  controlled  element  exerts  a  large  influence. on  the  amount 
of  variability  exhibited  by  the  pilot.  Wow  we  can  proceed  to  vary  the 
parameters  of  the  form  to  determine  whether  these  are  the  source  of  the 
variation.  This  was  done  for  Pilot  8  for  a  gain  range  of  50  to  1  in  Kc/jcu 
and  a  range  of  25  to  1  in  Ko/(ja>)^,  as  shewn  in  Figs.  35  and  36.  In  general, 
these  data  indicate  that  changes  in  controlled  element  gain  are  offset 
by  the  pilot,  thereby  constituting  further  verification  of  Adjustment 
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Figure  36.  Effect  of  Kq  Variation; 
Yc  -  *L  -  2-5 


Rule  5  (Chapter  II).  There  are  -variations  with  K*.  but  these  are  more 
or  lase  random  and  appear  to  amount  to  the  same  thing  as  the  run-to-run 
-variations  discussed  above.  That  is,  if  the  data  are  considered  only 
to  be  representative  of  successive  runs,  then  the  general  trends  of 
variability  with  frequency  Indicated  in  Figs.  35  and  36  are  very  similar 
to  those  present  in  lign.  31  and  33.  On  this  basis  the  conclusion  noted 
above  that  repeatability  is  a  function  of  Y0  tom  can  be  refined  to  the 
statement  that  the  run-to-run  variability's  dependent  on  Yc/K 5. 

fl.  Iatersuhjeot  Variability 

By  compiling  data  for  the  most  numerous  subject  situations  the  impact 
of  subjects  on  the  variability  of  the  describing  functions  can  be  examined. 
This  Is  done  In  Figs.  37 -to,  where  data  similar  to  those  In  Figs.  31  —  314- 
are  shown  but  with  subjects  replacing  runs  as  the  key  variables.  The  data 
presented  in  these  figures  are  obtained  by  averaging  two  successive  runB 
for  each  pilot,  save  for  Pilot  8  for  whom  three  runs  were  averaged.  AIbo, 
the  data  baBe  is  made  larger  for  the  two  unstable  con  ’-.rolled  elements  by 
using  all  K0'8  (Figs.  38  and  to).  This  procedure  is  unlikely  to  introduce 
much  more  variation  because  the  effect  of  Kc  changes  are  largely  offBet 
by  the  pilot j  but  even  more  dramatic  evidence  to  this  end  are  the  extremely 
low  ranges  apparent  In  the  phase  data  of  Fig.  38»  and  also  in  the  ampli¬ 
tude  ratio  data  when  the  different  K0  values  are  taken  into  account. 

Again  we  notice  the  same  general  trends  as  already  observed  for  the 
run-to-run  and  Kc  lntrasubject  changes.  That  Is: 

a.  For  Kc/Cj©)  and  Kc/(j«)2  fairly  wide  ranges  of  variation 
at  frequencies  much  smaller  than  or  much  greater  than  a*,, 

and  tighter  limits  on  the  range  in  the  crossover  region. 

b.  For  the  unstable  controlled  elements  the  range  of  vari¬ 
ation  In  amplitude  ratio  across  the  entire  frequency 

range  is  generally  small. 

c.  There  is  a  rather  wide  variation  in  the  phase  for  the 
case  involving  the  critically  difficult  controlled 

element  Kc/j©(j©-1 .5).  Again  this  is  attributed  to  the  time 

variations  in  behavior  implicit  jn  the  brief  periods  in  which 

control  was  nearly  lost. 
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An  additional  fact,  whiuh  i»  mant  easily  seen  m  Mg.  o'i,  should 
also  be  noted.  There  le  a  great  deal  of  indifferent  behavior  before 
and  after  crossover,  but  each  pilot  appears  to  be  indifferent  in  bis 
own  way*  Thus  individual  pilots  may  express  their  tracking  style  where 
it  isn’t  critical  for  closed-loop  performance,  but  everyone  behaves 

V  ■ 

nearly  identically  under  the  constrained  conditions. 

0.  (WASD-AVXXMI  mmSM  VUMOVXGM  AID  MATXJTIGAl  ASALSIH 
1 .  Describing  function  Averages 

The  central  aspects  of  manual  control  in  this  study  are  the  nature 
of  pilot  adaptation  and  optimalization  induced  by  controlled  element 
and/or  forcing  function  changes.  The  most  vital  data  to  show  these 
effects  are  averages  of  |YpY0|dfe  for  fixed  forcing  function  condi¬ 
tions.  Averages  of  this  nature,  together  with  ±1 o  bands,  are  given 
in  Pigs.  4l ,  42,  and  4?  for,  Y0  -  Kq/Jo,  Ko/(jo>-2),  and  Kc/(jcu)2, 
respectively.  Parts  a,  b,  and  o  of  these  figures  correspond  to  values 
of  ay,  »  1,5,  2.5,  and  4,0  rad/sec,  respectively.  Figure  44  shows  a 
similar  set  for  Ye  -  Kc/jo{jc-l/T),  a*  •  1 .?,  l/4"j  these  are  segre¬ 
gated  as  Farts  a,  b,  c,  and  d  corresponding  to  values  of  l/T  -  0.0,  0.5, 
1.0,  and  1.5  rad/sec,  respectively. 

The  data  in  Figs.  4l  -44  constitute  major  findings  of  the  study. 

In  Chapter  VII  they  will  be  extensively  analyzed  and  interpreted .  Here 
it  is  appropriate  only  to  summarize  certain  salient  features  of  these 
figures  as  data  per  se  in  the  context  of  the  analytical-verbal  model. 

One  of  the  moat  Interesting  general  aspects  of  the  describing  function 
data  is  the  tendency  for  the  amplitude  ratios  to  approximate  -20  db/daoade 
slopes  throughout  the  measurement  range.  This  tendency  is  adhered  to 
quite  well  for  the  Y0  ■  Kq/Jo  and  Ko/(jo>)2  data  (Figs.  4l  and  44),  whereas 
for  the  unstable  controlled  elements  the  tendency  is  most  prevalent  only 
in  the  immediate  region  of  crossover.  In  fact,  near  the  controllability 
limits  the  slope  becomes  considerably  shallower  than  -20  db/decsde.  All 
of  this  is  in  general  agreement  with  the  pre-experiment  predictions. 
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Another  general  feature  of  great  interest  is  the  universal  presence 
of  lev  frequency  phase  lags.  Xu  edmout  all  oases  these  lagging  phases 
are  not  associated  with  amplitude  ratio  break  points  within  the  measure¬ 
ment  band.  Such  breals  points  oust  accordingly  occur  at  frequencies 

I 

somewhat  lower  than  those  in  the  forcing  function.  The  analytical-verbal 
model  makee  no  provision  whatsoever  for  the  low  frequency  phase  lags, 
and  consequently  should  be  modified  to  account  for  their  presence. 

Summary  open-loop  describing  function  data  of  supreme  importance  are 
crossover  frequencies  and  phase  margins.  Figure  45a  presents  crossover 
frequency  as  a  function  of  forcing  function  bandwidth  for  the  three  con¬ 
trolled  elements  for  which  was  a  variable,  and  Fig.  45b  provides  the 
same  information  for  the  variable  divergence  Yc,  ■  1.5,  1  /V',  Con- 

a 

side ring  the  data  and  curve  reading  tolerances,  the  crossover  frequencies 
for  the  three  variable -mj,  controlled  elements  (Fig.  45a)  each  appear  to 
be  essentially  oonstant.  A  fair  case  can  be  made  for  a  slight  increase 
in  tuc  as  increases,  but  the  maximum  deviation  from  the  mean  value  is 
less  than  5  percent.  The  most  significant  deviation  from  a  oonstant 
crossover  occurs  for  Yc  ■  K0/(J«)2  where  the  »  4  point  has  regressed. 
All  of  these  data  are  in  substantial  accord  with  two  of  the  %  Invari¬ 
ance  properties,  Adjustment  Buies  5b  and  c,  and  hence  offer  dlreot  ■ 
verification  and  validation  for  the  analytical-verbal  model. 

The  average  crossover  frequencies  shown  in  Fig.  45a  for  Kc/jcn  and 
Kc/(j®)2  occur  in  the  relative  order  and  have  values  within  the  probable 
ranges  predicted  in  the  pre-experiment  analyses.  The  experimental  ci^'s 
are,  however,  about  20  percent  lower  than  the  estimates  of  the  center  of 
the  probable  crossover  range  (see,  e.g.,  Fig.  7)* 

Phase  margin  data  derived  from  Figs.  4l  -44  are  summarized  in  Fig.  46. 
In  general,  it  appears  that  phass  margin  increases  directly  with  forcing 
function  bandwidth.  Initially,  for  small  a*'*  tha  trend  Is  linear, 
rounding  off  as  increases.  As  promised  previously,  these  phase  margin 
data  will  result  in  significant  improvements  in  Adjustment  Rule  4b.  The 
phase  margins  for  Y0  -  Ko/jm(jm  -  l/T)  shewn  in  Fig.  46b  are  all  very  low. 
This  is  due  to  the  smfell  tiy.  (1 .5  rad/sec)  used  for  these  cases  and,  for 
the  higher  l/T  values,  the  marginal  controllability.  The  phase  margin  and 


Figure  4l  •  Open-Loop  Describing  Functions  for  Yc  =.  Kg/j a> 


fe2.  Open- Loop  Describing  Functions  For  Yc  =  !£<./(  ja>  —  2) 


Crossover  Frequency  ,  uc  (rad /sec)  Crossover  Frequency  ,  wc  (rad /sec) 


^  »  4.9  ,Kc/(jw«2) 

- Oc  -  4.75,  Kc/jw 
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Forcing  Function  Bandwidth  ,  euj  (rad/sec) 


Crossover  Frequency  Variation  with  Fording  Function  Bandwidth 
for  Yc  -  Kc/j<o,  Kc/( Jco)2,  and  Kq/CJoj  —  8) 


0  0.5  i.O  1.5  2.0 

Divergence  Pole,  I  /  T  (sec-1) 

(b)  Crossover  Frequency  Variation  with  1  /T  for  Yc  -  K0/jo>(jco  -  l/T), 

*  a*  -  1  •%  1  A" 

* 

Figure  45  .  Crossover  Frequency  Variations 
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(a)  Phase  Margin  Variation  with  forcing  Function  Bandwidth 
for  Yq  -  Ko/jto,  Ke/(jaj)  ,  and  ^/(jco-  2) 
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(b)  Phase  Margin  Variation  with  l/T  for  Y0  «  Kc/jo)(ja>  -  l/T) 

a*  -  1.5,  1/V 

Figure  46,  Phase  Margin  Variations 
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crossover  (in  Pig.  45)  shown  for  this  Y0  with  l/T  »  0  differs  slightly 
from  those  for  Yc  ■  Kc/(d»)®  because  different  populations  were  used 
(see  Pig.  27). 

8.  Statistical  Comparisons  of  the  Data 

Although  the  nature  of  this  investigation  was  such  that  detailed 
statistical  comparisons  of  data  were  not  indicated  in  general,  aggregations 
such  as  those  shown  in  Pigs.  4l  and  43  for  Yp  -  K0/jcs  and  Yp  «  Kq/ ( Joo)2 
containing  between  15  and  21  individual  runs  make  some  simple  statistical 
comparisons  feasible.  While  the  numbers  in  the  sample  a-s  such  as  to  be 
on  the  ragged  edge  for  the  applicability  of  omall  Bample  statistics,  the 
following  assessments  can  be  made  for  the  aggregated  data: 

a.  Are  the  phase  and  db  amplitude  data  for  YpYc  normally 
distributed?  If  they  are  normal,  and  past  experience 

indicates  that  they  will  be  (Refs.  15  and  50),  the  Btatistl 

cal  comparisons  of  the  means  can  be  made  more  easily. 

b.  Is  the  se" active  variability  characteristic  of  Yn, 
which  1b  obscured  somewhat  in  the  averaged  YpYc  data, 

a  statistically  significant  observation? 

c.  Are  the  mean  values  of  YpY0  in  fact  different  for 
different  forcing  function  valueB? 

These  <iuestions  are  examined,  in  the  order  stated,  belcw. 

The  small  number  of  data  points  made  it  inadvisable  to  use  a 
chi-squared  goodness  of  fit  test  to  determine  whether  the  distributions 
of  4-YpYo  and  | YpYc | ^  were  Gaussian.  Instead,  the  data  were  converted 
to  standard  scores,  i.e.,  u  ■  (X“X)/o  ®o  that  they  could  be  expressed 
in  termB  of  a  Gaussian  distribution  of  zero  mean  and  unity  standard  devia¬ 
tion.  These  standard  values  were  then  aggregated  by  intervals  which  were 
constant  for  any  one  set  of  runs,  and  cumulative  percentages  computed. 

The  cumulative  values  were  then  plotted  against  the  standard  scores  on 
cumulative  probability  paper.  To  put  the  scatter  about  the  desired 
straight  line  in  perspective,  95  percent  limits  for  random  variations 
about  the  desired  straight  line  were  calculated  for  the  various  sample 
sizes  used  (Ref.  23,  Chap.  6).  This  test  of  the  distribution  of  l?pYeldb 
was  carried  out  for  averaged  |YpYcldb  data  generated  with  Yj  -  K0/jo>  and 


Y0  m  Kc/(jtD)^.  Because  of  the  large  amount  of  effort  involved,  the 
plots  were  ON.lnu1ji.ted  nt  a  limited  number  of  input  frequeneiea  which 
bracketed  the  crossover  frequency  range.  These  were  <%  ■  0,969,  2.54, 
and  15,8  rad/sec.  All  of  the  36  distributions  studied  were  close  to 
Oaussian  as  defined  by  the  95  percent  limits.  A  typical  plot  is  shown 
in  Fig.  47.  This  plot  is  for  »  2.5,  Yc  ■  ^/(jo)®,  -  2.J4  rad/sec, 

for  both  |Yp|db  and  4  Yp. 

The  seoond  question  raised  above  deals  with  the  extent  to  which  the 
different  variabilities  shown  in  the  averaged  data  as  functions  of  u> 
(opecifically,  the  <%'  s  in  the  forcing  function)  were  statistically 
significant  effects.  To  obtain  a  better  appreciation  for  thiB,  succes¬ 
sive  comparisons  were  made  of  a2  and  of  selected  frequencies  at  either 
side  of  the  crossover  frequency  with  a$  values  at  frequencies  bracketing 
the  crossover  frequency  jefa  ■  a|Xp|dp(03n)  or  ^Yp^^]  *  These  variances 
were  conpared  with  each  other  by  the  ratio  of  variances  erf/ of  -  v2, 
also  known  as  the  F-teBt  (Ref.  23).  The  ratio  of  the  larger 
variance,  of,  to  the  smaller,  a§,  was  always  UBed.  Tables  for  this 
test  dependent  on  N1  and  Ng  for  various  p  levels  are  readily  available. 
The  p  level, which  indicates  the  probability  of  exceeding  the  tabulated 
v2  values  by  chance,  was  arbitrarily  set  at  or  below  the  5  percent 
level.  Table  VI  presents  the  results  of  conparing  the  variances  for 
and  4Yp,  for  both  Yc  ■  Kq/Jcu  and  Kc/(J<b)2,  at  input  frequencies 
0.157,  7*57,  and  13*8  rad/sec  with  the  variances  at  each  of  the  six 
frequencies  from  O.969  to  13.6  rad/aee.  Note  that  each  comparison  box 
is  divided  into  six  segments  —  three  for  | Yp |  and  three  for  4*p  for  each 
of  the  three  o&  values  1 .5,  2.5,  and  4.0  rad/sea. 

In  the  aanparisons  we  note  that  for  the  controlled  element  K/(j<n)2 
in  Table  Vl-b  there  is  significantly  more  variability  in  Yp,  for  botb 
amplitude  ratio  and  phase,  at  the  lover  and  higher  frequencies  (<%  -0.969, 
7.57,  and  13.8  rad/seo),  than  there  is  at  the  frequencies  afc»2.34  and 
4.03  rad/sec  whioh  bracket  the  crossover  region.  Thus,  as  we  depart 
from  the  region  of  crossover  in  either  direction,  the  funnel  pattern 
of  increasing  variability  whicji  appears  In  Figs.  43a,  b,  and  c  becomes 
evident  in  the  statistical  analysis  of  Table  Vl-b  as  well. 
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Figure  47.  Typical  Cumulative  Distributions  of  Amplitude  Ratio  and  Ftiase 
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in  the  column  is  significantly  smaller  than  the  variance  for  |Yp|dt 
at  the  frequency  indicated  in  the  row.' 
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The  picture  for  Table  Vl-a  Is  not  as  clear-cut.  There  is  statistical 
svidsncs  rue z*-  ‘/sriiitbility  st  >»twgwi>g  4-w«> **  **  cf 

crossover,  which  for  Fig.  4la,  b,  and  c  is  from  4  to  5  rad/sec.  This 
evidence  is  Strong  for  the  phase  of  Yp  and  somewhat  spotty  for  the 
magnitude.  This  finding,  of  course,  reflects  the  visual  impression  one 
gets  from  examining  Fig.  4la,  b,  and  c. 

The  third  question  relates  to  the  describing  function  difference  as 
a  function  of  forcing  function  bandwidth,  a^.  Figures  48  -  JO  present 
the  YpYc  data  for  Yc  ■  Kc/jts,  Kc/(jc3-»2),  and  Kc/(jui)2  plotted  with  u>j. 
as  a  parameter.  The  Yc  ■  Kq/Juj  and  Kq/CJcb^  data  show  some  minor  dif¬ 
ferences  in  the  amplitude  ratio  outside  the  crossover  region,  and  major 
differences  in  the  phase  throughout  the  frequency  range.  The  Kq/Qcu—  2) 
data  show  only  the  phase  differences,  and  even  these  are  much  lees  than 
those  for  the  other  two  controlled  elements.  A  "t"  test,  using  bilateral 
confidence  limits  since  the  direction  of  differences  was  not  assumed, 
was  performed  to  determine  the  statistical  significance  of  apparent 
differences  in  the  YpYc  means  for  Ye  -  Kc/jto  and  Kq/(Jo>)2  at  oi  ■  0.969, 
2.54,  and  13.0  rad/sec  (Ref.  23,  Chap.  15).  In  eBBenoe,  these  teBt 
results  indicate  that  the  reasonably  large  differences  between  the 
curves  are,  in  fact,  real  effects,  whereas  the  smaller  differences  are  not. 

Finally,  the  data  for  Yc  -  Kq/Ju>(Jo)— l/T)  are  presented  in  one  plot 
with  l/T  as  a  parameter  in  Fig.  51 . 
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Figure  49*  Averaged  Open-Loop  Describing  Functions  for  Y0  ■  Kq/(Joj-  2) 
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The  describing  function,  Yp,  comprises  only  one  part  of  the  quasi- 
linear  system}  the  remnant  is  an  equally  important  component.  Both  the 
magnitude  and  the  form  of  4*33(0})  are  of  interest —  the  magnitude  for 
its  influence  on  c?  and  e^,  and  the  form  because  of  the  insights  it  may 
provide  into  the  detailed  nature  and  sources  of  the  remnant. 

The  remnant,  $03(0) ,  can  be  measured  directly  by  making  determinations 
of  0cc(<c)  at  other  than  forcing  function  frequencies  and  over  such  inter¬ 
vals  that  the  energy  in  <6CC  at  the  input  frequencies  does  not  contaminate 
the  measurement,  ^(o]  can  also  be  found  at  forcing  function  frequencies 
by  using  measurements  of  p  and  4>cc  in  connection  with  Eq  6,  l.e., 

®rm(“5  ■  [l  ~  (58) 


These  methods  are  compatible  and  complementary. 


Considered  as  a  noise  injection  signal,  the  point  of  application  of 
the  remnant  can  be  moved  from  the  pilot's  output  to  other  locations  in 
the  loop  as  long  as  no  nonlinear  elements  are  passed  In  the  process.  In 
other  words,  the  remnant  may  he  considered  to  be  Injected  at  the  pilot's 
output,  input,  or  someplace  In  between  if  such  pilot  nonlinearities  as 
the  indifference  threshold  are  ignored.  Referring  to  Pig.  12  and  denot¬ 
ing  open-loop  remnant  at  e  and  c  injected  into  the  closed-loop  as  *33 
and  *nnc,  respectively,  the  remnant  forms  are  related  by 


®nn0  ♦nnJXpl2 

M  Vo?  |l  +YpYc?  ’ 


(59) 


119 


In  principle ,  remnant  data  could  be  computed  for  all  the  elements 
in  th*  experimental  plan  shown  in  Pig.  H?.  However,  these  computations 
are  both  time  consuming  and  expensive,  so  only  a  relatively  few  runs 
ware  completely  reduced.  ®b  •  selection  was  based  on  limiting  values  of 
the  relative  remnant,  pft,  and  on  desired  trend-establishing  controlled 
element  and  forcing  funotion  combinations.  pa  data  are  presented  in 
the  second  section  and  the  results  of  the  selection  process  are  given 
it  the  third  section  of  this  chapter. 

The  major  contents  of  the  third  section  are  remnant  and  correlation 
coefficient  data.  The  effects  of  controlled  element,  controlled  element 
gain,  and  forcing  function  on  the  remnant  data  are  explored  and  come 
features  indicating  nonstationary  behavior  as  a  likely  remnant  source 
are  revealed.  The  fourth  and  fifth  Bectlone  present  additional  Informa¬ 
tion  tending  to  isolate  the  remnant  sources.  Amplitude  distributions  of 
signals  in  the  loop  are  given  in  the  fourth  section,  with  particular 
emphasis  on  those  of  a  non-Oaussian  nature.  The  fifth  section  presents 
fine-grained  output  power  spectral  density  measurements  intended  to 
reveal  the  presence  or  absence  of  nonlinear  or  sampling  behavior. 

I.  MLAffXVl  HDttAia  DATA,  ^ 

The  relative  remnant,  p^,  is  a  measure  of  the  ratio  of  oherent 
output  to  total  output  of  the  pilot.  The  value  of  pa  can  range  from 
nearly  zero  to  nearly  one.  It  grossly  reflects  nonlinearities,  non- 
Btationary  behavior,  and  effective  noise  injection,  although  it  is  not 
a  measure  of  any  one  of  these  unless  the  others  are  negligible.  Thus 
low  values  of  do  not  necessarily  imply  nonlinearity,  or  nonstation- 
arity,  or  noise  injection — although  they  do  imply  that  at  least  one 
of  these  effects  is  present  in  significant  amounts,  is  defined  by 
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Presuming  forcing  functions  made  up  of  sinusoids,  this  becomes 


s  *?<%> 

n*1 


XI  <p§K)  +  J  f  W“)  ^ 

n»1  Jo 


(Si) 


where  n<p^(<%)  is  the  amplitude  of  that  portion  of  the  output  power  spectra 
linearly  correlated  with  the  forcing  function  at  frequency  <%.  can  be 
computed  readily  from  data  derived  using  either  the  watt-hour  meter 
analyzer  or  the  spectral  and  cross-spectral  analyzer.  For  data  from 
the  watt-hour  meter  analyzer, 


i  K 

T  £ 

n»1 


h 


+  Vc 


cpiK) 


(6a) 


or,  since 


(S3) 


(64) 


When  the  spectral  and  cress-spectral  density  analyzer  is  used, 
appropriate  formula  is 


Pa 


the 


(65) 


p 

Using  various  of  these  formulas  the  relative  remnant  measure,  p£,  was 
computed  for  almost  all  of  the  data.  The  results  are  presented  in 
Figs.  52  -57  and  are  discussed  below. 
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It  will  be  recalled  from  the  detailed  e^erluieatal  uutrix,  Fig,  27b, 
that  Pilot  8  wan  used  In  an  extensive  exploration  of  controlled  element 
gain  effects  for  Kq/Joj  and  KqAjciJ2.  Averaged  data  (three  rune  for  each  Kq) 
for  this  series  are  shown  in  Figs,  52  and  55,  Pilot  7  was  used  to  investi¬ 
gate  at  least  four  controlled  element  gains  for  Yc  ■  Ke/j»  for  all  forc¬ 
ing  function  bandwidths.  Averages  (two  runs)  from  this  series  are  also 
presented  in  Fig.  52,  In  general,  Figs,  52  and  55  Indicate  a  decrease  in 
pjj*  as  controlled  element  gain  is  increased.  In  fact,  Pilot  8  demonstrated 
one  of  the  very  largest  (p|  *  0.81)  and  one  of  the  smallest  (p|  =  O.O65) 
values  of  relative  remnant  in  the  entire  experimental  program  in  this 
Kg/jm  series. 

Figures  54-56  present  individual  pilot  values  as  functions  of 
forcing  function  bandwidth  for  Ye  ■  K q/2(U,  Kq/(Ju}-2),  and  Kc/(jo)2, 
respectively.  All  of  these  values  are  averaged  over  two  rune  except 
those  for  Pilot  8  for  whom  three  runs  were  used.  It  Is  interesting  to 
note  that,  acroee  the  figures,  Pilot  6  generally  exhibits  the  hlgheBt 
p^  values,  whereas  Pilot  4  sturdily  maintains  his  position  on  the  lower 
side.  The  exceptionally  high  values  of  p^  exhibited  by  Pilot  6  for  many 
of  the  conditions  noted  indicate  a  close  correspondence  between  his 
actions  and  those  of  an  equivalent  constant-coefficient  linear  Byatem. 
Finally,  Fig.  57  presents  averaged  p2  data  for  five  pilots  in  control 
of  Y0  -  Ko/jaitJw-l/T),  a*  -  1.5,  lA".  For  FigB.  54-57  the  values  of 
Kq  for  the  various  cases  can  be  found  by  reference  to  Fig.  27. 

0.  naiMR  vown  uioxral  nanism  ass  oomkatiof  cooncnmw 

From  Figs.  52  -57,  high,  low,  and  typioal  conditions  as  defined  by 
the  average  p2  vere  picked  out  for  further  examination.  These  data 
were  culled  to  find  compatible  sets  which  could  exhibit  the  ranges  of 
Yq,  Kq,  and  aij_  of  interest  in  determining  the  effects  of  these  variables 
on  remnant.  The  conditions  surviving  the  seleotion  process  are  shown 
flagged  in  Figs.  52  -57.  Typical  single  runs  were  then  selected  from 
the  data  represented  by  the  flagged  points.  These  runs  are  listed  in 
Table  VII,  where  each  p^  shown  is  for  a  specific  run. 
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figure  54 .  Relative  Remnant  Versus 
Forcing  Function  Bandwidth; 

Yc  -  Kc /Jcd 


A 


Figure  Jfj.  Relative  Remnant  Versus 
Forcing  Function  Bandwidth; 

Y0  -  KC/(JW  “  2) 


Figure  56.  Relative  Remnant  Versus 
Forcing  Function  Bandwidth; 

Y0  -  Kc/(jo))2 


124 


TABLE  VII 


CONDITIONS  SELECTED  FOB  DETAILED  ANALYSIS 


i 

*0 

INPUT  a* 

OPERATOR 

/ 

0.11 

50/(jco)2 

2.5 

8 

0.65 

2/(jo>)2 

2*5 

8 

0.81 

1/J» 

2*5 

8 

.  0.065 

50/j(c 

2.5 

8 

0.55 

6/jto 

1.5 

7 

0.59 

b/jm 

2.5 

7 

0.39 

6/joc 

1.0 

7 

0.3*+ 

2*5/(jo>)2 

1.5*  lA" 

5 

0.34 

2*5/ja>(jo>  -  1 .5) 

1.5s  iA" 

5 

O.36 

5/(j<0“  2) 

1.5 

3 

0.67 

Cj* 

8 

1 

IO 

£.5 

3 

0.72 

cu 

1 

8 

^r\ 

4.0 

3 

0.35 

5/(jw)2 

1.5 

2 

0.50 

5/(jcn)2 

2.5 

2 

0.69 

5/(j«)2 

4.0 

2 

The  remnants  for  these  selected  runs  were  measured  with  the  spectral 
and  cross- spectral  analyzer*  To  obtain  the  best  accuracy,  the  power 
spectrum  was  measured  at  the  pilot's'  output,  ®Q0,  at  high  frequencies, 
and  at  the  controlled  element  output,  ftm..  at  low  f  re  queue  lee.  The 
uncorrelated  part  of  the  power  spectrum  was  found  by  multiplying  the 
total  power  spectrum  by  (1  -p2) .  To  relate  the  remnant  as  measured  at 
the  controlled  element  output  to  the  pilot* f  output,  these  values  were 
divided  by  the  squared  magnitude  of  the  controlled  element  amplitude 
ratio, 

®nn  ■  T"TS  0  “  P^mm  at  lcw  fluencies  (66) 

l^el 

Otherwise  Eq  58  applies  (e.g.,  at  high  frequencies). 
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Ojm  is  the  closed-loop  remnant  referred  to  the  pilot's  output.  If 
the  measurement  frequency  was  set  to  one  of  the  input  frequencies,  a 
relatively  high  value  of  p  was  generally  obtained,  due  to  the  large  con¬ 
tribution  of  the  correlated  portion  (essentially  a  line  spectrum) •  At 
frequencies  other  than  the  input  frequencies,  the  value  of  p  is  very 
nearly  zero,  so  that  the  remnant  is  measured  directly.  Values  of  the 
remnant  computed  at  the  input  frequencies  generally  fit  a  smooth  ourva 
through  those  measured  between  input  frequencies,  indicating  the  gener¬ 
ally  continuous,  power  spectral  density-like  (l.e.,  significant  line 
spectra  absent)  nature  of  the  remnant.  To  indicate  the  continuous  nature 
of  the  remnant  power  spectral  density,  it  will  hereafter  be  denoted  as 
cpnjj,  following  our  previously  established  convention. 

After  examining  the  remnant  data  in  its  closed-loop  form  as 
and  in  the  open-loop  forma,  Cf^  and  (pnne,  which  result  in  equivalent 
closed-loop  effects,  it  was  found  that  the  highest  degree  of  similarity 
among  remants  for  the  Table  VII  conditions  exists  if  the  remnant  le 
viewed  as  an  open-loop  quantity  injected  at  the  operator's  input,  l.e., 
as  dVme.  This  is  related  to  by  Eq  59,  or 


H 


% 


m 


(67) 


The  values  can  be  further  normalized  by  dividing  by  the  mean-squared 
o 

value  of  the  input,  o|,  so  that  the  result  has  the  dimensions  of  one 
over  frequency,  (rad/seo)”^ .  Zero  db  is  taken  as  1.0  (rad/sec)-1 . 

The  open-loop  remnants  referred  to  the  pilot's  input,  for  the  oases 
noted  in  Table  VH,  constitute  the  basic  data  base  for  considering  the 
effeote  of  controlled  element  gain,  forcing  function  bandwidth,  and 
controlled  element  form  on  the  remnant.  These  data  are  presented  in 
Figs.  58-65.  Associated  with  each  remnant  plot  are  the  corresponding 
correlation  ooeffloient,  p,  values.  It  will  be  recalled  that  values  of 
this  parameter  leas  than  one  indioate  time-varying  behavior  over  the 
measurement  run  length. 
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1 .  Iff  cot  of  Kq  YkrUtloa 

Tiro  sets  of  data,  for  Ye  ■  Kq/Jo  In  Fig.  58  and  for  Ko/(joi)2  la 
Fig.  59#  Illustrate  the  general  dependence  of  open-loop  remnant  referred 
to  pilot ' s  input  on  controlled  element  gain.  Remnant  increases  with 
controlled  element  gain,  but  the  variation  is  not  as  extrema  as  that  of 
the  gain.  There  are  also  substantial  differences  In  p  between  the  two 
controlled  elements  considered,  Indicating  that  the  pilot  is  fluctuating 
more  in  the  course  of  a  measurement  run  for  Y0  ■  Kq/ (jco)  than  for 
Y0  -  Kq/Joj. 

S.  If  feat  of  Toraiag  Function  Bandwidth 

The  effect  of  foroing  function  bandwidth  on  remnant  spectral  density 
and  the  associated  linear  correlation  is  examinod  fcr  three  different 
controlled  elements,  Yc  ■  6/jo,  5/(jm)2,  and  5/(jto-2),  in  Figs.  60, 

61,  and  62,  respectively.  The  Yc  ■  6,/Joj  data  (Fig.  60)  indicate  an 
increase  in  remnant  with  increase  in  forcing  function  bandwidth. 

A  correlation  coefficient  decrease  with  increase  is  also  present  in 
this  Instance.  The  close  correspondence  of  these  trends  and  their 
direct  connection,  if  time  variation  is  presumed  to  be  the  dominant 
cause  of  remnant,  is  a  strong  Indication  that  Just  such  time-varying 
remnant  causes  are  Indeed  dominant  here. 

The  remnants  and  p  values  for  Pilot  2  in  control  of  Yc  ■  5/ ( J<o)2 
are  somewhat  different  for  the  three  foroing  function  bandwldths, 
although  no  orderly  variation  appears  to  be  present  (Fig.  6l ) .  The 
spread  present  is  probably  as  indicative  of  run-to-run  variability  as 
it  is  of  any  (pnnp  dependence  on  the  forcing  function  bandwidth. 

Finally,  Fig.  62  shows  both  the  open-loop  remnant  and  the  correlation 
coefficient  for  1*16  unstable-dlvergenoe  controlled  element,  Yc  »  5/(30 -2), 
to  be  invariant  with  forcing  function  bandwidth. 

In  summary,  the  data  examined  above  indicate  that  the  effect  of 
forcing  funotion  bandwidth  on  the  remnant  can  vary  from  minor  to  none. 
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J.  Bffeet  of  Controlled  Zlaouit  Variation  on  Bemnant 

The  Pilot ; $  data  presented  previously  la  Figs •  38  ttfid  59  were  seen 
there  to  be  dependent  on  controlled  element  gain  for  the  two  controlled 
elements  considered.  To  determine  the  effeot  of  controlled  element  par  se, 
these  data  were  averaged  for  each  controlled  element,  and  the  averages 
obtained  are  compared  in  Fig.  63.  When  data  from  other  pilots  are  added 
to  those  of  Pilot  8,  as  In  Pig.  6k,  It  Is  seen  that  the  major  effect  of 
variation  on  the  remnant  Is  as  much  inter subject  as  inter-controlled- 
element . 

An  exceptionally  Interesting  and  crucial  example  of  a  more  substantial 
remnant  variation  due  to  controlled  element  changes  is  shown  in  Fig.  63. 
Here  the  two  limiting  cases,  l/T  -  0  and  1 .5,  for  the  critical  controlled 
element,  Y0  -  2.5/jo(ja)-l/T),  are  shown.  The  remnant  for  l/T  -  0  is 
considerably  less  than  that  for  l/T  -  1 .5,  while  the  correlation  coeffi¬ 
cient  varies  In  an  opposite  fashion.  The  l/T  "1.5  critical  case,  In 
fact,  has  the  lowest  correlation  coefficients  of  tne  entire  experimental 
series.  These  low  p's  Indicate  a  good  deal  of  describing  function  time 
variation  in  the  course  of  the  run.  As  noted  previously,  the  pilots 
were  often  on  the  verge  of  losing  control.  It  will  be  recalled  from 
Fig.  ik  that  the  amplitude  ratio  was  quite  tightly  held,  whereas  the 
phase  exhibited  wide  run- to- run  variations.  These  data,  taken  In  con¬ 
text  with  those  of  Fig.  63,  are  very  strong  evidence  Indeed  that 
describing  function  time  variations  are  a  major  remnant  contributor. 

S.  AMPLZTUDI  DIWWJUTIOHI  03P  IZQKU4  Df  TO  LOOP 

One  can  attempt  further  Insights  into  the  origin  of  the  remnant  by 
examining  the  amplitude  distribution  for  c(t)  and  e(t)  for  the  runs 
listed  In  Table  VII.  Since  l(t)  has  a  nearly  Gaussian  distribution,  a 
non-Gauasian  distribution  at  c  or,  more  unexpectedly,  at  e  would  Indi¬ 
cate  a  strong  source  of  nonlinear  behavior  and  a  possible  basis  for 
remnant  generation.  The  amplitude  distribution  functions  were  examined 
as  follows. 

The  data  on  the  runs  listed  in  Table  VII  were  sampled  once  per  second 
for  a  total  of  24-0  observations  for  each  quantity  considered,  l.e., 
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either  e(t)  or  e(t).  The  mean  and.  variance  vere  then  confuted,  for  each 
data  set.  Since  the  first  eight  deciles  for  the  standardised  variable  x/uv 
in  a  normal  distribution  are  x/ox«  ±0 . 2J ,  ±0.53,  ±0.83,  and  ±1.28,  multi¬ 
plying  these  standard  values  by  the  measured  and  adding  the  mean  value 
defines  the  actual  bounds  for  all  ten  deciles.  Everything  net  included 
within  x/dg  •  ±1 .28  is  in  ths  plus  or  minus  last  decile.  If  the  data  Vere 
Gaussian,  an  equal  number  of  data  points  would  be  found  in  each  decile. 
Beaause  the  actual  data  are  not  necessarily  Gaussian,  the  subdivisions 
are  not  necessarily  deoils  for  the  actual  data.  Such  possible  differences 
between  expected  and  observed  values  provide  the  basis  for  a  chi -squared 
test  (Ref.  23)  for  goodness  of  fit.  Using  this  test  (for  goodness  of  fit 
between  expected  and  observed  data  in  each  decile),  at  the  0.05  level 
(Ref.  23)  six  c(t)  distributions  and  one  e(t)  distribution  were  determined 
to  be  non-Gaussian. 

In  order  to  approximate  a  probability  density  function  from  these 
data,  bar  graphs  were  created  such  that  the  amplitude  of  each  bar  was  the 
relative  frequency  (number  of  measured  points  within  the  decile  divided 
by  the  total  number  of  points)  divided  by  the  width  of  the  decile.  This 
approximate  distribution  is  symmetrical  about  c/fcfe  -  0  or  e/ce  «»  0,  and 
has  eight  amplitude  bars  since  the  amplitudes  of  the  last  two  deolles 
become  zero  beoause  their  widths  are  infinite. 

The  amplitude  distributions  shown  in  Figs.  66b,  f,  and  g  are  all 
Gaussian  as  indicated  by  the  chi- squared  test — a  somewhat  surprising 
result  by  the  look  of  the  distributions.  All  figures  except  Fig.  66d 
are  measurements  of  o(t)j  Fig.  66d  is  the  one  e(t)  distribution  that  was 

e 

not  Gaussian.  The  numerical  coding  on  each  element  of  Fig.  66  represents 
pilot,  controlled  element,  and  ay.,  In  that  order. 

From  visual  observation  of  tracking  records  and  subjects,  we  can 
assart  that  two  kinds  of  oontrol  action  are  likely  to  be  major  causes 
of  non-Geusslan  output  behavior.  One  is  an  occasional  tendenoy  for 
bang-bang-like  limit-to-limlt  control  activity.  The  other  is  a  tendenoy 
to  oontrol  oontrolled  elements  similar  to  X0  ■  Kq/(J<o)2  by  output  pulses 
having  areas  roughly  proportional  to  the  stimulus.  Idealized  bang-bang 
action  will  result  in  bimodal  distributions  wherein  the  relative  frequencies 
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are  large  at  either  extreme.  The  pulsing  behavior  alao  results  In 
distributions  with  a  relatively  large  proportion  of  amplitudes  near  the 
fringe*)  although  Intemedlate  amplitudes  are  also  likely  to  be  presents 
Several  distributions  shown  in  Fig.  66  are  compatible  with  this  pulse¬ 
like  activity,  although  these  distribution  plots  are  at  best  only  roughly 
Indicative  of  possibilities. 

Of  all  the  non-Oausslan  output  amplitude  distributions)  the  strongest 
indication  of  blmodal  operation  is  shown  in  Fig.  66a  for  Pilot  2 
[Ke/(j£u)®)  “  I*?]*  Further  evldenoe  of  nonlinear  operation  Is  seen 

in  the  low  value  for  of  0.35  for  this  condition  (see  Table  VII).  As 
the  Input  cutoff  f regency  is  increased)  p|  becomes  larger,  0.50  for 
■  2.3  and  0.69  for  ■  4.0.  At  the  same  time  the  amplitude  distri¬ 
bution  becomes  more  nearly  Oausalan,  showing  much  less  tendenoy  for 
blmodal  behavior  (Figs.  66b  and  c).  The  results  of  the  chi- squared 
tests  for  cuj.  ■  2.3  and  4.0  indicated  a  Gaussian  distribution  for 
uy,  m  2.3  but  not  for  tiy.  -  4.0)  which  Is  contrary  to  what  would  be 
expected.  The  wider  Input  bandwidth  should  require  a  greater  number  of 
email  corrections  by  the  operator,  which  would  force  him  to  exhibit  a 
more  linear  characteristic.  However,  the  failure  of  the  ohl-squared 
teat  to  verify  the  expected  behavior  In  this  case  is  not  serious,  since 
the  number  of  classes  on  which  it  was  based  was  small.  The  chi-squared 
test  beoames  more  sensitive  as  the  number  of  classes  on  which  it  is 
based  is  increased.  No  large  difference  between  the  distributions  for 
ojj,  -  2.3  and  4.0  is  evident  from  Figs.  66b  and  c. 

The  effect  of  controlled  element  gain  on  the  tracking  behavior  of 
Pilot  6  is  illustrated  in  Figs.  66e  and  f  for  Kq/(Joi)®  and  in  Figs.  66g 
und  h  for  Ko/ja>.  It  le  interesting  to  note  that  a  non-Gausslan  output 
is  produced  by  the  operator  for  the  email  value  of  Kq  when  controlling 
Kq/(J<o)s,  and  for  the  large  value  of  Kq  when  controlling  Kq/Jui.  However, 
there  is  a  distinct  difference  in  the  shape  of  the  distributions  for 
these  cases.  Figure  66h  (Kq/Jo))  indicates  some  degree  of  blmodal  opera¬ 
tion,  whereas  Fig.  66f  [fcc/(J«j)^]  shows  no  indication  at  all.  Ref arenas 
to  the  p|  values  given  in  Table  VIX  shove  a  much  lower  value  for  large 
values  of  Kq  for  both  controlled  elements  than  for  small  values  of  Kq . 
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Thus,  while  the  large  values  of  Kq  generally  result  in  Increased  remnant, 
the  amplitude  distribution  of  the  remnant  does  not  necessarily  become 
non-Qaussian. 

wBrvs  fUHIIl  IMML  ammrnWm 

Still  another  way  to  investigate  remnant  sources  is  to  examine  0oe(o>) 
In  great  detail  There  have  been  repeated  suggestions  that  the  remnant 
oan  be  explained  by  means  of  a  model  which  includes  certain  pronounced 
nonlinearities  (e.gi,  Ref.  36)  or  by  means  t>f  a  model  which  Includes  a 
nearly  periodic  sampler  (Refs.  5,  6,  and  60).  These  suggestions,  if 
followed  through,  place  requirements  on  what  is  to  be  expected  in  the 
spectrum  of  the  operator’s  output. 

To  some  extent  the  expectations  depend  on  the  analytical  description 
of  the  forcing  function,  l.e.,  whether  $n(o)  should  be  considered  as  it 
really  is  —  a  sun  of  deterministic  sinusoids  —  or  as  an  equivalent  to  a 
random  process.  Either  presumption  may  have  merit,  depending  on  the 
problem  at  hand,  and  sometimes  the  results  desired  are  insensitive  to 
the  viewpoint  taken.  Here,  however,  a  potential  conflict  does  exist — 
the  type  of  output  power  speotrum  due  to  operator  nonlinear  or  sampling 
aotions  will  be  quite  different,  depending  on  whether  the  forcing  function 
spectra  are  a  series  of  lines  or  continuous  In  nature.  To  avoid  any 
such  confusion,  the  only  spectral  data  examined  here  will  be  for  single 
runs.  The  appropriate  foraing  funotion  description  will  then  be  as  a 
series  of  line  spectra. 

With  the  line  spectrum  forcing  funotion,  nonlinearities  in  the 
operator  would  be  expooted  to  result  in  output  spectrum  peaks  which  are 
harmonically  related  to  the  forcing  function  frequencies.  Constant-rate 
sampling  on  the  part  of  the  operator  will  also  tend  to  produce  reourrlng 
peaks  and  valleys  in  the  output  spectrum.  If  the  sampler  is  precisely 
periodic  at  a  frequency  u^,  output  spectral  lines  would  be  expected  at 
frequencies  <%  ±  mu%,  a  ■  0,  1,  2,  3,  Slight  variations  in  stapling 

rate  over  a  measurement  run  would  tend  to  slur  the  lines  into  peak*. 

A  search  for  peaks  of  either  or  both  natures  was  undertaken  by  a 
painstaking  spectral  analysis  of  several  output  spectra  from  c>|  up  to 
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40  rad/see.  Figures  67-69  show  a  survey  conducted  for  a  pure  gain 
controlled  element..  Pilots  2,  6,  and  3  were  examined  for  Y0  -  Kq  ■  1, 

5,  and  10#  respectively.  Tease  data  shew  evidence  neither  for  sampling 
nor  for  nonlinear  behavior.* 

To  determine  if  higher  levels  of  forcing  function  power  at  high 
frequencies  would  modify  the  above  results;  the  forcing  function  was 
changed  from  the  augmented  rectangular  form  to  a  pure  rectangular  shape 
wherein  each  sinusoid  in  the  forcing  function  had  the  same  amplitude. 

The  resulting  forcing  function  spectra  are  referred  to  as  R1 4 
(R  •  rectangular;  1 4  rad/see  ■  approximate  bandwidth).  Figure  70 
indicates  that  this  change  in  forcing  function  spectral  shape  fails  to 
evoke  any  indication  of  sampling  or  nonlinear  behavior. 

Because  the  pilot  sometimes  appears  to  exert  control  by  meanB  of  a 
pulsing  motion  when  Ye  -  ^/(joj)2  (e.g.;  Pilot  2  in  Fig.  66),  the  output 
spectral  density  for  some  of  these  cases  was  examined  for  conditions 
which  might  make  the  presence  of  output  sampling  more  likely.  Figures  71 
72,  and  73  are  for  Pilots  2,  7,  and  8  controlling  ay,  -  1.5,  i/2"  for 
*c  ■  Kq/(Ju>)2  with  Kc  ■  5;  5;  and  10;  respectively.  Again  there  is  no 
evidence  for  either  sampling  or  strong  nonlinear  behavior. 

In  summary;  the  evidence  presented  in  this  section  strongly  indicates 
that  neither  periodlo  sampling  nor  nonlinear  transfer  iharaoteristics  are 
dominant  remnant  sources.  Although  Section  D  presented  some  amplitude 
distribution  data  tending  to  support  a  pulsing  type  of  control  action 
for  oertaln  K3/(jo))2  oases,  this  does  not  appear  to  result  in  substantial 
remnant  relative  to  other  sources.  So  much  for  what  the  remnant  1b  not. 
However,  enough  data  and  Interpretation  have  already  been  presented  to 
Indicate  that  the  major  source  of  remnant  is  nonstationary  behavior 
during  measurement  runs;  and  to  tend  to  isolate  the  major  kind  of  non- 
steady  activity  as  a  time-varying  phase  shift. 


*Th«  scale  is  such  that  1 .0  in2/rad/,sec  on  the  graph  is  zero  db.  The 
tape  noise  level  is  indicated  by  a  broken  line. 
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a.  anacownic* 

The  describing  function  data  of  Figs.  4i  -44  contain  the  quaei- 
linear  transfer  characteristics  of  the  human  pilot  for  the  experimental 
situations  examined.  The  dependence  of  the  operator's  transfer  charac¬ 
teristics  on  the  forcing  function  and  controlled  element  task  variables 
is  explicitly  revealed  by  comparisons  between  these  data.  For  the 
experimental  conditions  considered  vs  can  qualitatively  describe  the 
trends,  features,  differences,  similarities,  etc.,  exhibited  by  the  date 
for  and  aoroas  the  several  conditions.  Such  a  qualitative  description  pro¬ 
vides  on  a  verbal  level,  with  boos  generalisation,  the  information  explic¬ 
itly  exhibited  numerically  by  the  data  themselves,  nils  procedure  has 
already  been  used  in  the  chapters  presenting  the  describing  function  data 
and  will  be  invoked  again  later. 

Qualitative  description,  vhils  necessary,  does  not  provide  a  sufficient 
basis  for  the  abstractions  needed  to  achieve  simplification  and  generality. 
This  oan  best  be  acoonplished  using  quantitative  models  which  In  some 
sense  or  other  approximate  the  actual  data.  When  the  experimental  results 
are  approximated  by  analytical  models,  ths  numsrioal  valuss  of  the  model 
parameters  serve  as  quantitative  measures  of  the  experimental  data.  Then 
we  aan  develop  some  "laws"  of  operation  for  ths  approximate  analytical 
models  in  whioh  the  variations  in  these  parameter  values  exhibit  trends 
similar  to  those  observed  experimentally.  In  this  way  verbal  description 
is  rsduoed  to  simple  quantitatively  stressed  characterisations.  To  the 
extent  that  the  analytical  model  approximates  ths  actual  data  and  ia 
reasonably  simple  in  form,  we  may  obtain  much  greater  Insight  into  the 
"laws"  whioh  actually  govern  the  operator's  behavior.  To  aohlave  this 
end  vt  shall  evolve  in  this  chapter  eoma  simple  analytical  models  for  the 
data  of  Figs.  4l  -44. 
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The  purposes  which  the  analytical  models  are  Intended  to  serve 
elude: 

1 .  Characterization  of  past  observations  for  fixed 
experimental  conditions  with  simple  analytical  . 
forms,  thereby  achieving  economies  in  description 
of  the  experimental  data 

2.  Quantitative  description  of  experimentally  observed 
trends  in  terms  of  analytical  relationships  between 
model  parameters  and  task  variables 

3.  Establishment  of  the  basis  for  a  rationale  which 
can  be  used  to  better  "understand"  bhe  ways  in 
which  the  pilot  behaves  as  a  control  system 
component 

4.  Provision  of  a  basis  for  quantitative  extrapolation, 
thereby  making  possible  the  prediction  of  pilot 
behavior  in  novel  situations 

With  these  purposes  in  mind,  turn  next  to  the  possible  types  of  models 
which  might  be  used  to  characterize  the  data.  Clearly,  a  highly  desir¬ 
able  feature  would  be  some  form  in  which  a  minimum  number  of  parameters 
are  capable  of  representing  the  phenomena  selected  for  emphasis  in  that 
particular  model.  Equally  important  is  the  desire  that  the  particular 
model  form  be  one  which  is  analytically  simple  and  tractable.  These 
desires  can  conceivably  conflict  with  the  requirement  that  the  analyti¬ 
cal  model,  be  reasonably  descriptive  of  the  data.  Fortunately,  both 
desires  and  requirements  can  be  met  by  using  describing  function  forms 
made  up  largely  of  ratios  of  rational  polynomials  with  the  addition  of 
a  transport  lag  (pure  time  delay)  term.  The  uee  of  such  analytical 
models  for  human  operator  describing  function  description  has  been 
extensively  discussed  and  applied  elsewhere  (e.g..  Refs.  1J,  34,  and  53) 
and,  in  fact,  the  existing  analytical-verbal  model  is  of  precisely  this 
nature.  Our  major  task,  therefore,  is  to  examine  the  data  presented 
here  using  such  forms  rather  than  to  Justify  the  particular  variety  of 
form  selected. 

In  the  following  sections  analytical  models  shall  be  derived  by 
curve-fitting  the  experimental  describing  function  data.  In  the  process 
three  baBic  levels  of  approximation  shall  be  used.  At  the  first  level 
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are  extremely  simple  models  which  exhibit  characteristics  similar  to  the 


actual  data  in  the  crossover  region.  •The  forma  yaert  •mount  to  the  cross¬ 
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When  both  o^,  and  Ta  ai-  allowed  to  '/ary,  this  equation  has  enough  flexi¬ 
bility  to  characterize  the  baslo  crossover  region  trends  for  the  Kc/jm 
and  Ko/fjcu)2  data.  With  somewhat  less  accuracy  it  is  also  suitable  for 
the  Kq/(Jcc— 2)  cases.  With  this  first  model  the  variations  of  paramount 
importance  are  those  of  Te  and  with  controlled  element  and  foraing 
function  bandwidth.  In  the  course  of  the  discussion,  Adjustment  Rules  4b 
(system  phase  margin)  and  5b  (o^-cuu,  dependence)  will  be  modified. 

Modified  versions  of  the  crossover  model  are  more  appropriate  tc 
systems  involving  extensive  pilot  lead  and/or  those  having  controlled, 
element  dynamics  with  break  points  in  the  general  crossover  region.  Hhese 
include  such  open-loop  describing  functions  as 
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for  Yc  ■  Kc/.lu^Jco  -  1/t)j  and 
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for  Yc  -  Ko/fjm—  1  /T) .  These  direct  extensions  of  crossover  models  are 
the  simplest  forms  possible  whlah  are  compatible  with  even  a  gross  char¬ 
acterization  of  the  data  over  the  measurement  bandwidth. 
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Both  the  crossover  nodal  and  its  direct  extensions  fail  to  acoount 
for  the  low  frequency  phase  lags.  These  are  handled  by  refining  the 
extended  crossover  model  to  inolude  a  catchall  increment  to  the  lav 
frequency  phase. 

With  the  extended  crossover  model,  all  of  the  major  trends  supported 
by  several  data  points  over  a  portion  of  the  measurement  bandwidth  are 
described  reasonably  veil.  There  are,  however,  certain  isolated  data 
points  which  require  further  explanation  or  which  can  serve  as  a  basis 
for  further  elaboration.  These  isolated  points  occur  at  the  extremes  of 
the  measurement  bandwidth,  where  data  variability  tends  to  be  largest. 
However,  for  the  Y0  ■  !£<>/( joi  -  2)  case  the  constraining  nature  of  the 
controlled  element  dynamics  acts  to  keep  the  variability  low,  so  these 
data  offer  admirable  possibilities  for  further  elaboration  and  refine¬ 
ment.  This  is  accomplished  using  "precision  fits"  for  the  data  vhloh 
support  the  so-called  precision  model.  This  most  refined  model  is  the 
only  one  which  has  a  basic  form  which  is  not  necessarily  restricted  to 
use  over  only  a  finite  bandwidth,  l.e.,  it  satisfies  physical  limita¬ 
tions  as  o>  approaches  zero  or  infinity  while  characterizing  the  data 
exceptionally  well  over  the  measurement-  band. 

The  three  baaic  levels  of  model  precision  introduced  here  will  be 
dlsoussed  in  the  three  Immediately  following  sections.  An  «  fin*l  pre¬ 
view  to  help  keep  things  straight  for  the  reader,  a  t  itual 

model  types  are  used.  These  are  summarized  below  »(.f:  .  -  of  the 

parameters  vhlch  are  allowed  to  vary  In  order  to  f  '  ■<  ■  a  given  Y 

Crossover  Models 

.  u0rJayr« 

VpY0  i  - rj —  ,  u*j,  variable  (72) 


Vo  1  —  »  *e  variable  (75) 


•  5W*rAU..r.«  V 


“  t.  tdAiilc 


:»i|i 


Extended  Crossover  Models 

First  Approximation! 


Yp  -  Kp  (jd - |  Kp,  TL,  Tj,  Te  variable  (7*0 

\Txd»  +  1/ 

Ssoond  Approximation:  - 

YP  ;  Kp(^3^7r)e~J  (uyre  +  J  Kp,  TL,  Tj,  «reJ  a  variable  (75) 


Precision  Model 


The  final  section  of  the  chapter  discusses  the  mean- squared  error 
performance  measure,  crossover  frequency  regression,  and  related  aspects 
of  the  data.  In  this  section  the  nature  of  the  gain,  i.e.,  crossover 
frequency,  and  phase  margin  adjustments  and  their  relationships  with 
average  performance  are  explored  in  terms  of  the  crossover  model. 


>.  m  CKMOVnt  MOBIL 


Visual  examination  of  the  data  shown  in  Figs.  48-50  indicates  that, 
to  a  first  order  at  least,  the  crossover  model  discussed  in  connection 
with  Bq  10  is  likely  to  have  a  fair  range  of  applicability.  As  noted 
when  it  was  introduced,  the  model  would  appear  to  be  a  better  deeorlpti  u, 
of  amplitude  ratio  than  of  phase  characteristics.  Also,  the  ampT 
ratio  for  Yo  ■  Ko/(jo>-2)  tends  to  be  somewhat  flatter  than  -SO  -  /('■■■  -  d 
Still,  the  elementary  crossover  model  is  adsquate  to  describe  --.'i 

in  the  crossover  region  using  only  0^  and  %  as  parameters. 
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Although  the  basic  YpYQ  approximate  form  is  the  same  for  all  these 
data,  the  numerical  values  of  crossover  frequency  and  re  are  not.  In 
order  to  separate  the  effects  of  the  tvo  task  variables,  Ye  and  ajj.',' a 
perturbation  treatment  is  appropriate .  Tor 


Vo  " 


o>ce 


rd«ffe 


(77) 


the  crossover  frequency  is,  of  course,  and  the  phase  margin  is 


%  "  5  “  Vte  (70) 

If  ts  and  cue  are  now  divided  into  tvo  components,  one  depending  on  Y0 
alone  and  the  other  on  ci^  and  Y0,  then 

<^(Yo>  «i)  -  «o0Cra)  +  *%(“*>  *c)  (79) 

*e(¥o>  “h)  -  W  -  AT(co1,  Y0)  (80) 


The  "o"  eubsoript  in  Bqs  79  and  80  indicates  values  taken  vhen  a*  ■  0. 
Using  these  equations  vith  the  phase  margin  expression,  and  assuming 
that  the  product  AoqAt  viU  be  negligible,  the  phase  margin  becomes 


•  f  -  *o“b0  +  -  Vkh,  (81 ) 

An  important  consideration  in  carrying  this  development  further  is  the 
value  of  the  phase  margin  vhen  0^-0.  The  phase  margin  data  presentad 
in  Fig.  46a  indicate  a  trend  toward  zero  phase  margin  as  ^  decrease*. 
Although  this  trend  to  the  origin  can  only  be  demonstrated  dlreotly  for 
<0).  0,  it  is  a  common  observation  that  signals  circulate  throughout  the 

manual  oontrol  loop  without  any  foroing  function  as  long  as  the  operator 
la  in  active  oontrol.  In  the  absenoe  of  other  inputs  or  disturbances, 
the  presenoe  of  these  signals  ivpllss  an  on- the -average  condition  of 
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zero  phase  margin.  So,  both  the  dlreot  and  Indirect  evidence  indicates 
that  phase  margin  Is  zero  when  ay,  is  zero.  Because  cifc0  and  t0  are  eq.ual, 
by  definition,  to  on  and  t«,  respectively,  for  «y  ■  0,  then  the  phase, 
margin  for  this  -case  becomes  ...  r  . 

-  f  ~  '*cfa0o  ■  0  (82) 

or 

T°  "  sfe-  .  .  <85>  ■ 

Equation  83  indicates  that  t0  and  u\,Q  are  not  independent  entitles,  and 
that  only  one  need  be  specified.  Since  these  quantities  are  the  only 
ones  in  Eqs  79  and  80  which  depend  on  the  controlled  element  alone,  either 
cjQq  or  t0  must  be  used  to  subsume  controlled  element  differences  as  repre¬ 
sented  in  the  crossover  model  if  and.  At  variations  with  Yc  are  small. 

When  the  above  Considerations  are  taken  into  account,  tie  phase  ' 
margin  expression  simplifies  to 


-  Mc0At  -  ToAlfc 

(64) 

At 

.  ^  +  TA 

(85) 

w 

The  nature  of  the  variations  of  Ai*>  and  At  with  uy  and  controlled  element, 
and  the  variation  of  ay0  or  t0  with  controlled  element,  can  be  explored 
using  crossover  frequency  and  phase  margin  data  extracted  from  Figs.  48-30. 
These  are  presented  in  Table  VIII.  Included  in  the  table  are  data  for 
Y0  •  Kq  tests  which  are  not  documented  in  detail  in  this  report  although 
they  are  indicated  in  the  experimental  plan  of  Fig.  27b.  For  these  tests 
only  three  subjects  with  a  total  of  six  runs  per  forcing  function  condi¬ 
tion  ware  available. 

The  crossover  frequency  is  plotted  as  a  function  of  uy  in  Fig.  74. 

Them  is  a  slight  lnorease  in  uy  with  ay  which,  as  already  noted  in 
Chapter  V,  deviates  from  mean  values  by  only  3  percent  or  so  (with  the 
exception  of  Y0  -  Ky)  For  the  purposes  of  this  model  we  ohall  consider 
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TABUS  VIII 


PHASE  MARGINS  AND  CROSSOVER  FREQUENCIES 


*  1 .3  rad/seo 

mi  ■  2.3  rad/sec 

(Oj.  ■  4.0  rad/sec 

*c 

<PM 

(rad) 

(rad/sec) 

% 

(rad) 

<“e 

(rad/seo) 

(»d) 

(rad/seo) 

*c 

B 

5.1 

0.73 

0.75 

6.7 

Ko/Jo 

E3 

■4.6 

0-73 

0.94 

5.0 

Kc/(J0)-  2) 

0.35 

it. 6 

0.66 

0.70 

5.2 

Kc/(jm)2 

0.26 

3.2 

0.51 

KB 

0.73* 

1.0* 

^Regressive 


the  aotual  variation  indicated  by  the  data.  This  vlll  change  the  a^-u^ 
adjustment  rule  slightly ,  although  to  a  somevhat  lover  level  of  approxi- 
nation  the  adjustment  rule  is  adequate  as  it  new  stands. 

One  of  the  primary  purposes  of  Fig.  74  is  to  obtain  an  estimate  of 
ttb0.  For  Ko/Jcd,  Ko/(j<u  -  2),  and  Ko/(j<o)2  this  is  done  by  fairing 
straight  lines  through  the  data  and  extrapolating  to  a*  ■  0.  For  Y0  ■  K<j 
a  straight-line  extrapolation  vould  result  in  an  ouo0  lees  than  that  for 
either  Y0  -  Kq/Jo  or  Kq/( Jo>  —  2).  Because  all  of  the  existing  data 
points  for  Yq  ■  Kq  are  larger  than  those  for  the  other  controlled  elements 
ouoh  an  extrapolation  does  not  appear  warranted,  so  the  Y0»  K<j  points 
are  oonneoted  with  a  curve  vhioh  for  lew  a*  is  roughly  parallel  to  the 
ourves  for  the  other  controlled  elements.  The  <ifc0  and  associated  t0 
values  are  given  in  Table  R.  This  table  also  oontains  a  recapitula¬ 
tion  of  the  averaged  crossover  frequency,  S^,  from  Fig.  4ja .  When  o^0 
is  subtracted  from  the  <i^j  values  given  in  Table  VIII,  the  resulting 
data  points  coalesce  fairly  veil  vhen  plotted  versus  .  As  can  be 
appreciated  from  an  examination  of  Fig.  74  >  the  major  deviation  from  thiB 
is  Y0  ■  Ko*  The  data  appear  to  be  substantially  independent  of 
controlled  element,  and  a  linear  fit  is  adequate  (except  for  the  larger 
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TABLE  IX 

SUMMARY  OT  T0  AND  a*0 


T0 

(see) 

% 

(rad/eec) 

(rad/sec) 

Ko 

0.35 

4.8 

3-8 

Kc/Jm 

O.36 

4.4 

4.73 

Kc/(J«-  1/T) 

0.3  € 

4.4 

4.9 

Kc/(J<i>)2 

O.51 

3*0 

3.25 

o>L  data  for  Y0  ■  Kq) •  This  ia  given  approximately  by 

Aa^(<J^)  ■  0.l8<uj,  (86) 

Suffloient  information  la  now  available  to  compute  the  incremental 
effective  delay,  At,  using  Eq  8?.  The  reault  la  preaented  in  Pig.  73  , 
which  shove  a  remarkable  lack  of  dependence  of  At  on  controlled  element. 
In  the  low  foroing-funotlon-bandwldth  linear  range, 

At  (ay,)  A  o.oStui  (8?) 

This  reduction  in  effective  time  delay  aa  forcing  function  bandwidth  la 
inoreaaed  la  the  principal  factor  Involved  in  the  phase  margin  increase. 
A  very  almilar  result  la  obtained  if  the  oroaaover  frequency  la  preaumed 
constant,  with  a  value  given  by  3^,  although  the  At  data  points  obtained 
for  <!>(,  m  4  rad/sec  exhibit  a  wider  spread  than  that  shown  in  Fig.  75* 

In  summary,  the  basla  trends  reveals!  by  the  simple  crossover  model 
indole  an  increase  in  the  phase  margin  aa  forcing  function  bandwidth 
la  Inoreaaed,  which  is  aooomplished  by  a  reduction  of  the  high  frequency 
lag  characteristics  subsumed  in  tb.  A  small  part  of  thia  reduction  In 
effective  time  delay  ia  used  to  lnorease  oroasover  frequency,  although 
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this  eas  be  considered  a.  second-order  effect?  Beth  chenge*  *»**  i"  *■ 
direction  which  tends  to  compensate  for  the  inorease  in  mean- squared 
error  due  to  increase#  in  fora lag  function  bandwidth  —  a  topic  covered 
in  detail  in  Section  7.  Finally;  in  the  simple  erfissover  model  the  effect 
of  controlled  element  is  totally  given  by  t0  or,  alternatively,  by  U3qq. 

0.  ID  BOINDID  OtOCSOVll  MODZtL 

The  sing? lifted  crossover  model  discussed  above  is  deficient  in 
several  respects.  First,  for  the  controlled  elements  with  nonzero 
poles  the  open-loop  describing  function  in  the  region  of  crossover 
tends  to  be  somewhat  less  than  -20  db/decade.  For  these  cases  the 
data  are  better  fitted  with  open-loop  describing  functions  which  con¬ 
tain  the  controlled  element  dynamics  explicitly.  Seoond,  the  seoond- 
order  controlled  elements  require  lead  equalization.  This  la  implicit 
In  the  crossover  model;  but  for  more  refined  models  should  be  made 
explioit.  When  these  two  major  deficiencies  in  the  simplified  cross¬ 
over  modal  are  oorraoted;  the  open-loop  describing  function  models  are 
adequate  to  describe  most  of  the  closed- loop  system  characteristics. 

The  extrema  phase  lags  at  low  frequencies  usually  do  not  substan¬ 
tially  affect  the  aloeed-loop  uhareoterlstics  because  magnitude  |YpY0| 
it  much  greater  than  unity  at  the  frequencies  where  the  lags  are  present. 
However,  the  ubiquitous  nature- of  the  low  frequency  phase  lags,  and 
their  occasional  importance  on  closed-loop  dynamioe,  demands  attention. 
This  it  accomplished  here  by  adding  a  oatohall  increment  to  the  low 
frequency  phase  which  takes  into  aooount  dynamlos  having  amplitude  ratio 
break  points  below,  or  in  the  lowest  frequency  portion  of,  the  measure¬ 
ment  bandwidth.  The  basis  for  this  approximation  derives  from  consider¬ 
ing  that  the  low  frequency  phase  lag  is  due  to  an  equal  number  of  lags 
and  leads  having  break  points  generally  below  the  lower  measurement 
frequencies.  Thus,  for  M  leads  and  lags  occurring  at  1  /T]_aai^  and 
l/Tljlg1  the  phase  will  be 

APlow  "  t»n“1  Tleadtw  "  i  tan”1  (88) 
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When  all  Of  wifi  break  points  occur  below  tlia  uimmm urwnmai  bandwidth,  the 
phase  anglt  within  the  measurement  bandwidth  will  be  approximately 


**"  4  I,  {t'ZS&iX'  I  (*'^1 

“  I  (®3.g  Tlead)i 


where 


a 

ao 


“  ‘  ill  (jj1«  '  ^lead)j_ 


(89) 


(90) 


The  effective  time  o one tent  i/a  describee  the  effeot  within  the  measure¬ 
ment  band  of  leads  and  lags  below  measurement  frequencies.  In  this  sense 
It  is  analogous  to  Te,  which  lumps  high  frequency  phenomena  into  a  simple 
low  frequency  approximation  suitable  within  the  measurement  bandwidth. 

In  the  desorlbing  funotion  the  low  frequency  tern  is  represented  ns 

In  the  following  the  various  open- loop  desorlbing  functions  are 
curve-fitted  in  two  ways  which  differ  only  in  whether  or  not  they  include 
the  low  frequency  phase  correction.  To  restrict  somewhat  the  total 
number  of  parameters  varied  for  a  given  controlled  element  form,  the 
crossover  frequencies  are  held  to  be  invariant  with  o^,  for  the  Y0  ■  K0/jo>, 
Kc/( Joi  -  2),  and  Kq/(J(o)&  (with  the  exception  of  the  regressive  ay.  ■  4.0) 
cases.  The  second-order  diverging  controlled  element,  -  l/T), 

situations  amount  to  a  speoial  set  in  which  is  not  constant.  Conse¬ 
quently,  the  first  thres  controlled  elements  and  the  second-order 
diverging  oase  are  separately  treated  below  as  two  grand  families. 

1 .  Extended  Crossover  Models  for  Y0  ■  Kq/Joj,  Xq/(Jc»  -  8),  and  Xo/(jw)a 


Examination  of  the  averaged  YpY0  data  of  Figs.  48-50  reveals  the 
following  general  trends t 
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a*  Crossover  fru<tuauay  is  nearly  independent  of  «j.  (#»** 

Fig.  45a) 

b.  Fhas*  aarginlsnearlypropoFtion^^  46a) . 

o.  The  inaremsntal  lew  frequency  phase  lag  generally 

increases  with  ay  ' 

Two  exception*  to  these  generalisations  need  to  he  noted.  One  is  thq 
ay  •  4.0  regressive  case  for  Y0  ■  Kq/(Ju>)2  and  the  other  is  the  ordering 
of  the  low  frequency  phase  lag  for  to  ■  Kq/Jcj,  where  the  ay  -  4.0  condi¬ 
tion  has  less  lag  than  the  other  two  uy  conditions.  With  these  noted 
exceptions,  the  trends  listed  are  exhibited  by  the  analytical  models . 

a.  Ya  -  The  simplest  models  are  those  for  Kq/Joj.  The 

appropriate  first  and  second  approximation  forme  for  the  analytical 
model*  are: 


Beoond  approximation:  YpY0  ■ 


.  .  vjw* 

Vo  -  Jo 

(91) 

■  3  e-J(«We  +  «*/“) 

Vo  *  °  j. 

(9a) 

The  first  approximation  model  ie  Identical  In  fom  to  the  oroaeover 
model,  except  that  %  is  presumed  to  be  the  oonatant  value  of  Sq.  The 
general  adequacy  of  these  models  is  graphically  illustrated  by  Figs .  76a 
through  760*  The  constants  are  summarised  in  Table  X.  Note  here  that 
1Q)  whioh  is  given  by  x/2B^,  oorrespor.de  to  the  value  of  t*  that  would 
be  present  for  ay.  -  0. 

TABLE  X 

SUMMARY  OF  DESCRIBING  FUNCTION  CONSTANTS  FOR  Y0  -  Kc/ja) 

2\j  ■  4.75  rad/eec  j  r0  ■  0.33  sec 


a 

AT* 

(To  "  U) 

rad/aeo 

eec 

0.14 

0.11 

0.20 

0.088 

0.15 

0.21 
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*•  Ya -Ko/(i»- a).  The  flrat  arid  second  approximation  analytical 
model  forma  for  Y0  “  Ke/ ( J»  “  2)  are  r 


First  approximation t ~ 

(— Jtb/S  + 1)  >. 

(9J) 

Second  approximation: 

YpY0 

(“"e  +  ®/<») 

(-J«/2  +  1 ) 

(94) 

where  K  ie  the  total  open-loop  D.C.  gain.  Figures  TTa  through  77c  illus¬ 
trate  these  fits,  which,  together  with  the  analytical  model  data,  are 
summarized  in  Table  XI. 


TABLE  XI 

SUMMARY  OF  DESCRIBING  FUNCTION  CONSTANTS  FOR  Yc  «  Kc/(ja>  -  2) 

■  4.9  rad/sec 


.. 

% 

Te 

K 

a 

(T0  -  Te) 

deg 

rad 

sec 

rad/ sec 

sec 

1*5 

24.0 

0.42 

0.16 

2.65 

0.2 

0.086 

2.5 

35-0 

0.58 

0.13 

2.65 

0.19 

0.12 

4.0 

39.0 

0.68 

0.10 

2.65 

0.18 

0.14 

For  this  case,  ?0  will  be  given  by 


T0 


J-tarT^ 
<4-  £ 

0.24  BBC 


(95) 


As  can  be  appreciated  by  examining  Fig.  77 1  the  extended  model  is  quite 
adequate  to  describe  all  of  the  data  except  for  the  highest  and  lowest 
frequency  amplitude  ratio  points  for  all  three  u^_  values.  These  dis¬ 
crepancies  can  be  corrected  by  adding  a  high  frequency  lead  and  some  lags 
and  a  low  frequency  lag-lead,  while  making  appropriate  modifications  to 
Te  and  a.  For  instance,  a  is  reduced  to  zero  by  adding  the  low  frequency 
lag-lead  pair  (Tj^Joj  +  1)/(TkJu>  +  1 )  to  the  describing  function  form. 
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Consideration*  such  as  these  are  used  in  the  next  section  to  develop  a 
precision  model. 

a.  Te  >  The  two  forms  for  Yn  -  Ko/(.1cu)2  are: 

First  approximation: 

y 

Second  approximations 

Just  as  with  Y0  -  Xo/j<jo,  the  first  approximation  fits  are  similar  to  the 
simple  crossover  model,  but  with  crossover  frequency  held  constant  at  t^.. 

As  a  prelude  to  the  final  fitting  procedure  for  the  second  approxi¬ 
mation  form,  a  large  number  of  variations  were  used  at  the  low  frequency 
end.  A  comparison  of  these  for  the  three  taj,  conditions  indicated  that 
an  excellent  case  could  be  made  for  a  nearly  constant  Tl  having  a  value 
of  5  sec  or  greater.  In  Figs.  78a  through  78c  the  minimum  value  is  used. 
Table  XU  summarises  all  the  data  for  Yu  ■  Kq/( Joo)s  models. 


v  v 


TABLE  XXI 

SUMMARY  OF  DESCRIBING  FUNCTION  CONSTANTS  FOR  Y0  -  Kq/CJco)2 
■  j.aj  rad/tec  j  Tc  «  0.483  sec 


*L 

a 

ATft 

(To  "  fa) 

deg 

rad 

rad/sec 

sec 

1.5 

5.25 

16.0 

0,279 

0.585 

5 

0.2 

0.098 

a.? 

5.25 

25.3 

0.445 

0.335 

3 

0.333 

0.148 

4.0 

1.8 

42.0 

0.753 

0.260 

5 

0.533 

0.223 

4.  Nummary  Data.  The  key  variations  of  sffectlve  time  delay,  t e, 
with  forcing  funotion  bandwidth,  a*,  are  shown  in  Fig.  79.  The  trends 
shown  are  generally  compatible  with  those  that  would  be  exhibited  if  the 
simple  crossover  model  data  were  presented  in  a  similar  way.  This  oan 


157 


' ■  -W.-M  .... , ,  . . .  . . — . . . . . m*''  . .  “  '*•>•>  "*•  . 


perhaps  best  be  appreciated  by  comparing  the  variation  of  the  incremental 
tine  delay  for  the  single  crossover  model  with  foroing  function  bandwidth 
shown  in  Fig.  7?  with  the  similar  plot,  Fig.  80,  nade  for  the  extended 
models.  The  major  differences  between  these  stem  from  the  absence  of  the 
correction  for  crossover  frequency  increase  with  a&  in  Fig.  80.  In  this 
respect  the  simple  crossover  model  results  are  superior  to  those  of  the 
extended  versions.  If  the  decreasing  slops  tendency  which  is  especially 
evident  for  the  Kc/(J<o  -  2)  controlled  element  is  to  be  taken  seriously, 
an  asymptotic  lower  level  for  Te  would  be  indicated.  This  would  appear 
to  have  a  value  qpmewhere  near  0.1  sec  or  Issb. 

The  variation  of  the  incremental  low  frequency  phase  lag  measured  by  a 
with  forcing  function  bandwidth  is  shown  in  Fig.  81  .  Two  limiting  cases 
appear.  For  the  highly  constrained  Ke/(ja>-  2)  condition,  a  is  substanti¬ 
ally  Independent  of  forcing  function  bandwidth.  At  the  other  extreme,  for 
Yc  ■  K0/(j«)2,  a  is  approximately  proportional  to  a*.  The  simple  integrator 
controlled  element  has  a  trend  of  a  versus  cuj.  which  ranges  between  the  two 
extremes.  The  low  frequency  Incremental  pbaae  lag  has  no  precedent  in  the 
existing  analytical-verbal  model.  However,  it  is  necessary  to  reoall  that 
such  terms  are  not  new.  They  have  appeared  before  and  ware  carefully  con¬ 
sidered  in  Ref.  34  before  being  thrown  out  as  major  off eats.  At  least  one 
of  these  "reasons"  (i.e.,  that  the  points  were  inconclusive  due  to  varia¬ 
bility  in  computational  problems)  is  no  longer  as  strong  as  it  once  was. 
Other  "reaeons"  auch  as  the  relatively  inconsequantlal  nature  of  the  low 
frequency  phase  lags  in  terms  of  closed-loop  characteristics  are  still 
pertinent  for  some  controlled  elements,  but  not  all. 

Both  the  variation  or  the  lack  thereof  of  the  low  frequency  incremental 
phase  lag  and  the  effective  time  delay  with  forcing  function  bandwidth  are 
major  experimental  facts  whioh  must  bt  explained  in  any  adequate  model  of 
the  human  operator.  This  has  been  accomplished  in  another  phase  of  the 
research  program  of  whioh  this  report  documents  a  portion,  although  it  la 
beyond  our  present  soope.  It  turns  out  that  essentially  all  the  variations 
noted  are  compatible  with  a  variable-gain  adaptive  model  of  the  neuro¬ 
muscular  system  which,  in  turn,  is  generally  compatible  with  the  known 
physiological  characteristics  of  the, same  system. 


Low  Frequency  Phase  Lag  Coefficient 

a  (sec1)  Incremental  Time  Delay  , 


Figure  80.  Dependence  of  Incremental  lime  Delay 
on  Forcing  Function  Bandwidth 
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Figure  M  .  V*  rations  of  Low  Frequency  Incremental  Phase  Lag 
with  Forcing  Function  Bandwidth 
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Examination  of  Fig.  44  and  its  comparison  with  the  pre-experiment 
analyst  a  shown  in  Figs.  6  sad  9  indicates  that  s  minima  font  for  the 
open- loop  describing  function  should  be 


KpKo^o  +  l)e~3tiyte 
Jm(jo>  -  l/T) 


Suitable  parameters  o&n  be  selected  in  this  form  such  that  the  amplitude 
ratio  is  reasonably  well  described  over  the  measurement  band,  while  the 
phase  is  accurately  characterized  only  in  the  region  of  crossover.  In 
fact,  some  of  these  data  absolutely  require  the  low  frequency  phase  lag 
correction  in  order  to  be  fitted  at  all  well,  even  at  frequencies  near 
crossover.  Consequently,  the  more  refined  eeoond  approximation  form, 


«^KC(TU»  * 

4<o(j«  -  l/T) 


will  be  used  at  the  outset. 

From  comparison  of  the  several  plots  it  is  noted  thnt  the  major  trends 

ar*!  «.  a*,  increases  with  l/T  (see  Fig.  4jb) 

b.  Phase  margin  is  generally  very  low  and  decreases  as 
l/T  increases  (see  Fig.  46b) 

i 

o.  The  low  frequency  phase  lag  Increases  with /an 
increase  in  l/T 

As  a  consequence  of  the  first  of  these  trends,  Ta  will  deorsaee  as  l/T 
increases,  up  to  the  point  where  control  is  lost. 

Just  as  with  the  Y0  *  K^/ ( Joi)2  case,  a  large  number  of  variations 
were  tried  in  fitting  the  low  frequency  data.  As  it  turned  out,  the 
lead  time  constant,  Tl,  can  be  considered  Independent  of  l/T,  with  a 
value  of  about  5  see .  Then  a  and  t*  will  carry  the  primary  variation 
with  controlled  element  change.  The  actual  curve  fits  are  given  in 
Figs.  82a  through  82d  and  tabulated  in  Tablo  XIII.  The  uee  of  a  alone  to 
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approximate  effect*  which  are  probably  due  to  a  low  frequency  lag-lead 
ha*  a  tendency  to  overestimate  the  low  frequency  phase  lag  at  the  very 
lowest  frequencies  of  measurement.  Shi*  is  especially  noticeable  in 
Figs.  88c  and  8Sd.  She  effect  can,  of  opurse,.  ht.ellBdJMted  by  using 
the  lag-lead  instead  of  a,  with  the  consequent,  addition  of  another  ... 
parameter. 

SABLE  XIII 

SUMMARY  OF  DESCRIBING  FUNCTION  CONSTANTS 
FOR  Y0  -  Xo/jtu(jto  *•  l/T),  a*  .  1  .5 


As  can  beat  be  appreciated  by  examining  Fig.  8?,  the  increase  in  low 
frequency  phase  lag  and  decrease  in  effective  time  delay  behave  In  essen¬ 
tially  a  parallel  faahion  as  the  controlled  element  divergence  is  increased. 
la  fan-c,  it  la  evident  that  the  conflict  between  the  low  frequency  phase 
lag  and  the  high  fn^uency  phr  e  decrease  accompanying  rc  decreases  is 
the  significant  factor  in  the  pilot's  ability  to  control  diverging  con¬ 
trolled  elements.  Here,  beoaue*  of  the  intimate  oonneotion  between  Te 
and  a,  the  deareaee  oan  only  be  effective  as  long  a*  its  effect  la  not 
overpowered  by  the  incremental  low  frequency  phase  lag.  The  source  of 
the  close  tie  between  and  a  again  involves  the  details  of  the  neuro¬ 
muscular  subsystem,  and  h*noe  is  again  beyond  our  current  scope. 

S.  YKDOnXOR  MODlXi 

In  its  moat  refined  version,  the  extended  crossover  model  given  by 
%  1  f’OC) 
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Figure  83.  Variation  of  Effective  Time  Delay  and 
Low  Frequency  Z&oremental  Phase  lag  with 
1/T  for  Y0  -  Ko/jaUa}  ~  l/t) 
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is  sufficient  to  characterize  almost  ill  the  data  for  all  the  controlled 
elements  tested.  Exceptions  occur  at  very  low  and  very  high  frequencies. 
For  instance,  in  the  discussion  completed  just  above,  values  of  a  based 
on  the  lower  variability  mid-frequency  and  low  frequency  data  tend  to 
overestimate  the  phase  lag  at  the  very  lowest  frequencies.  Also,  for  the 
Y0  »  Ko/( Jcu  -  2)  data  the  very  lowest  and  vary  highest  amplitude  ratio 
points  are  not  compatible  with  curve  fits  using  the  extended  crossover 
model  alone.  These  deviations  between  data  and  model  offer  opportunities 
for  model  extension. 

A  more  subtle  defect  of  the  extended  crossover  model  of  Eq  TOO  is  the 
possible  confusion  resulting  when  this  form  is  used  at  frequencies  for 
which  it  was  never  intended.  Neither  of  the  exponential  term's  components 
can  be  extrapolated  beyond  the  measurement  band  without  the  possibility 
of  conceptual  error.  It  must  always  be  remembered  that  the  e”^a'^  term 
represents  lags  and  leads  having  break  points  below  or  In  the  lowest 
portion  of  the  measurement  bandwidth;  and  similarly  that  at  least  part 
of  the  includes  the  low  frequency  effects  of  leads  and  lags  having 

break  •noints  above  or  in  the  highest  regions  of  the  measurement  band. 

Thus,  extrapolation  of  Eq  100  either  to  very  low  frequencies  where  the 
ph_se  will  tend  to  minus  infinity  or  to  very  high  frequencies  where  the 
amplitude  ratio  approaches  a  constant  is  unwarranted. 

The  primary  intent  of  this  section  is  to  take  advantage  of  the 
additional  scope  offered  for  model  extension  by  the  data  which  are  not 
quite  characterized  by  the  refined  crossover  models  while  also  removing 
the  subtle  limiting  frequency  difficulties  of  the  extended  form.  To  do 
this  a  precision  model  which  contains  many  more  components  is  formulated. 
Moat  cf  those  additional  components  have  break  points  outside  the  measure¬ 
ment  bandwidth  which  must  be  estimated  by  their  ">ften  small  effects  within 
this  band.  Consequently,  only  the  lowest  variability  data  available  are 
appropriate  for  use  as  an  initial  data  base  in  model  elaboration.  For 
this  reason  the  precision  model  developments  here  will  use  only  the 
exceptionally  low  variability  data  obtained  for  the  highly  constraining 
Kq/ ( jar  —  2)  controlled  element. 
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Sven  with  these  rather  remarkable  data,  the  development  of  a  precision 
model  con  be  only  academic  unless  a  good  deal  of  additional  information 
is  considered  in  the  model  construction  process.  When. this  is  doae,  how¬ 
ever,  the  precision  model  evolved  is  compatible  not  only  with  the  data  of 
the  current  program,  but  with  the  other  findings  introduced  into  the 
structure  as  well. 

The  model  development  process  can  best  be  described  as  a  sequence  of 
operations,  each  based  on  one  or  more  fundamental  considerations.  Each 
of  the  steps  included  are  detailed  below. 

1 .  As  a  minimum  to  replace  the  a  and  to  better  fit  the  lowest 
frequency  amplitude  ratio  data,  a  lag-lead  form  is  indicated 
wherein  the  lead  is  Just  within  the  measurement  band  and  the  lag  occurs 
at  much  lower  frequencies.  Careful  consideration  of  the  | Yp | /Kp  data  at 
low  frequencies  for  the  three  forcing  function  oases  indicates  that  these 
are  nearly  identical.  This  feature  is  also  reflected  by  the  values  previ¬ 
ously  found  for  a,  which  are  essentially  the  eame  for  all  three  conditions. 
Therefore,  the  lead-lag  added  will  be  the  same  tor  all  three  forcing  func¬ 
tion  conditions.  This  very  low  frequency  lag-lead  has  no  amplitude  ratio 
and  negligible  phase  effoct  at  the  higher  frequencies,  so  the  break  points 
can  be  determined  independently  of  other  parameters.  The  values  selected 
are: 


2.  Data  and  interpretations  from  this  and  other  experimental  programs 
have  increased  our  knowledge  of  the  dynamics  of  the  neuromuscular 
system.  As  already  noted  several  times,  the  details  of  these  are  beyond 
the  present  scope,  but  two  findings  related  bo  the  precision  model  devel¬ 
opment  task  are: 

a.  A  fairly  complete  description  of  the  high  frequency 
neuromuscular  system  dynamics  requires  a  third-order 
system,  i.e., 

_ 1 _ 

+ ’)  [(c§)2  +  4?  +  1 

All  these  terms  are  shown  in  Eq  1 .  Their  presence  there  was  based 
largely  on  high  frequency  neuromuscular  system  dynamics  as  revealed 
by  the  dynamic  portion  of  step  response  data.  More  fundamentally, 
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thi  third- order  form  constitutes  a  minimal  analytical  description 
of  the  muscle  impulse  response. 

b.  The  effective  neuromuscular  system  high  frequency  lag-, 


+  %  Oos) 

tends  to  vary  in  the  same  fashion  as  a.  In  this  instance,  for 
Y0  m  Kq/(Jcj  -  2),  a  is  aonstant  with  ay..  (Therefore,  Tjf  will  also 
be  essentially  constant. 

3.  For  the  data  being  considered,  the  components  of  the  effective 
time  delay,  re,  will  be 


Te  4  T  +  ®Ni  +  ajj-  “  TL  4  (t+Tn)-Tl  '  (TO?) 

The  pure  time  delay,  t,  is  due  to  latencies  in  the  visual  process,  delays 
due  to  neural  conduction  and  coding,  etc.  Accordingly,  its  value  is  not 
likely  to  change  with  forcing  function  bandwidth.  Thus,  for  these 
Yq  m  Kq/(J<u  -  2)  data  the  term  T  +  Tjj  will  be  independent  of  the  forcing 
uituaction  bandwidth.  In  Eq  103  this  leaves  the  lead,  T^,  which  is  generally 
adjustable,  as  the  likely  source  for  the  observed  changes  in  re  with  oij. . 

if.  Minimum  values  of  l/T j,  can  be  found  by  considering  that: 

a.  The  mid-frequency  regions  for  the  |Yp|/Kp  data  are  the 
i*ame  for ’all  ay,  thereby  indicating  that  1  /Tl  influences 

only  the  two  highest  frequency  amplitude  ratio  data  points. 

b.  Because  the  neuromuscular  system  dynamics  do  not  change 
with  cm.  in  this  particular  case,  the  departure  of  the 

highest  frequency  amplitude  ratio  point  (o^  ■  13.3  rad/ sec)  from 
the  partial  amplitude  ratio  curve  fit  (which  includes  the  lead, 

I’L,  but  excludes  the  neuromuscular  system  lags)  should  be  the  same 
for  all  three  forcing  function. conditions.  With  only  a  small 
amount  of  cut-and-txy,  lead  values  can  be  found  which  are  compati¬ 
ble  with  these  factors.  These  are: 


<4 


Tl, 


max 


Te  +  Tl 


0.156 

0.125 

0.104 


1.5 

2.5 
4.0 


13.0 

9.0 

7.8 


0.077 
0.111 
’  0.128 


0.233 

0.236 

0.232 


Note  that  t*  +  Tr.  Is  essentially  a  constant,  as  le  required  if 
t  +  %  la  to  be  Invariant  with  c^. 

5,  At.  -thie  5ai%t  the  third- order.  neuromacu^r  aad  pure  tins  4S3*y ; 
effects  must  be  estimated  from  their  influence  on  about  fiv*  or  . 
■lx  phase  points  and  two  amplitude  satio  points.  A  unique  determination 
of  these  parameters  is  not  to  be  expected,  but  good  fits  which  are  gen¬ 
erally  compatible  with  associated  data  oan  be  found.  As  a  first  step, 
the  pure  time  delay,  t,  can  be  estimated  from  physiological  data,  its 
components  include:* 


Flash  stimulus,  shortest  path  to  cortex .  0.04  to  0.07  sec 

Conduction  within  the  cortex  to  convergence 
center  of  visual,  propri oceptive,  and  kines¬ 
thetic  information  to  formulate  motor  output  0.005  sec 


Cortex  to  spinal  cord . .  0.01  to  0.015  sec 

Spinal  cord  to  periphery  (via  alpha 

motoneuron  a) . . .  0.005  to  0.01  sec 

0.06  to  0.10  see 

For  the  preeent  data  a  value  of  0.09  sec  shall  be  used.  This  Is  selected 
because  it  Is  near  the  upper  end  of  the  range  given  above  and,  inciden¬ 
tally,  because  it  worked  somewhat  better  than  other  values  within  this 
range  in  the  aut-and-try  fitting  efforts. 


6.  The  residual  amplitude  ratio  and  phase  left  over  when  the  terms 
thus  far  considered  are  removed  from  the  data  are  shown  in 
Fig.  84.  These  data  are  to  be  fitted  with  * 


_ 1 _ 
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After  considerable  aut-and-try,  representative  numerical  values  were 
determined  to  be 


1ST 


These  correlate  reasonably  well  with  representative  third-order  system 
approximations  which  characterize  the  so-called  dynamic  portion  of  the 
step  response  (see  Ref.  34)  for  Ye  •  Kc. 


*Personal  communication— -Dr.  Oeorge  Moore,  Aset.  Prof,  of  Physi¬ 
ology,  UCLA  Medioal  School. 
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■me  final  forma  found  to  fit  the  three  tc^  caeea  Tor  Y0  ■  Kq/Uo)  -  2) 

are 


(104) 


for  <a^  m  1 .5,  2.5,  and  4.0,  respectively.  The  over-all  merit  of  this 
detailed  precision  model  can  be  judged  from  Figs.  85a  through  85c,  which 
present  the  data  and  the  precision  model  forms  of  Eq  104.  The  g«nerel 
adequacy  of  these  curve  fits  io  exceptionally  good  evidence  for  the  exten¬ 
sion  of  the  general  trodel  form  shown  in  Eq  1  to  that  given  by 


(105) 


In  this  more  general  model  the  gain  end  equalization  terms  are  identical 
to  those  previously  used,  whereas  the  neuromuscular  nystem  characteristics 
have  been  expanded  by  the  addition  of  the  low  frequency  lag-lead. 

An  interesting  aspect  of  the  precision  model  is  Its  reduction  of  pure 
time  delay  terms  to  minimum  or  near-minimum  values.  All  the  high  frequency 
leads  and  lags  whioh  were  previously  lumped  into  low  frequency  approxima¬ 
tions  represented  by  e”^®  have  now  been  removed,  leaving  only  the  bare 
bones  of  a  pure  time  delay  whioh  is  compatible  with  that  expeoted  on 
physiological  grounds. 
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!to  this  section  some  of  the  adjustments  revealed  by  the  describing 
function  date  will  be  rationalised  in  terms  of  the  simple  noxmliESd 
mean- squared-error  models  derived  in  Chapter  XXX. 

1 .  Regression 

The  first  topic  is  that  of  regression.  ThiB  phenomenon  tree  first 
noticed in  Elklnd’s  rectangular  forcing  function  spectra  data.  Moat  of 
these  had  open-loop  describing  function  characteristics  reasonably  com¬ 
patible  with  those  given  in  Eq  1 0,  so  the  normalized  mean- squared- error 
model  of  Eq  35  forms  an  appropriate  analytical  framework  with  which  to 
"explain"  these  data  and  the  regression  phenomenon. 

In  Table  XIV  are  presented  data  based  on  Elkind'a  original  results, 
Ref.  13.  The  normalised  mean- squared  error  is  only  that  component  which 
has  frequencies  within  the  forcing  function  bandwidth.  This  la  a  olose 
approximation  to  ef/of .  Because  the  rectangular  forcing  function  spectra 
used  by  Elklnd  had  no  power  in  the  region  of  crossover,  the  crossover 
frequency  for  the  lower  forcing  function  bandwidths  was  not  direotly  meas¬ 
urable.  However,  extrapolation  of  the  available  low  frequency  describing 


TABLE  XIV 

NORMALIZED  MEAN-SQUARED-ERROR  DATA 
DERIVED  FROM  ELKIND'B  Ye  -  1  EXPERIMENTS 


FORCING 

FUNCTION 


(rad/sea) 


tI 

a* 

«c 

(rad/ sec) 

*o“i 

4 

*1 

0.0081 

7.6 

1.57 

0.21 

0.0058 

0.011 

7.6 

1.57 

0.31 

0.013 

0.041 

7.6 

1.37 

0.32 

O.036 

0.098 

7.6 

1-57 

0.83 

0.092 

0-31 

6.3 

1  -30 

1  .24 

0.30 

1.16 

2.0 

0.4l4 

2.07 

R.16 

R.24 

R.40 

R.64 

R.96 

R1.6 


1 .0 

1.5 

2.5 
4.0 
6.0 


10.0 
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function  data  to  the  crossover  region  reveals  that  the  <%  vas  probably 
the  same  for  all  of  these  oases  (this  extrapolation  vas,  in  faot,  the 
primary  Justification  for  Adjustment  Rule  5b) .  Bikini’s  B-6  forcing 
funotion  s^otra  oan  be  used  to  estimate  a  likely  crossover  frequency 
for  the  lover  aaj.  data.  This  is  the  basis  of  the  7.6  rad/sec  shown  for 
R.16  through  R.64  (see  Tig.  29).  The  T0  corresponding  to  this  oroesover 
frequency,  when  it  is  oonsidered  as  an  is  0.21  sec.  The  rest  of 
the  data  in  Table  XXV  follows  directly.  Figure  86  shows  the  data  from 
Table  XXV  superimposed  on  the  normalized  mean-squared  error  versus 
normalized  crossover  frequency  families.  This  plot  actually  involves 
three  quantities,  i.e.,  ef/of>  and.  ^o^d  fa118  constant  ?0  is  a 

normalizing  parameter) .  The  actual  variation  being  illustrated  1b  e^/c^ 
versus  (normalized  ),  with  ay.  (normalized)  as  a  parameter.  The  data, 
are  entered  on  the  plot  using  the  first  two  as  abscissa  and  ordinate.  The 
arrow  leading  from  each  data  point  has  its  head  at  the  appropriate 
Thus,  .the  R.96  data,  for  instance,  show  a  normalized  mean-squared  error  of 
O.JI  versus  normalized  crossover  of  1 .30  with  the  arrowhead  touching  a 
family  7^  ■  1  .24.  If  the  arrow  lengths  were  zero,  the  experimental  data 
would  coincide  precisely  with  all  three  parameters  of  the  plot. 

The  approximately  oonotant  oroesover,  low  forcing  function  bandwidth 
data  and  the  R1.6  large  #n  bandwidth  point  all  correspond  quite  closely 
to  the  theoretical.  The  lover  Ud  data  are  sIbo  compatible  with  mini¬ 
mization  of  the  mean-squared  error,  whereas  the  high  014  point  (R1.6) 
Illustrates  the  regression.  The  R.96  data  point  departs  a  good  deal 
from  the  idealization.  This  too  is  as  it  should  be,  for  the  open- loop 
describing  funotion  here  is  transitlonnry  and  departs  considerably  from 
the  simple  -20  db/decade  constant- slope  amplitude  ratio.  This  data  point 
could  have  been  left  off  the  plot  for  this  reason;  Instead,  it  is  entered 
to  lndloate  how  closely  the  theoretical  ay.  <  0.8a>c  condition  is  likely 
to  apply. 

In  the  present  program  only  one  condition  resulted  in  an  regression 
This  was  for  Yc  «  ^/(jo)2,  <14  ■  4.0  (see  Pig.  45a).  The  idealized  mean- 
squared-error  plot  of  Pig.  11  (or  Pig.  86)  is,  unfortunately,  not  as 
quantitatively  appropriate  for  the  Y0  ■  K0/(jco)2  data  as  it  is  for  the 
Y0  ■  K0  or  Kq/Jo  cases  beoauoe  the  open-loop  phase  differs  drastically 


from  that  presumed  in  the  theoretical  calculations .  Still,  the  actual 
data  hear  a  qualitative  resemblance  to  the  theoretical  behavior,  I.e., 
is  essentially  constant  for  lew  forcing  function  bandwidths,  and  the 
crossover-  regresses  when  %  approaches  this  value. 

t.  fevforaa&oe  Measures 

In  Section  0  of  Chapter  Ill  estimates  of  the  trend  of  normalized  mean* 
squared  error  with  ujj.  were  made  presuming  the  crossover  model  of  Eq  10. 
These  estimates  can  now  be  compered  with  the  actual  results.  The  predic¬ 
tions  were  made  in  terms  of  the  mean- aqua red- error  component  due  to  the 
forcing  function  alone,  so  only  that  portion  will  be  considered  here. 

Included  in  the  data  for  Elkind's  yq  ■  1  experiments  sunorarized  in 
Table  XIV  are  predicted  mean-squared  errors  based  on  the  one-third  law 
(Eq  37).  These  predictions  are  compared  with  the  actual  data  in  Fig.  87. 
The  correlation  between  the  theory  and  experiment  indicated  in  this  figure 
is  quite  satisfactory. 

Similar  data  based  on  the  Y0  ■  Kq/Jw  aeries  reported  here  are  presented 
in  Fig.  88.  Three  types  of  correlations  are  shown.  The  first,  similar 
to  that  already  shown  for  Elkind'e  data,  Is  with  the  one-third  law.  For 
the  second  correlation,  denoted  in  Fig.  38  as  sif/afj  oontinuous*  the  14oal“ 
ized  data  of  Fig.  11  are  used  for  the  predicted  values.  The  "continuous" 
subscript  here  ref ere  to  the  continuous  nature  of  the  foroing  function's 
rectangular  power  spectral  density  used  in  the  calculations  leading  to 
Fig.  11.  In  the  actual  experimental  series  the  foroing  function  vae  com¬ 
posed  of  discrete  sinusoids  rather  than  continuous  functions.  This  makes 
a  small  difference  In  the  normalized  mean. squared-error  prediction,  as  is 
shown  by  the  third  correlation  given  In  Fig.  88,  i.e.,  the  aotual  normal¬ 
ized  mean-squared  error  compared  with  *f/a?],iiBorete‘  Estimates  made  on 
all  three  bases  compare  favorably  with  the  aotual  results  aohleved. 

For  the  final  comparison  of  estimated  versus  actual  performance 
measures,  the  Y0  ■  Kc/( Jai)2  situations  shall  be  examined.  Figure  69 
presents  correlations  for  this  oass  which  parallel  those  given  In  Fig.  88 
for  the  single-integrator  controlled  element.  Here  the  correlation  is  by 
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Figure  89.  Comparison  of  Actual  Mean-Squared  Error 
with  Estimated  Values  Based  on  Various  Assumptions, 

Yc  ■  Kc/fjo,)2 
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□o  means  as  good  as  that  for  the  two  previous  oases  examined.  The  primary 
reason  for  this  is  the  departure  of  the  actual  open- loop  phase  angle  from 
that  of  the  simple  crossover  model.  To  Illustrate  the  magnitude  of  this 
effect)  calculations  of  a|/ °f]  discrete  W8re  ^  the  phase  for 

only  the  highest  frequency  point  within  the  large  amplitude  portion  of 
the  forcing  function  was  changed  from  its  actual  value  to  one  compatible 
with  the  simple  crossover  model.  This  solitary  modification  would  change 
the  ®f/of]4iecrete  points  to  the  locations  indicated  by  the  arrowheads  In 
Fig.  89.  These  points  would  be  considered  very  good  correlation  indeed, 
eo  they  imply  that  even  the  one-third  law  could  be  used  for  the  Yc  ■  Kc/(Ja))2 
case  with  simply-applied  special  corrections. 

The  data  correlations  given  above  show  how  well  the  one-third  law  works 
even  though  the  actual  forcing  function  and  YpYc  differ  considerably  from 
the  idealized  model  on  which  the  law  is  based.  This  Implies  that  lower 
limit  values  of  steady-state  performance,  i.e.,  normalized  mean-Bquared 
error  due  to  forcing  function,  can  be  estimated  using  extremely  simple 
idealizations  for  conditions  where  the  crossover  model  is  a  reasonable 
approximation.  The  main  source  of  difficulty  in  such  an  estimation  pro¬ 
cedure  is  in  the  determination  of  the  crossover  frequency.  For  this  task, 
the  techniques  illustrated  in  Chapter  III  can  be  used,  with  the  Te  values 
modified  to  be  consonant  with  the  new  data  of  this  report.  The  discussion 
above  has  bypassed  any  problems  of  a*;  estimation  by  considering  trends 
only,  and  by  using  experimentally  determined  valueB  of  o\,. 

3.  tg  Trajeotoris*  tad  Phase  Marcia  Adjustment 

The  connection  of  performance  measures  and  mean- squared-error  minimiza¬ 
tion  with  crossover  frequency  1b  but  half  the  story  of  the  crossover  model 
adjustment.  The  other  parameter  involved  is  the  adjustment  of  te  or  its 
exact  equivalent  for  the  crossover  model,  the  phase  margin.  An  apprecia¬ 
tion  for  the  possible  rationale  behind  this  adjustment  can  be  obtained  by 
considering  Fig.  90 ,  which  is  a  crosBplot  of  Fig.  11 .  In  this  figure  the 
absaissa  is  the  effective  time  delay,  Te,  normalized  by  the  mean  crossover 
frequency,  The  families  shown  are  for  constant  forcing  function  band¬ 
width,  u>£,  also  normalized  by  cq- .  Thus,  the  variation  being  illustrated 
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here  1b  ef/of  versus  re  (normalized)  with  cm  (normalized)  as  a  parameter. 

On  this  basis  Fig.  90  can  be  thought  of  as  a  seotion  of  a  hyptrsurface 
for  whioh  Fig.  11  is  another  section.  It  will  be  observed  that 
for  small  values  of  normalized  mean- squared  error  the  cm/ty;  family  lines 
are  nearly  horizontal,  indicating  that  normalized  mean-squared  error  is 
nearly  independent  of  tb.  This  region  corresponds  to  values  for 

which  the  one-third  law  is  an  adequate  approximation  to  the  normalized 
mean-squared  error.  As  forcing  function  bandwidth  is  increased,  the 
families  tilt  toward  the  left  and  the  mean-squared  error  becomes  a  much 
stronger  function  of  the  time  delay.  The  minimum  mean-squared  error, 
which  always  corresponds  to  Te  *  0,  becomes  a  much  sharper  minimum  as 
cui/ZEfc  increases. 

For  the  human  controller  with  random  forcing  functions,  Te  can  never 
become  zero,  but  the  Te  data  do  Bhow  distinct  reductions  as  functions  of 
cuj,.  The  nature  of  tblB  adjustment  in  terms  of  the  theoretloal  mean-squared 
error  versus  a^,Te  plot  is  illustrated  by  the  data  superimposed  on  the  theo¬ 
retical  grldwork  of  Fig.  90  •  The  data  used  for  this  purpose  are  for 
Yc  ■  Kc/j<o  and  Kc/(  Jo>  -  2).  These  caseB  were  selected  because  they  corre¬ 
spond  most  closely  to  the  simple  crossover  model.  The  actual  data  points 
are  entered  on  the  plot  using  the  phase  margin  and  .  The  nature  of  a 
linear  phase  margin  variation  with  cu^  is,  in  these  coordinates,  indicated 
by  the  qijfl  ■  1  >3(^/3^)  curve  shown.  A  phase  margin  adjustment  of  this 
kind  places  sufficient  emphasis  on  phase  margin  Increase  with  cpi  to  assure 
that  the  Bystem  with  time  delay  has  a  mean-squared  error  reasonably  close 
to  the  physically  unrealizable  absolute  minimum  for  t  ■  0.  The  human  oper¬ 
ator  approximates  this  adjustment  for  £  0.3,  but  tends  to  fall  off 

somewhat  for  larger  values,  although  phase  margin  still  Increases. 

An  alternative  view  to  the  same  adjustments  .noted  above  can  be  developed 
using  the  sketch  of  Fig.  91 .  This  figure  shows  the  open-loop  and  the  closed- 
loop  error/input  transfer  functions  for  a  feedback  system  based  on  the  simple 
crossover  model.  The  maximum  amplitude  ratio  of  the  error/input  transfer 
function  is  very  nearly  the  departure  at  the  break  point,  which  coincides 
with  the  crossover  frequency,  ofc.  As  noted  in  the  plot,  this  is  a  function 
of  the  phase  margin  alone,  and  for  small  phase  margins  the  peak  magnitude 


181 


Figure  91.  Sketch  of  Closed-Loop  Error/Forcing  Function 
Transfer  Function  for  System  Eased  on  Simple  Crossover  Model 

iu  approximately  equal  to  the  phase  margin  in  db.  Nov  consider  a  rectangu¬ 
lar  forcing  function  spectrum  with  bandwidth  ay..  The  resulting  system 
mean-squared  error  will  be  proportional  to  the  area  under  the  crosshatched 
portion  of  the  O^e  amplitude  ratio  curve  vhen  this  is  transformed  to  linear 
units  and  squared.  When  ay.  is' small  compared  with  tug,  the  actual  and 
asymptotic  | Oi* | curves  are  nearly  coincident.  Under  these  circumstances 
any  peaking  effect  at  toe  <*ue  to  small  phase  margins  has  very  very  little 
effect  on  the  mean-squared  error.  This  is  the  region  wherein  the  one-third 
law  applies.  As  ay  gets  closer  to  ao^,  the  effect  of  the  amplitude  ratio 
departure  from  the  asymptotes  becomes  more  significant.  Thus,  the  appro¬ 
priate  adjustment  as  cy/o^  becomes  larger  is  for  the  phase  margin  to  be 
increased,  thereby  decreasing  the  peak  and  the  positive  departure  of  the 
actual  curve  from  the  asymptote,  and  hence  the  rms  error. 
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It  will  be  recalled  from  Chapter  II  that  the  primary  purposes  of  the 
experimental  series  were  the  validation  of  the  existing  analytical-verbal 
model  and  the  extension  of  this  model  in  aacuracy  and  detail.  The  experi¬ 
mental  data  analyses  and  interpretations  reported  here  have  accomplished 
these  purposes  and,  in  addition,  have  revealed  many  other  facets  of  human 
pilot  dynamics.  The  total  effort  can  be  conveniently  summarized  under 
four  major  headings:  (l)  The  Data  in  Qeneralj  (2)  Status  of  the  Exist¬ 
ing  Analytical-Verbal  Describing  Function  Model)  (3)  Extensions  to  the 
Analytical-Verbal  Describing  Function  Model;  (4)  Status  of  Bemnant  Data. 
Conclusions  of  this  program  relating  to  these  topics  are  summarized  in 
outline  and  tabular  form  below.  The  ranges  of  known  validity  are  restricted, 
of  course,  by  the  limits  inherent  in  the  data  and  analyses  presented  earlier. 

A.  TEE  SATA  IE  CDUBRAIi 

I.  General  nature.  The  data  accumulated  have  shown  consistent, 
repeatable  results  which  can  be  fitted  into  and/or  extend  the  general 
context  of  past  results  and  theories  derived  therefrom. 

fl.  Tie-In.  The  data  taken  for  Yc  -  1  are  close  enough  to  ElMnd's 
to  constitute  our  tie-in  with  past  data  and  to  enable  hie  data  to  ba 
considered  ae  a  subset  of  oure. 

3«  Quasi-Linearity.  Conventional  quaei-linear  constant-coefficient 
describing  funotion  models  can  be  evolved  to  characterize  the  data.  The 
primary  variables  which  fix  the  conditions  for  quasi-linearity  in  this 
study  are  controlled  element  and  forcing  funotion. 

4.  oi  Independence.  There  is  no  evidence  of  nonlinear  dependence  of 
the  describing  function  measurement  on  forcing  function  amplitude. 

J.  leleotivs  Variability.  Human  controllers  exhibit  a  remarkable 
capaoity  to  suppress  all  sources  of  variability  and  to  operate  with  high 
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repeatability  and  intersubject  uniformity  in  frequency  region*  where 
necessary,  e.g.,  at  crossover  or  over  wider  frequency  bands  in  the 
control  of  conditionally  stable  systems.  In  other  frequency  bands  where 
variations  In  pilot  dynamics  do  not  materially  affeot  closed- loop  per¬ 
formance,  there  is  muah  more  intra-  and  interpilot  variability,  and 
pilots  exhibit  individual  styles. 

6.  Normality  of  Describing  Function  Data.  Both  amplitude  ratio  and 
phase  data  appear  to  be  distributed  in  a  Gaussian  manner. 

7-  Time-Varying  Behavior.  In  general,  pilot  dynamic  characteristics 
are  reasonably  stationary  over  the  measurement  run  lengths  as  indicated 
by  the  typical  p  ranges  over  the  main  power  portion  of  the  forcing  func¬ 
tion  bandwidth  given  below  for  uy.  »  2.5. 


*c 

p  Range 

Kc/Joj 

Ke/(jffl  -  2) 
K0/(Jo>)2 

4 

0.97  -  0.99 

O.99  - 

0.85  -  0.98 

a.  Selective  variability.  There  is  more  time  variation  In 
frequency  rangee  wherein  variability  is  relatively  unimportant  to  closed- 
loop  performance  than  over  more  crucial  frequency  ranges. 

b.  Task  ooqplsxlty,  o^.  Stations rity  decreases  as  ay.  increases 
(e.g.,  a  change  in  p  from  O.98  to  0.95  as  ay.  changes  from  1 .5  to  4.0  for 
Yq  ■  Kq/Joj)  .  This  effect  is  not  so  pronounced  when  the  controlled  element 
ia  made  more  constraining,  e.g.,  Y0  ■  K0/(jui  -  2). 

0.  Task  complexity,  Yq.  As  the  controlled  element  dynamics 
become  more  difficult  to  control,  the  time  variation  of  the  pilot's 
dynamics  increases  substantially.  In  this  sense,  p  could  conceivably  be 
used  as  a  measure  of  task  difficulty. 
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B.  KATUB  or  m  mKEM  (CZRGA  1$60) 
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1 .  In  almost  all  cases  tbs  general  describing  function  form 


Kpe“J<tfr«(TLjCD  +  1) 
(Tjjco  +  1 )  (TNJ<n  +  1 ) 


(106) 


can  be  suitably  adjusted  to  provide  a  satisfactory  description  for  cross¬ 
over  region  characteristics. 


2  There  is  striking  evidence  that  for  a  variety  of  forcing  func¬ 
tions  and  controlled  elements  the  slope  of  YpYc  at  crossover  is 
-20  db/decade.  Under  these  conditions  the  model  of  Eq  106,  when  taken 
in  company  with  the  controlled  element  dynamics,  reduces  to  the  especially 
simple  form 


where 


>*c 

.  o^e-J^e 

(107) 

jo) 

“c 

i  0*jo  +  0.l803j_ 

T0 

■  t0  -  0.08o>l 

(108) 

T0 

■  «/2“b0 

Values  of  T0  for  conditions  other  than  those  tested  can  be  estimated  using 
the  data  given  below  and  interpolation  based  on  the  slope  of  | Yc I ^  over 
the  likely  crossover  region. 


Yc 

dl*ddb 

d  In  m 

(db/dec 

OJ-OJq 

jade) 

To 

(sec) 

Kc 

0 

0.33 

Kc/jo) 

-20 

0.36 

c> 

K> 

*4o 

0.52 

18? 
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Tor  Yc  -  Ko/ja>,  ^/(jto)2,  Ko/(ja>-  2),  and  Ko/jo>(Jo>  -  l/T), 
which  constitute  critical  experiments  for  the  olroa  i960 
analytical-verbal  model;  the  data  obtained  support  «  priori  predictions 
based  on  this  model  as  to  Yp  form.  Pra- experiment  estimates  of  cross* 
over  frequencies  were  properly  ordered  and  reasonably  accurate,  although 
generally  somewhat  higher  than  the  experimentally  determined  values. 

e|  performance  measure  trends  based  on  the  crossover  model  for 
Ye  -  Kc  and  Kq/Jcu  are  supported  by  the  experimental  data. 

The  essentially  constant  crossover  frequency  for  uy  <  0.8a*. 
and  the  phenomenon  of  a*,  regression  for  ay,  >  0.8a*  can  be 
explained  using  the  crossover  model  of  Eq  10?  and  the  presumption  that 
the  pilot  tends  to  minimize  the  mean- squared  tracking  error. 

g  The  adjustment  rules  for  the  existing  analytical-verbal  model 
are  generally  adequate,  although  Adjustment  Rules  4b  (7a rams tar 
Adjustment  -  system  phase  margin)  and  5b  and  5°  (flfe-Oi  Independence  and 
0*  Regression)  require  minor  modifications  in  the  nature  of  refinements 
(see  below).  In  terms  of  the  adjustment  rules  the  major  impact  of  the 
experimental  program  has  been  to  validate  with  concrete  data  those 
"rules"  which  had  previously  been  based  primarily  on  extrapolations  and 
rationalization.  In  this  sense  all  of  the  rules  given  in  Chapter  II, 
with  the  modifications  noted  below,  are  now  solidly  based  on  experiment. 

a.  mnmicM  to  m  muriaurmuL  dhoti aim  ru»orirar  model 
1 .  Adjustment  Rules 

a.  Adjustment  Rule  4b  should  now  be  stated  as: 

"System  phase  margin,  qXp  1b  directly  proportional  to  ay,, 
the  forolng  function  bandwidth,  for  values  of  a*,  lees  than 
about  2.0  rad/sec.  The  strong  effect  of  forcing  function 
bandwidth  on  the  phase  margin  is  associated  with  the  varia¬ 
tion  of  T[}  with  the  same  task  variable." 

h.  A  new  adjustment  rule,  4c,  should  be  added: 

"Equalization  time  constants  or  Tj,  when  form  selec¬ 
tion  requires  l/Tj,  or  l/Tj  «  a*,  will  be  adjusted  such  that 
low  frequenoy  response  will  be  essentially  insensitive  to 
slight  changes  in  Tj,  or  Ti  (ay.  «  a*), 
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A  .  IkA  JUfl f  Pljl  Ap  ShOliii  HO*?  wiia^n^A^  ah  « 

3b.  Independence:  System  oroaeover  frequency  depends 

.  «inly  Slightly  ok  forcing  function  bandwidth  for  ay<  0.6  0bo 

(ttfe0  is  that  value  of  a^  adopted  for  ay.  «  a*j). 

5o.  cq,  Regression:  When  ay.  nears  or  beoomes  greater  than 
0.8a^o  ,  the  crossover  frequency  regresses  to  values 
much  lower  than  u*j0. 

2.  Low  frequency  These  Lags.  A  lagging  phase  angle  at  low  frequen¬ 
cies  appears  in  all  the  data.  The  general  low  frequency  describing 
function  form  of  Eq  106  is  not  suitable  at  the  lower  frequencies.  To 
adequately  characterize  the  data  it  must  be  modified  to 


P  "  (TiJoj  +  1 )  (%!»  +  1 ) 


(109) 


The  parameter  a  added  here  is  a  "high  frequency"  approximation  to  low 
frequency  leads  and  lags  occurring  below  the  measurement  bandwidth.  In 
higher  order  approximations  to  the  low  frequency  phase  lag,  a  lag-lead 
largely  outside  the  bandwidth  of  the  measurements  is  usually  adequate, 
a  iuareaaes  both  with  ay  and  with  controlled  element  order,  except  for 
the  constraining  controlled  element,  Yc  ■  Kq/CJoj  -  2).  Although  not 
described  in  this  report,  the  nature  of  a  and  its  variations  oan  be 
explained  by  properties  of  the  neuromuscular  system. 

3-  Input-Adaptive  Neuromuscular  Characteristics.  The  first-order 
lag,  Tn,  approximation  to  the  high  frequency  neuromuscular  system,  dynamics 
is  input  adaptive.  Its  decrease  is  ordinarily  tae  basis  for  the  varia¬ 
tion  of  Te  with  ay.  indicated  by  Eq  108  [the  data  for  Yc  *»  Kc/(jcu  -  2) 
are  an  exception  to  this  statement]  . 

4.  Higher  Order  Neuroausoular  lags.  The  previous  strong  evidence 
for  a  third-order  neuromuscular  lag  has  now  been  added  to  by  describing 
function  measurements.  These  lags  replace  essentially  all  of  the  residual 
previously  assigned  to  a  Te  or  to  a  r  and  Tjj  in  more  approximate  models. 
The  net  effective  lag,  Tn1  +  2{;n/u$j,  is  essentially  equal  to  the  first- 
order  approximation,  Tf],  and  is  variable  with  ay  in  a  similar  fashion. 
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J.  Residual  t.  The  third-order  neuromuscular  system  plus  a  pure 
time  delay  of  about  0,09  sec  are  adequate  descriptors  of  the  high  fre¬ 
quency  describing  function  data.  The  pure .time-:  delay  is  of  the  same  : 
order  of  magnitude  as  neural  delays  In  the  vleu&l  modality,  lie. , 
latencies  in  the  visual  process  plus  conduction  delays,  etc . 

6.  (tenesmi  Describing  Function  Fern.  The  most  general  describing 
function  model  form,  replacing  Eq  1 ,  is  given  by 

(110) 


[aw®+i\ 

- 

1 

“ 

) 

Larj  VTjjcn  +  lf 

|W>+i/ 

(Tn,  Jco+l) 

7rf+  !5a  jB+i 

\«*J/  +  o*j 

The  major  elements  of  complication  introduced  In  this  precision  model 
are  those  in  the  braces.  These  all  arise  from  the  neuromuscular  system, 
vhich  has  both  very  high  and  very  low  frequency  effects.  For  almost  all 
practical  cases  Eq  110  can  be  simplified  to  Eq  106  or  even  Eq  107 •  When 
the  controlled  element  is  such  as  to  make  the  syStom  conditionally  stable, 
Eq  109  should  be  used  for  safety's  sake,  Bince  the  low  frequency  phase 
lag  can  be  Important  for  this  kind  of  system. 

s.  KAMI  or  mourn  data 

1 .  Values  of  the  remnant  computed  at  the  forcing  function  frequen¬ 
cies  generally  fit  a  smooth  curve  through  values  measured  between 
and  above  forcing  function  frequencies.  This  Indicates  that  the  power 
spectral  density  of  the  remnant  is  generally  continuous  and  that  line 
spectra  indicating  periodicities  are  absent. 

g  At  very  low  frequencies  the  remnant  data  for  a  wide  variety  of 
controlled  elements  coalesce  best  when  all  the  remnant  is 
reflected  to  the  pilot's  input. 

.  Remnant  increases  with  controlled  element  gain,  with  forcing 
function  bandwidth,  and  with  control  order.  For  extreme  con¬ 
trolled  element  forms  such  as  Y0  -  Kq/Jc^Joi  -  1 .5)  the  remnant  increases 
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greatly,  primarily  because  of  the  pilot's  time-varying  behavior  induced 
by  his  attempts  to  retain  control  ovgr  this  drastically  unstable  controlled 
element . 

Some  evidence  for  pulflihg  behavlor  in  control  of  second- order 
controlled  elements  is  present  from  output  amplitude  distribu¬ 
tions.  These  indicate  a  tendency  for  the  pilot's  output  to  be  pulse  areas 
roughly  proportional  to  the  stimulus  amplitude. 

.  Careful  examination  of  the  output  power  spectral  density  Indicated 
no  evidence  for  periodic  sampling  or  significant  nonlinear  behavior 

g  Partly  by  process  of  elimination  and  partly  by  direct  evidence, 
it  appears  that  the  major  source  of  remnant  is  nonstationary 
pilot  behavior,  i.e.,  time-varying  components  in  the  effective  time  delay 
and  gain.  For  the  second-order  controlled  elements  the  pulsing  nature  of 
the  pilot's  output  contributes  an  additional  remnant  source. 
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